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VULCAN - A DEPLETION PROGRAM FOR USE WITH INDEPENDENT 
ONE-DIMENSIONAL SPATIAL CALCULATIONS 
by John L. Anderson, Jr. 

Lewis Research Center 

SUMMARY 

A digital computer depletion program, VULCAN, which treats the time behavior of 
nuclide concentrations in a nuclear reactor, is described. This depletion program is 
suitable for use with any multigroup one -dimensional spatial solution; that is, it is 
independent of the actual generating source of cross sections or fluxes. This flexibility 
allows any space -energy calculational sequence to be extended, without alteration, to 
include a depletion stage. 

The depletion calculation is performed for each mesh interval, and an averaging 
process is used to reduce the number of regions for which new atom densities are 
required. The depletion equations for several classes of nuclides are shown and cast 
into forms solvable by a computer. Data input instructions, a sample problem, and a 
FORTRAN IV listing of the program are included. VULCAN, written for an IBM 7094, 
relies entirely on fast memory computer storage (32K). 


INTRODUCTION 

In the face of various economic and operational restrictions the life expectancy of a 
nuclear reactor must be considered as a part of its design. This life expectancy or 
lifetime is essentially the length of time a reactor core or fuel element may operate 
before reactivity or structural considerations force a shutdown. For example, the 
structural consideration may be the governing one because of gaseous fission products 
and consequent pressure buildup. 

In the interest of economics the degree of fuel utilization within a spent element or 
core should be as large as possible. Budgetary requirements and, for commercial 
power reactors, subsequent consumer costs make necessary the availability of a cost 
per unit time figure which depends on the fuel utilization. In another vein, space ap- 
plications of nuclear reactors, for which maintenance is rather inconvenient, demand 


lifetime and performance foreknowledge . 

The lifetime determination based on reactivity considerations may be made by con- 
sidering the time behavior of fissionable and parasitic nuclide concentrations. The 
behavior of individual nuclides is governed by depletion (productive and parasitic ab- 
sorption, decay and physical losses) and by buildup (decay and fission sources). The 
differential equations describing the various isotopic populations are typically solved by 
programs for digital computers; generically, these programs are called depletion 
programs. 

Depletion programs may emphasize any of several facets: speed, intervention con- 
venience, tape handling convenience, and calculational detail. Usually, the programs 
include the spatial (generally diffusion theory) calculations and cross section libraries in 
one all-encompassing computer program. These programs may treat reactors whose 
symmetry allows their calculation to be performed in one (refs. 1 and 2), two (refs. 3 
and 4), or three dimensions. They allow several depletion cycles to be performed with- 
out any other input or intervention. However, such codes typically have energy group 
and diffusion theory limitations because of insufficient computer memory capacity. These 
limitations may be overcome by auxiliary storage on peripheral devices but at the ex- 
pense of computer execution time. A somewhat more comprehensive review of the 
literature and programs may be found in reference 5. 

This depletion program, VULCAN, may be used with any one -dimensional space - 
energy solution. The cross sections and fluxes (mesh interval averages) that VULCAN 
requires may be from any source (constant, analytical, or even experimental) in any 
format. Consequently VULCAN is independent of the actual generating source of cross 
sections or fluxes. This flexibility allows any space -energy calculational sequence to be 
extended without alteration to include a depletion stage. A particular advantage result- 
ing from the independent status of VULCAN is the availability of considerably more fast 
memory computer storage since the multigroup energy and spatial solutions to 
Boltzmann’s equation have been obtained from other programs. This additional storage 
enables pointwise (at each discrete mesh interval of the spatial calculation) rather than 
zonal (average over several mesh intervals) nuclide concentrations to be followed. A 
representative 7 energy group-14 mesh interval-13 isotope-1 time step problem re- 
quired 0. 15 minute of IBM 7094-H computer execution time. 

The VULCAN program is comprised of about 1200 source statements and relies 
entirely on fast memory computer storage. However, because no peripheral storage is 
used, the following restrictions are imposed on this version of VULCAN: 
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Nuclides 

o 

VI 

Energy groups 

o 

VI 

Mesh intervals 

<20 

Zones 

<190 

Materials 

=2190 


Some of the restrictions may be relaxed but at the expense of tightening others (see 
the section VULCAN DESCRIPTION). VULCAN is written in FORTRAN IV language for 
use on an IBM 7094 with 32K memory. 


PREPARATION OF DEPLETION EQUATIONS 


General Solutions and Approximations 


The general equation governing the concentration of a particular isotope N as a 
function of space r, time t, and energy E is 

NF y*oo r°° 

— = 2 r^N.(r , t) J $(F, E, t)aj(E) dE + X p N p (r , t) + N r (r , t) J $(r , E, t)o^ y (E) dE 


- X N N(r,t)-N(r ,t) / 3> (r , E, t)o^(E) dE - £ N N(r,t) (l) 


where 


<h(r,E,t) 


0 ^, 0 * ,o* 

f’ n, y’ a 


N k 

NF 


neutron flux as function of space, energy, and time, neutrons/(b) (sec); 
neutrons/(10“ 28 m 2 ) (sec) 

microscopic fission, (n,y), and absorption cross sections for isotope k, 
b; 10- 28 m 2 

decay constant for isotope k ? sec”'*’ 

on q 

atom concentration of isotope k, atoms/(b) (cm); atoms/(10 m ) 
number of fissionable isotopes considered 


The first term on the right side of the equation is the source of isotope N from the 
fission yield of each fissionable isotope i considered. The second term is the 
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source from beta decay of the precursor nuclide N^. The third term is the source from 
radiative capture (n, y reaction) by N r . The fourth term is the loss resulting from beta 
decay of the isotope N itself. The fifth term is the loss from depletion or burnup of N. 
The sixth term represents losses of isotope N such as expulsion of material fragments 
or diffusion of gaseous constituents or products (£ N is the loss rate). 

There are two simplifications to be considered before this generalized equation may 
be used in a depletion scheme. The first is the reduction from the continuous variables 
r , E, and t to their discrete counterparts. In particular, the GAM-GATHER (refs. 6 
and 7) multigroup energy structure and the TDSN (ref. 8) spatial interval structure may 
establish the discrete E and r. There remains only the discrete t to assign, and 
this is done as part of the VULCAN input. 

The second simplification begins with the realization that the form of equation (1) and 

its consequent solution will depend on the nuclear behavior of the isotope. Contrast, for 

235 

example, the form of the equation for uranium 235 (U ) which is fissionable and may be 
formed by radiative capture in U and the form for xenon 135 (Xe 1 ) which is a fission 
product and is both created and depleted by beta decay. 

Now, before enumerating the forms of the depletion equation to be considered, it is 
convenient to discuss some methods for their solution. There are two assumptions in 
general use which facilitate the solutions: (1) constant flux or (2) constant power over a 
time interval, each of which will be expressed mathematically hereinafter. As the time 
interval decreases, the two approximations approach each other. However, the constant 
power approximation (CPA) seems to be appropriate since the usually understood objec- 
tive of the reactors considered is that of providing a constant power output. Since most 

232 233 

of the power in a reactor is produced by any of six isotopes (thorium 232 (Th ), U , 
U 235 , U 238 , plutonium 239 (Pu 238 ), and Pu 2 ^),the CPA need only be applied to 
equations which dictate their behavior. The equations for all other isotopes may be 
treated by the constant flux approximation and then solved analytically. 

At this stage it is appropriate to cast equation (1) in its discrete form, that is, 

$(?, E, t)o£ (E) dE 

becomes 



NG 

V (t)cr N 
/Lj jm v ' x, 

j J 

for any isotope N, mesh interval (spatial location) m, cross section cr , and time t 
and where NG is the total number of energy groups. Thus, 
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dN m (t) 
dt 


NF NG 


NG 


'm 


S r^N. (t) V $. (t)oj + A N m (t) + N (t) y (t)a^ 
i im w Z-/ ' i* p pm 7 rm w l~j jm w n , y • 

i j J j J 


NG 


- W‘> - N m<« Z - Vm® (2 > 


In the CPA, N. m (t) 


NG 


L, ^jm°f. * s a extant f°r any mesh interval m where i cor- 


J 


J 


responds to the six isotopes mentioned previously. In the constant flux assumption, 
is considered constant over the time interval, and N m (t) is thus free to be combined 


jm 


with dN m (t) leading to an analytic solution 


„ and 

Jin 


X i 


Of course, in either case, 

are implicitly constant over the energy j and space m intervals as is required by the 
respective codes. For the finite difference solution to the CPA, the following two sub- 
stitutions must be made: 


dN(t) N t+At ~ N t 
dt At 


and 


N(t) - 


N t+At + N t 


Formulation of Specific Equations 

It now becomes necessary to determine the isotopes which will be explicitly con- 
sidered. A fixed base of twelve fissionable isotopes and five assorted fission products 
was established, and the basic depletion equations determined. These equations and two 
other general forms are considered adequate for any foreseeable reactor application. 

The twelve fissionable isotopes used as a base are the constituent nuclides in two 
chains Th^^ — and |j238 __ p u 242. fi ve ka se isotopes are Xe*®®, 

samarium 149 (Sm^), iodine 135 (1*^), promethium 149 (Pm 1 ^*), and a fission 

135 149 

product aggregate; I and Pm are considered merely because they are precursors 
135 149 

of Xe and Sm , which are major parasitic absorbers. Four additional categories 
of isotopes may also be treated. The first of these categories is labeled "Special 
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Materials” and may include additional fission products and aggregates and neutron 
capture products. These isotopes are formed from fission and/or capture and decay 
(see the section INPUT INSTRUCTIONS). The analytic solution is used and so is appli- 
cable to rapidly saturating nuclides. The second is headed "Burnable Poisons” and 
treats isotopes depleted by absorption and those which are not continuously formed (e. g. , 

boron 10). The third category, called "Repeated Fissionable Nuclides" allows any 

235 

fissionable isotope to be included more than once. If, for example, U occurs in two 
locations in a reactor where the flux spectra differ, the microscopic cross sections may 
be sufficiently different to make separate depletion desirable. The fourth category 
is headed "Nondeple table Nuclides" and it allows nondeple table isotopes such as hydrogen, 
oxygen, or aluminum to be carried through the VULCAN calculation unaltered. 

The following two sections contain the specific depletion equations considered. Their 
solutions, in a form suitable for computer analysis, are shown along with some exemplary, 
but condensed, derivations. In the equations that follow, any (/<!> implies a type of sum- 
mation convention, that is, 


NG 



1 


u 3>. - 

X j 3 


cr* 

x 


Fissionable nuclides . - The following two simplified heavy element chains provided 
the twelve basic fissionable nuclides: 


Th 232 


T t236 , . . 

U chain: 


U' 


238 


— Pu 242 chain: 



- U 


236 




* 




* 


- Pu 


242 
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(n, y) process followed by ]3 decay 
/3 decay 
(n, y) reaction 

fission (all isotopes fission but virtually all power is produced by the six indicated) 

The two chains appear to obey identical reaction mechanisms. To the extent of the con- 
siderations in this depletion analysis, this is the case, with but one exception: Pu 2 ^* 
undergoes /3 decay with a nonnegligible decay whereas U , for these depletion pur- 
poses, is stable. Consequently, seven equations suffice to treat the twelve fissionable 
nuclides. 

In the following equations the superscripts represent the VULCAN internal identifica- 
tion (VID) number (see the section Nuclide List) for the first chain considered 
(Th 232 _ u 228 ). ^he isotope in the second chain (u^S — Pu 2 ^ 2 ) to which each equation 
applies is also indicated. The equation designation (a) will denote the differential equa- 
tion and (b) will denote the finite difference form used in the program. 


where 

\ 

I 

S 


Th 232 (VID = 1) or U 238 (VID = 7) 


— = - nM* 

dt a 


(3a) 


N h4t- N t_ K At +N t\l e 

a 


At 


1 i/2 - At 

i+At = N t( 

^2 + At , 


(3b) 


Pa 233 (VID = 2) or Np 239 (VID = 8) 


dN‘ 

dt 


:=n 1 (^-^)®-n 2 (a 2 + ^®) 


(4a) 
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in N 2 (x 2 + cr 2 ®) - NMo* - <7 f l 


X 2 + n 2 ®)at 


»* - <At- N t j,_ a *p 

2 x 2 + n 2 *1 


- (x 2 + CT 2 $j At |>+ N 2 exp - (x 2 + CT 2 $j At 


U 233 (VID = 3) or Pu 23y (VID = 9) 


3 N 3 - N 3 

= - N 3 ct 3 $ + N 2 X 2 1 

dt a At 


N Lt * + n 2 x 2 

2 / a 


N t 3 ( 2 - "a* at ) + Rat + N ?V at 


2 + At 

a 


If VID = 2 (VID = 8) is not included, then the following equation is used: 


N ? ( 2 - a a* At ) + K+At + "t) \°a ' °f)* At 


2 + cr At 

a 


U 234 (VID = 4) or Pu Z * U (VID = 10) 


dN _ ,4 , -m-2/ 2 ~2\ A (rfi _ rr 2 


-*X* + KK-°lF + lt * 


. . -cr 4 <I> At [/ o 9 

N t+at = N t e a + ( N t+at + N t 


♦to**** UL- 


-cr^<I> At 

a 


ii i «■■■■« III i 


ll III l 



u 235 (VID = 5) 


n: 


t+At 


_ N 5 a^$ + N 4 (a 4 - cr 4 j<h 

= N t ( 2 - j > &t) + (< At + N?)(°a - 4Y At 


2 + <£® At 

a 


Pu 241 (VID = 11) 


dN 


11 


dt 


= - N 44 ^a 44 $ + X 44 ) + N 40 ^ 40 - cr 40 )<& 


N 


11 

t+At 


N t n [ 2 - ( x11 + °r *H + fet + N t'% 


40 - o - 40 |$ At 

3, I 


2 +[ X 11 + a 44 <h)At 

a 


U 236 (VID = 6) or Pu 242 (VID = 12) 


£[N_ = _ N 6 ct®$ + N 5 
dt a 


(°l - * f 5 > 


N 


t+At 


. N? ( 2 - "a* at V ( N ta At + N tK 


cr^ - of ) $ At 
3. I 


2 + cr $ At 

3 


(7a) 


(7b) 


(8a) 


(8b) 


(9a) 


(9b) 


1 i qp; 1 4Q 1 4 q 

Poisons: Explicit fission pro duct s (Xe 1 , I , Sm ^_and Pm ) . - Each of 
these explicit fission products is treated under a separate heading. The analytic solu- 
tions enable the rapidly saturating fission products to be followed within a single time 
step. For each of the four products an equilibrium concentration equation is given, and 
in the case of Xe 435 the buildup after shutdown is treated. 

i op- 135 

Xe . The equation governing the time behavior of the Xe nuclide concentration 
X(t) is 
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(10a) 


dX(t) 

dt 


NF NG 


- v (t > + E E 4 ‘ s 


i=l 


J 


J 


A x X(t) - <£*X(t) - ^X(t) 


where 


Tji 1 

Sj = N <j£ 



is the average macroscopic fission cross section for isotope i, energy group j ; £ x 
represents the removal rate for xenon other than by a? and X ; and kr is the 

a X 1 

corresponding quantity for iodine. 

In the following derivations, in the interest of brevity, the time dependence will not 
always be shown and the summation over all energy groups will be implied. The Xe 00 
concentration at time t + At is then 



NF NG 

in which * x = \+ £ x , x j = x j - ^ = E E , and I t and X^ represent 

i j J 

concentrations of 1^^ and Xe*^ at the start of the time interval. In the event I*®® is 

not considered, equation (10b) is altered internally by letting X. = 0 and = T^ + rf 

135 1 

The equilibrium or steady-state Xe concentration is given by 


X ss = 


x'i + r* 
I ss 

X x + 


(10c) 


10 



with F* defined as in the preceding sentence and I gg representing the I 1 ** 5 equilibrium 
concentration. 

135 

The equation for Xe concentration at a time At after shutdown may be obtained 
from equation (10b) by letting <f>, r 1 , and r* equal zero. Then, I t and X t become 
*sd and X sdV the respective concentrations at shutdown. The time after shutdown of 
maximum Xe**^ concentration is 



(10d) 


135 135 

I : The isotope I does not have a decay precursor, thus its equation differs 
slightly from that of Xe 1 0 (eq. (10a)). Another loss mechanism £j (e.g. , diffusion or 
expulsion) may be incorporated into Aj, that is, Aj = Aj + £j. The 1*35 nu clide concentra- 
tion obeys the equation 


MU - A*I(t) - Ktjo 1 + 
dt 1 a 


£ 


i=l 


r? 

i 


NG 


j=l 




(11a) 


In finite difference form, equation (11a) becomes 


t+ At 


= I t exp [-(^+<^*)At] + —I^ 


X T + 


a 


i - 



(lib) 


The steady -state concentration is given by 


I 


ss 


r 1 

* I 
X I + CT a* 


(11c) 


149 149 

Pm : The isotope Pm has the same creation and destruction mechanisms as 
135 

I except that its effect on reactivity is small enough so that no additional loss 
mechanisms are considered. Thus the finite difference equation is 



and the steady state concentration is given by 


P ss = 


rP 


V 0 ?* 


(12b) 


149 135 

Sm : This isotope is similar to Xe except that it is stable and no extra loss 


mechanisms are considered necessary. The equations for Sm*^ are 

NF NG 

_?(t) + 

dt p i t-J i 4 - “'-'-a 


V (t> - Z *f S 4 * - 

i j=l 3 


(13a) 


and 


3 t+At 


1i - [-( x p + - exp (- °s* &t ) 




- X - a P 4> 
a p a 


r S (x + ct P $)+ T P X r / o \1 / a \ 

7" P"\/ S \ P L 1 “ 6XP (" °a* At )J + S t 6XP (‘ °2* At ) (13b) 

( X p + <%*) (of 1 ) V ' 


149 

In the event Pm is not considered, equation (13b) is altered internally by letting 
A = 0 and r® = r® + r P . The steady state concentration is 

r 

X p + r s 

S = P s s 

SS Q 

of® 


(13c) 


Special materials: These materials, which include fission products and aggregates, 
obey the following depletion equation: 


n: 




t+At 




- exp o^ P <I> Atj + N® P exp o^ P $ Atj 

(14) 


SP 

where 6 = 1 if material SP is a capture product (e.g. , (n,y)) of the preceding 

material (SP - 1) and 6 =0 if it is not (see the section INPUT INSTRUCTIONS). 


op 

For the fission product aggregates (VID =13), 6 = 0 
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Burnable poisons: The concentrations of the nuclides considered as burnable 
poisons obey the simple analytic form 

N f + At = N f exp(- <£<& At) (15) 

Depletion cycle: These depletion equations may now be incorporated into a complete 
depletion program. Such a depletion program, coupled with the energy and space codes 
(e. g. , GAM-GATHER and TDSN, refs. 6 to 8) enables a depletion cycle to be performed. 
A typical depletion cycle proceeds according to the flow diagram in figure 1. 



figure 1. - VULCAN (low diagram. 

Before entering the depletion program, two time intervals must be chosen. The shortest 
time interval Tj will be that between renormalizations of the flux level (within VULCAN) 
to a given power. The other time interval Tg will be that between recalculations to the 
spatial flux distribution. VULCAN will treat the Tg/Tj normalizations sequentially as 
a single problem. 


SAMPLE PROBLEM 

This sample problem serves as a check on the solution of the depletion equations 
insofar as the original solution is concerned and as a device for displaying certain options 
of the program. To illustrate the flexibility of VULCAN, a problem previously solved by 
the depletion program CANDLE (ref. 2) but found in reference 9, was chosen as the 
sample problem. The problem has 2 energy groups and 14 mesh intervals comprising 
three regions in a slab geometry, and the depletion is over a single time step. Although 
the input to VULCAN illustrates the use of several options, only the output pertinent to 
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the sample problem in reference 9 is given. The input information (e. g. , cross sections 
and fluxes) listed in reference 9 was supplied to VULCAN. The first flux normalization 
factor in the VULCAN output (0. 999738 E-09) indicates the degree to which this input was 
duplicated. For exact duplication of the variables controlling the power, this factor 
should be 1. 0 E-09, with the exponent indicating the transition of flux level as used in the 
VULCAN input to neutrons per barn per second. The optional cross section input format 
used illustrates that particular option. Subroutine TABLE, as listed, contains current 
data; however, the following parameters assumed these values in the sample problem: 
THERNU(5) = 2. 46, THERNU(9) = 2. 88, THERNU(ll) = 3. 00, FIWATT = 3. 15xl0 10 , 
DECAY(14) = 2. llxlO" 5 , DECAY(16) = 2.88xl0' 5 , and DECAY(17) = 3.85X10" 6 . 

A comparative list of atom densities at the end of the depletion interval is given in 
table I. In reference 9, however, only region-averaged atom densities were printed out. 
Although VULCAN does not normally yield atom densities averaged over the original 
regions, an appropriate value of the averaging criterion e (see the section Subroutine 
OUTPUT) will yield these. 


TABLE I. - COMPARATIVE ATOM DENSITIES 


Nuclide 

Material 1 


CANDLE 
(ref. 9) 

VULCAN 

^35 

0. 138761-3 

0. 138763-3 

1/236 

.500194-6 

. 500687-6 

Pm 149 

.261273-7 

.261213-7 

jl35 

.297741-7 

.297672-7 

Sm 149 

.505137-8 

.505096-8 

Xe 135 

. 122126-8 

.122129-8 

FP2 

1/238 

Pu 239 

Pu 240 

Pu 244 

.258698-5 

.258238-5 










Material 2 


CANDLE 
(ref. 9) 

0.485625-4 
.937653-7 
.563627-8 
. 642296-8 
.136152-8 
.451111-9 
.558072-6 
.677028-2 
.655909-6 
.349750-8 
.227263-10 


VULCAN 

0.485632-4 


.561020-8 
.639333-8 
. 135575-8 
. 449453-9 
.554921-6 
.677051-2 
. 428466-6 
. 225741-8 
. 145501-10 


239 240 241 

The discrepancy in Pu , Pu , and Pu may be attributed primarily to the 

238 239 

absence of a resonance correction to the smooth cross sections of U and Pu ; 
VULCAN assumes the cross sections directly reflect any resonance structure. 
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Input 


l 

VULCAN TEST 

PKOBLEM 

FROM WAPO- 


2 

1 

18 



6 

14 

2 



1 

2 

4 



1 

2 

1 


6495, 

.33 1.0 

1.0 


50.0 

.01 

. 1 

.5 



1 

23 

HYD 



2 

24 

OXY 



3 

5 

1U-235 



4 

7 

IU-238 



5 

6 

U-236 



6 

17 

PM 149 



7 

16 

1-135 



8 

15 

SMI 49 



9 

14 

X E 1 3 5 



10 

13 

F P 2 



11 

21 15 

2-MD25 



12 

19 

810 1 



13 

20 

810 2 



14 

9 

1PU239 



15 

10 

PU240 



16 

11 

IPU241 



17 

18 

SPMAT 



18 

22 15 

X F I S2 



1 

2.54 

1 



1 

7.6? 

1 



1 

12.7 

2 



4 

35.9 

2 



1 

3H . 1 

2 



6 

68.1 

3 




13 1.0 

1.0 


1.0 


14 .003 

.003 


.003 


15 .0 

.0 


.0 


16 .06 

.06 


.06 


17 .014 

.014 


.014 


1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

( 36X » 2F 1 2 . 8 / J 
HYOKHGEN 


TM-95 (HAPD-TM-22U (CANDLE) KEF 9 
3 2 3 22 

0 

5 2 2 

0 0 1 


1.0 

.003 

.0 

.06 

.014 


.00444 


.2151 


OXYGEN 


.014 


U235 


U-238 


U-236 


PM— 149 


11.9416 H .8224 

418.2 872.32 

.36457 .12498 

1.79 0. 

.0 .0 

3.9 .0 

0. 


0 . 


0 



SM-149 


XE— 135 


FP2 


2ND25 


BIO 1 


BIO 2 


PU239 


PU240 


PU241 


SPMAT 


XFIS2 


0 . 

50000. 

0 . 

2700000. 

0 . 

65. 

12.1801 9.03777 

418.2 872.32 

60. 

2470. 

60. 

2470. 

0 . 0 . 

987. 1874.9 

0 . 0 . 

373. 0. 

0 . 0 . 

1080. 2415. 

60. 

2470. 

11.9416 

418.2 



7 

1 

2 

3 

12 

13 

17 


18 







.02738 


.01369 

.00014185 

.001 

.001 

.001 

.00014185 


4 

1 

2 

4 

11 



.023643 


.025464 

.006771 

.00004915 





2 

1 

2 





.051714 


.025857 






1.00431 


.922701 

.735115 

.528068 

.369309 

.252466 

.16132 

.103651 


.0646481 

.0314333 

.0152257 

.00725616 

.0032126 

.00090724 

.268374 


.263683 

.246155 

.202483 

.147527 

.104308 

.0801488 

.0830485 

.096636 

.085596 

.0556303 

.0315068 

.0154594 

.00458762 


2 







1.25 


1.25 







3 








5 

19 

20 





1.0 


1.0 






.99 


.80 






.99 


.80 







5 

1 

2 

3 

12 

13 


.02738 


.01369 

.00014185 

.001 

.001 




4 

1 

2 

4 

11 



.023643 


.025464 

.006771 

.00004915 





2 

1 

2 





.051714 


.025857 







3 








5 

19 

20 





1.0 


1.0 






.99 


.80 






.99 


.80 
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Output 


* MEMORY MAP * 


system 

FILE BLOCK ORIGIN 
FILES 1. UN l T 06 

2. UN I T05 

PRE-EXECUTION INITIALIZATION 
CALL ON OBJECT PROGRAM 
OBJECT PROGRAM 


00000 THRU 02717 
02 72 0 

(NO BUFF POOL ATT ACHED) 
(NO BUFF POOL ATT ACHEOI 
02750 
02757 

02764 THRU 76056 



DECK 

ORIGIN 

CONTROL SECTIONS 

(/NAME/=NON 0 LENGTH, ( LOC ) ’DEL ET ED, *=NOT REFERENCED) 

1. 

VULCAN 

02 7 64 

/SET1 / 02765 


/SF T2 

/ 23615 

/SET 3 

/ 37445 

/SET4 / 42235 

/SET5 / 42603 




/S=T6 / 43461 


/SET7 

/ 43471 

/SET 8 

/ 434 77 

/SET9 / 43507 

/SE T 10 / 43513 




/SETH / 43521 



51430 * 





2. 

TTABLE 

51 454 

/SET1 /( 02 7 6 5 1 


/ SE T2 

/(23615I 

/ S FT 3 

/ ( 37445 ) 

/SET4 / ( 42235 ) 

/SETS / { 42603 ) 




TABLE 51736 








3. 

TXENON 

51764 

/SET1 / ( 02 7 65 1 


/ SFT2 

/ (23615) 

/SET 3 

/ (37445) 

/SET4 /( 42235) 

/SFT5 / I 42603 ) 




/St T6 / ( 434 61 ) 


/SET7 

/ (43471 ) 

/SET 8 

/ ( 43477) 

/S ET 10 /( 43513 ) 

XENON 52470 

4. 

POUT 

52525 

/ SE T 1 / ( 02 7 65 ) 


/ S ET2 

/( 23615 ) 

/SET 3 

/( 37445 ) 

/SET4 / ( 42 2 3 5) 

/SETS / ( 42603 ) 




/SET6 / ( 43461 ) 


/SFT0 

/ (43477 ) 

OUTPUT 

54603 



5. 

FADE 

54623 

/ SE T l /( 02 765 1 


/ SET2 

/ ( 2 36 1 5 ) 

/ SET 3 

/ (37445) 

/SET4 /( 42235) 

/SETS /( 42603 ) 




/Sc T6 / (43 461 ) 


/SET7 

/ ( 4 347 1 ) 

/SETS 

/ (43477) 

/SET9 /( 43507 ) 

/SET10 /( 43513) 




FISEQ1 56026 


F I SEQ2 

56031 





6. 

A8FLUX 

56 034 

/ SE T 1 /( 02 76 5 1 


/ SE T2 

/ ( 2 36 15 ) 

/ S ET 3 

/ ( 37445 ) 

/SET4 / ( 42235 ) 

/SFT5 /( 42603 1 




/SET6 / (43461) 


/SE T9 

/ (43507 ) 

ABSPHI 

5 6464 



7. 

SEL SH 

56502 

/ SE T 1 / ( 02 765) 


/ SE T2 

/ (23615 ) 

/SET 3 

/ ( 37445 ) 

/ SE T4 /( 42235) 

/SET5 /( 42603 ) 




/ SE T 1 1 / ( 4 352 1 ) 


S SF 

56736 





8. 

BURN T 

56 7 64 

/ SE T 1 / ( 02 765 ) 


/ SET2 

/ (23615 ) 

/ S ET 3 

/ 1 37445 ) 

/SE 14 / ( 4 2235) 

/SE T 5 / ( 42 603 ) 




/ SE T6 / ( 43461 ) 


/ SE T7 

/ ( 4 3 4 7 1 ) 

/ S FT 8 

/ ( 43477) 

/S FT 9 / ( 43507 ) 

/ SE T 1 0 / ( 435 13) 




/SETH / ( 43521 ) 


BURNUP 

65766 





9. 

.L X CON 

66 006 

. L X STR 66006 


. LXSTP 

66011 

. LXOUT 

66033 * 

.LXERR 66042 * 

.LXCAL 66045 




.LXRTN 66045 

* 

I 0 E X I T 

66045 

. DBCLS 

66225 * 

.LXARG 66314 • 

.LO 66323 




• C L SE 66331 

* 

. LFBL 

66332 • 

. LUNB 

66333 * 

. OFQUT 66334 


10. 

-LXSL 

66340 

.LXSEL 66340 


. LXSL! 

66341 

. LXTST 

66344 

.LXOVL 66404 

.LXMOD 66446 




.LXIN3 66472 


. LXOI S 

66475 * 

.LXFLG 

66476 • 

. L TCH 66477 


1 1. 

.FPTRP 

66 5 05 

.FFPT. 66505 

♦ 

. FXEM. 

66507 

. FXEM 

66507 * 

.ETAG 66736 * 

.FMCRT 67052 




.:SOTO 67065 

* 

. FXOUT 

67071 * 

. FXARG 

67103 * 

.EXIT. 67105 * 

EXIT 67105 




.3VCHK 67112 

♦ 

OVFLOW 

67113 * 

SYS ONE 

67140 

.NOP 67141 

.BLANK 67154 

12. 

.ERAS. 

67 171 

E.l 67171 


E.2 

5717? 

E.3 

67173 

E . 4 67174 


13. 

.XCC. 

67 175 

CC.l 67175 


CC.2 

67176 

CC.3 

67177 

CC.4 67200 


14. 

FCNV 

67201 

• FC ON . 67201 


. ODPRE 

67211 * 

.FCNV . 

67217 

.ENDFS 67230 

.CNVSW 67232 




• F D XI 67236 


• FDX2 

67237 

. DBC 

67241 

. DBC 1 4 67304 * 

. STOP 4 67307 




. OB C 20 67323 


. OBCIO 

67337 

. 0BC99 

67352 * 

. DDSW 67375 * 

.FI XS W 67404 




-OXPSE 67420 

* 

. DOZE T 

67525 • 

. IC10 

67566 * 

. STOP J 67603 * 

• FC OUT 67613 




.FCAR3 67621 

• 

. ALCOD 

67623 * 

. ANPT 

67634 

. ON P T 67651 

.LNTP 67725 




.AOUT 67764 


. OOUT 

67776 

.LOUT 

70025 

.GOUT 70035 

. TGOUT 70044 




.FLT 70144 


. FXFL1 

70301 

• FXO 

70305 

. FXFL 2 70311 

*F XFL 3 70315 




. I NTS 70321 


. TOPAC 

70337 

.W IDT H 

70343 

•FPACK 70350 

.TEST 70351 




. DO SF 70402 

* 

. KOUNT 

70403 

.LIST 

70406 

.DONE 70415 

.OUTBF 70551 




.lHAR 71003 


. DE XP 

71012 * 

.TEN 

71016 • 

. FBD8F 71020 

.DATUM 71026 




.WORD 71041 


. MOD 

71042 

. PCX 

71043 

.FEXP 71044 

.DIG 71045 

15. 

FIOH 

71060 

.FIOH. 71060 


.DO 10 

71105 * 

• DDLET 

71243 * 

.DCPT1 71446 * 

. F XP T 71644 




.FFIL. 71675 


. DDF I N 

71713 • 

. FRT N. 

71722 

. DDR TN 71726 * 


16. 

FIOS. 

72030 

..FIOS 72 03 0 


. FCLS 

72052 • 

.. FIOC 

72077 

. . FSEL 72120 

..FBIK 72137 




. . F TC K 72141 


. I OFF 

72225 * 

. .X EM 

72233 * 

• . FCHK 72250 * 


i t • 

FR Du. 

tdi 53 

. .FRDD 72353 








18. 

FHRD. 

72400 

• • F WR D 72 4 0 0 








19. 

FBCD. 

72426 

..FBCD 72426 


. . FRCW 

72430 

• . FPCB 

72441 



20. 

UN IT 06 

72471 

. . UN06 724 71 








21. 

UN 0 5 

72472 

. UN05. 72472 








22. 

UN 0 6 

72473 

. UN 06. 72473 








23. 

FPUN 

72 4 74 

.FPUN. 724 74 








24. 

FLOG 

72531 

ALOGIO 72531 

* 

A LOG 

72532 





25. 

FXPF 

72 706 

EXP 72706 








26. 

• 10 E. 

73011 









27. 

. 1 0 E 56 

73023 









28. 

.RWDOE 

73031 









29. 

•BCREA 

73051 

BCREAD 73051 


BCREAD 

(73051 ) 





30. 

. BCR WO 

73131 

. . B C RD 73131 


. . BCWD 

73133 * 

. . BRD0 

73145 



31. 

FXP2 

73203 

• XP 2. -732 03 








32. 

PISTUG 

73 2 76 

/JOLO / 73277 


PI STUG 

(73276) 





33. 

PLOTXY 

74102 

/JOLO /( 73277) 


PLOTXY 

(74102) 





34. 

KHAR 

75761 

KHAR (75761) 








35. 

•LOGE. 

76033 

.LOGER 76033 








36. 

.tlex. 

76 041 

.TLEXP 76041 








37. 

•XFXP . 

76045 

.NGDEF 76045 


. XEXP 

76051 * 

.ZUOEF 

760 54 



UNUS 

ED CORE 


7605 7 

THRU 

77777 






BEGIN EXECUTION 



VULCAN fEST PROBLEM FROM WAPD-TM-95 l WAPD-TM-22 1 ) (CANDLE! REF 9 


NG 

2 


NFG NISOT NMAT 

110 3 


NT INC NFR ACT KCELL 
2 3 22 


N20NR NRM 

6 14 


M3R0UT 

2 


NUF 

0 


NSPMAT 

1 


KGEO 

1 


NBPOI 

NFIS 

NYDNUC 

NONDPL 

NFIRPT 

2 

4 

5 

2 

2 


KFLUX 

KSSF 

KHAI N 

KF AST 

KFORM 

2 

1 

0 

0 

1 


POWER 

IN 

WATTS 

. 64953E 04 


FIRST 
0 ISTANCE 
(CM) 
1-00000 


SECOND 

DISTANCE 

(CM) 

1-00000 


TIME 

XENON 

IODINE 


NC RE ME NT 

RE MOV AL 

REMOVAL 

DIFFER 

(HOURS) 

PROBABILITY 

PROBABIL ITY 


50.00 

-0. 

-0. 

-0. 


SEQUENT I AL 

VULCAN 

NRPT 

ISOTOPE 

ID 

ID 


LABEL 

1 

23 

-0 

HYO 

2 

24 

-0 

OXY 

3 

5 Y 

-0 

U-235 

4 

7 Y 

- 0 

U-238 

5 

6 

'0 

U-236 

6 

17 

-0 

PM149 

7 

16 

-0 

1-135 

8 

15 

-0 

SM149 

9 

14 

'0 

XE 135 

10 

13 

-0 

FP2 

11 

21 

15 

2N025 

12 

19 

- 0 

310 1 

13 

20 

*0 

810 2 

14 

9 Y 

-0 

PU239 

15 

10 

-0 

PU240 

16 

11 Y 

-0 

PU241 

17 

18 

- 0 

SPMAT 

18 

22 

15 

XF I $2 

NHIR( I ) 

R AR ( 

I) 

MI DRi 


1 2- 5400 
1 T.6200 
1 12.7000 
4 35.9000 
l 30.1000 
6 68.1000 


1 

1 

2 

2 

2 

3 


2 2 2 

3 3 3 3 3 


MESH VOLUMES 


0.254000E 01 0 * 508000E 01 0.508000E 01 0.580000E 01 0.580000E 01 0.580000E 01 0.580000E 01 0.220000E 01 

0.500000E 01 0-5 OOOOOE 01 0-500000E 01 0.500000E 01 0.500000E 01 0.500000E 01 0.681000F 02 


ZONE VOLUMES 


0 - 25400 OE 01 0-5C8000E 01 0.508000E 01 


0-232000 F 02 0.220000E 01 0.300000E 02 


FISSION 

YIELD 

YIELD 

YI ELO 

YI ELD 

YIELD 


PRODUCT 

FROM 5 
( U-235) 

FROM 7 
(U-230) 

FROM 9 
( PU239) 

FROM 11 
( PU241 ) 

FROM 

( 

17 

( PM 149) 

0.01400 

0.01400 

0. 01400 

0. 01400 


16 

(1-135) 

0.06000 

0. 06000 

0. 06000 

0.06000 


15 

( SMI 491 

0. 

0. 

0. 

0. 


14 

( X E 1 35 ) 

0.0030C 

0.00 300 

0.00300 

C. 00300 


13 

(FP2 ) 

I.COOOO 

1.00000 

I. 00000 

1. 00000 


18 

( SPMAT) 

0. 

0. 

0. 

0. 


GROUP NU FOR 

NU FOR 

NU FOR 

NU FOR NU 

FOR NU FCR 

NU FOR 


VID 1 

VIO 2 

V I 3 3 

VI 3 4 VID 

5 VID 6 

VID 7 

1 

0. 

0. 

0. 

0. 1. 

000 1.000 

1.000 

T 

2.600 

0. 

2.503 

0. 2. 

460 0. 

2.600 


YIELD YIELD 


YIELD YIELD 


NU FOR 
V ID 8 

0. 

0. 


NU FOR 
V t D 9 

1.000 

2.880 


NU FOR 
VID 10 

1.000 

0. 


NU FOR 
VID 11 

1.000 

3.000 


NU FOR 
VI) 12 

0- 

0- 


YI ELD 


18 



SEQUENTIAL ID = 1 VULCAN 10 - 23 

ABSORPTION BY GROUP I HIGH TO LOW) 
0.444000E-02 0.2L5100E 00 


SEQUENTIAL IC = 2 VULCAN 10 * 24 

ABSORPTION BY GROUP (HIGH TO LOW! 
0.140000E-01 -0. 


SEQUENTIAL ID = 3 VULCAN ID = 5 

ABSORPTION BY GROUP (HIGH TO LOW) 
0.119416E 02 0.41B200E 03 

NU *F I SS ION BY GROUP (HIGH TO LOWI 
0.882240E 01 0.872320E 03 

FISSION BY GROUP (HIGH TO LOWI 
0.8P2240E 01 0.3546C2E 03 


SEQUENTIAL ID = 4 VULCAN ID * 7 

ABSORPTION BY GROUP (HIGH TO LOWI 
0.364570E 00 0.179000E 01 

NU*F I SS ION BY GROUP (HIGH TO LOWI 
0.1 24980E 00 0. 

FISSION BY GROUP (HIGH TO LOWI 
0.1 2498 OE OC 0. 


SEQUENTIAL ID = 5 VULCAN 10 = 6 

ABSORPTION BY GROUP (HIGH TO LOWI 
0. 0.390000E 01 

NU *F I SS I ON BY GROUP (HIGH TO LOWI 

0 . 0 . 

FISSION BY GROUP (HIGH TO LOW) 

0 . 0 . 


SEQUENTIAL ID = 6 VULCAN ID = 17 

ABSORPTION BY GROUP (HIGH TO LOW) 

0 . 0 . 


SEQUENTIAL IC = 7 VULCAN ID = 16 

ABSORPTION BY GROUP (HIGH TO LOW) 

0 . 0 , 


SEQUENTIAL ID = 8 VULCAN ID = 15 

ABSORPTION BY GROUP (HIGH TO LOW) 
0. 0.5 OOOOOE 05 


SEQUENTIAL ID = 9 VULCAN ID = 14 

ABSORPTION BY GROUP (HIGH TO LOW) 
0. C • 270000E 07 


SEQUENTIAL ID = 10 VULCAN ID = 13 

ABSORPTION BY GROUP (HIGH TO LOW) 

0 . 0.6 50000E 02 


SEQUENTIAL ID = 11 VULCAN 10 = 21 

ABSORPTION BY GROUP (HIGH TO LOW) 
0.12180 LE 02 0.418200E 03 

NU*F I SS ION BY GROUP (HIGH TO LOW) 
0.903777E 01 0.872320E 03 

FISSION BY GROUP (HIGH TO LOW) 
0.903777E Cl 0.354602F 03 


SEQUENTIAL ID = 12 VULCAN ID - 19 

ABSORPTION BY GROUP (HIGH TO LOW) 

0 .-600000 E 02 0.247000E 04 


SEQUENTIAL ID = 13 VULCAN ID = 20 

ABSORPTION BY GROUP (HIGH TO LOW) 

0 .6 00000 E 02 0 • 247000E 04 


SEQUENTIAL ID * 14 VULCAN 10-9 

ABSORPTION BY GROUP I HIGH TO LOW) 
O. 0.987000E 03 

NU+FISSION BY GROUP (HIGH TO LOW) 

0. 0.187490E 04 

FISSION BY GROUP (HIGH TO LOW) 

0. C.651007E 03 



SEQUENTIAL ID = 15 VULCAN ID = 10 


ABSORPTION BY GROUP I HIGH TO LOW) 
0. 0.373000E 03 

NU*FISSI0N BY GROUP (HIGH TO LOW) 

0 . 0 . 

FISSION BY GROUP (HIGH TO LOW) 

0 . 0 . 


SEQUENTIAL 10 = 16 VULCAN ID = 11 

ABSORPTION BY GROUP (HIGH TO LOW) 
0. 0.108000E 04 

NU*F I SS ION BY GROUP (HIGH TO LOW) 

0- 0.241500E 04 

FISSION BY GROUP (HIGH TO LOW) 

0. 0.80500 OE 03 


SEQUENTIAL 10 = 17 VULCAN ID = 18 

ABSORPTION BY GROUP (HIGH TO LOW) 
0.600000E 02 0 • 247000E 04 


SEQUENTIAL ID * 18 VULCAN ID = 22 

ABSORPTION BY GROUP (HIGH TO LOW) 
0.1 1941 6E 02 0.418200E 03 

NU*F ISS ION BY GROUP (HIGH TO LOW) 

- 0 . - 0 . 

FISSION BY GROUP I HIGH TO LOW) 

- 0 . - 0 . 


VULCAN 

ID 

ISOTOPE 

LABEL 

MATERIAL 

1 

MATERI AL 

2 

23 

HYD 

0.27 380E-01 

0.23643E-01 

24 

OXY 

0. 13690E- 01 

0.2 5464E— 01 

5 

U- 235 

C. 14 185E-03 

0. 

7 

U— 238 

C. 

0. 67710E-02 

6 

U— 236 

0. 

0. 

17 

PM 149 

0. 

0. 

16 

I- 135 

0. 

0. 

15 

SH 149 

0. 

0. 

14 

X E 135 

0. 

0. 

13 

FP2 

0. 

0. 

21 

2ND25 

C. 

0.49150E-04 

19 

B1C 1 

C . 100006-02 

0. 

20 

BIO 2 

0.10000E-02 

0. 

9 

PU239 

0. 

0. 

10 

PU 240 

0. 

0. 

11 

P U 241 

C. 

0. 

18 

SP MAT 

C.10000E-02 

0. 

22 

X F IS2 

0. 14 185E-03 

0. 


MATERIAL 

MATERIAL 

MATERIAL 

MATERIAL 

3 

0 

0 

0 


0.51714E— 31 
0.25857E-01 
0. 

0. 

0 . 

0. 

0 . 

0. 

0 . 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


MATERIAL 

0 


0 , 


FLUXES BY GROUP AND MESH POINT 
FOR THE BASIC CELL 
GROUP 1 

0 • 10043 IE 01 C.9227C1E 00 0.735115E 00 0.528063E 00 0.369309E 00 0.252466E CO 0.161320E 00 0.103651E 00 

0 *64690 IE- 0 1 0 • 31 4333E- 01 0. 152257E-01 0.72561SE-02 0.321260E-D2 0.907240E-03 

GROUP 2 

0.268374E OC 0.2636E3E 00 0.246155E 00 0.202483E 00 0.147527E 00 0.1043C8F 00 0.8014B8E-01 0.830485E-01 

0.966360E-01 0.855960E-01 0.556303E-OI 0.315063E-01 0.154594E-01 0.458762E-02 


****************************************************************************************** 


FDR THE START OF TIME INTERVAL NUMBER 1 OF 2 


FLUX NORMALIZATION FACTOR 0.999738E-09 CONVERTS TO ABSOLUTE FLUX IN NEUT RONS /8 ARN*SEC 
FOR A POWER LEVEL OF 0.649533E 04 WATTS 

CFLL NUMBER 1 OF 2 


****************************************************************************************** 

THE POWER FOR EACH CORE REGION 
0.1 2 50 00 E 01 0.1 25000E 01 

THESE ARE THE RADIAL POWER FACTORS 
0.1 OOOOOE 01 0.1 OOOOOE 01 

THE POWER FACTOR FOR THIS CELL IS 1.0000000 
FLUXES BY GROUP AND MESH POINT 
FOR CELL NUMBER 1 
GROUP l 

0. 1 00405E- 0 8 0.922459E-09 0.734922E-09 0. 527 930E -09 0.369212E-09 0.252400E-09 0.161278E-09 0.103624F-09 

0.6 46311E— 10 0.314251E-10 0. 1 52217E-10 0. 72 5 42 S F- 1 1 0.321176E-11 0.907002E-12 

GROUP 2 

0.26 8304E- 0 9 0.263614E-09 0.246090E-09 0.202430E-09 0.147488E-09 0 . 10428 IE-09 0.801278E-10 0.830267E-10 

0 . 966 10 7E- 1 0 0.8 55 736E- 1 0 0.556157E-10 0. 3 1 4 9 85 E~ 1 0 0.154553E-10 0.458642E-11 


SaF-SH IELD ING FACTORS USED ARE (HIGH TO LOW ENERGY) 


VULCAN 10 * 5 1.000000 1.000000 

VULCAN ID * 19 0.990000 0.800000 

VULCAN ID * 20 C. 990000 0.800000 
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** *♦♦*** ** ** ** ** ** ** HESH INTERVAL NUMBER 


1 ** ** ** ♦* ** *♦ ♦♦ ** *♦ ** 


FISSION PRODUCT CONCEN TRA TIONS AT VARIOUS FRACTIONS OF THE TIME STEP 


FRACTION 

XE135 

SM1A9 

I 135 

PMIA9 


0.01 

0. 702 1A8E- 1 0 

0. 1 2607AE-1 1 

0. 1535A9E-08 

0.36637 6E-09 


0.10 

0* 50558BE- 09 

0.115001E-09 

C. 122953E— 07 

0. 3 5519 8E-C8 


0.50 

0. 1 1 A1A2E- 08 

0.195 72 0E-08 

0.231185E-07 

0* 15535 7E-07 


EQUILIBRIUM 

0 • 123286E- 08 

0.152252E-07 

0. 3039A1E— 07 

0.530516E-07 


TIME AF T ER SHUTDOWN UNTIL MAXIMUM XENON 

CONCENTRATION 

IS REACHED = 10.83 HOURS 

THE XENON- 

135 CONCENTRATION (PEAK) AT THIS TIME IS 0.13A919E-07 


ISOTOPE 

N ( T ) 

N ( T* 1 ) 

N(T+l)/N(T) 

POWER 

POISON 

DESCR. 




FRACTION 

FACTOR 

HYD - l 

0 • 273800F- 01 

0. 273800E-01 

0. lOOOOOE 01 

0. 

0.A9A032E-01 

OXY - 2 

0.I36900E-0I 

0. 136900E-01 

0.100000E 01 

0. 

0 • 55850AE-02 

U— 2 35- 3 

0. 1A1850E- 03 

0. 13871AE-03 

0. 97 78 92 E 00 

0. lOOOOOE 01 

0.A99992E 00 

11- 238- A 

0 . 

0. 

0. 

0. 

0. 

U— 236- 5 

0 . 

0. 5 0990AE— 06 

0. 

0. 

0. 15A853E-0A 

PM1A9- 6 

0. 

0.26521 9E-07 

0. 

0. 

0. 

1-135- 7 

0 . 

0. 3022 3 8E-07 

0. 

0. 

0. 

SM1A9- 8 

0. 

0. 509A7AE-08 

0. 

0. 

0. 198363E-02 

XE 135- 9 

0. 

0. 122601E-0Q 

0. 

0. 

0. 2577666-01 

FP2 -10 

0. 

0. 2621 93E-05 

0. 

0. 

0. 132710E-02 

2ND25-11 

0. 

0. 

0. 

0. 

0. 

BIO 1-12 

0. 100000E-02 

0. 8992 76E-03 

0. 8992 76 E 00 

0. 

0.153938E 02 

BIO 2-13 

0. 100000E-02 

0. 8992 76E-03 

0. 899275 E 00 

0 . 

0. 1 53938E 02 

PU ? 39- 1 A 

0 . 

0 . 

0 . 

0 . 

0 . 

PU2AO-15 

0 . 

0 . 

0 . 

0 . 

0 . 

PU2AI-I6 

0 . 

0 . 

0 . 

0 . 

0 . 

SPMAT-I7 

0. 100000E-02 

0. 877980E-03 

0. 87798 0 E 00 

0 . 

0.18A219E 02 

XF IS2-18 

0.1A1850E-03 

0. 1 3871 AE-03 

0. 977892 E 00 - 

- 0 . 

0.A999S2E 00 


** ** ** ** ** ** ** *+ ** ** HESH INTERVAL NUMBER 2 ** ** ** ** ** ** ** ** ** ** 


FISSION PRODUCT CONCENTRATIONS AT VARIOUS FRACTIONS OF THE TIME STEP 


FRACTION 

XF135 

SM1A9 

I 135 

PM1A9 

0.01 

0. 691652E- 10 

0.1232 37E -1 1 

0. 1 50070E— 08 

0 • 35 8076E-09 

0.10 

C* 502 168E- 09 

0. 1125A8E-C9 

0.120167E-07 

0.3A7150E-08 

0.50 

0. 113A77E-08 

0. 192A20E-C8 

0.27A315E-07 

0.1518376-07 

EQUILIBRIUM 

0 . 1 2257AE- 08 

0*1 51 A5 0£- C7 

0.297055E-07 

0. 5184966-07 


TIME AFT BR SHUTDOWN UNTIL MAXIMUM XENON CONCENTRATION IS REACHED = 10.83 HOURS 


THF XENON- 

135 CONCENTRATION (PEAK) AT 

THIS TIME IS 

0. 

1 31 95AF-07 


ISOTOPE 

DESCR. 

N ( T ) 

N ( T+ 1 ) 

N( T + D/N 

(T) 

POWER 

FRACTION 

POISON 

FACTOR 

HYD - 1 

0 • 273800E- 0 1 

0.273800E-01 

O.IOOOOOE 

01 

0. 

0.A9A573F-O1 

OXY - 2 

0.1369C0E-01 

0. 136900E-01 

0. 1 00000 E 

01 

0. 

0.525265E-0 2 

U— 235- 3 

0. 141850E-03 

0. 138787E-03 

0. 978A09E 

00 

O.IOOOOOE 01 

0.A999S2E 00 

U-238- A 

0. 

0. 

0. 


0. 

0. 

U— 236- 5 

0. 

0. A96078E-06 

0. 


0. 

0. 1 51 52 AG- OA 

PM1A9- 6 

0. 

0.2 5921 OE— 07 

0. 


0. 

0. 

1-135- 7 

0. 

0. 2 953 90E-07 

0. 


0. 

0. 

SM1A9- 8 

0. 

0. 50290 7E-08 

0. 


0 . 

0. 196936E-02 

XE 135- 9 

0. 

0.1218936-08 

0 . 


0. 

0. 2577576-01 

FP 2 -10 

0. 

0. 2 562 60E- 05 

0. 


0. 

0. 130A5 5E- 02 

2N D2 5- 11 

0. 

0. 

0. 


0. 

0. 

BIO 1-12 

0. lOOOOOE- 02 

0. 901 563E-03 

0. 90 1563 E 

00 

0. 

0.15A201F 02 

810 2-13 

0. 100000E-02 

0.901 563F-03 

0.901563E 

00 

0. 

0.15A2016 02 

PU239-IA 

0. 

0. 

0 . 


0. 

0. 

PU2A0-15 

0. 

0. 

0. 


0. 

0. 

PU2A1-16 

0. 

0. 

0. 


0. 

0. 

SPMAT-1? 

0, 10C000E-02 

0.880588E-03 

0. 8 8 05 8 8 E 

00 

0 . 

0*1 8A828E 02 

XF IS2-18 

0. 1A1850E-03 

0. 138787E-03 

0. 978A09E 

00 

-0. 

0.A999S2E 00 



** ** ** ** +* ♦♦ ** ♦* ** ** ME SH INTERVAL NUMBER 


3 ** ** ** ** ** ** ** ** ** ** 


FISSION PRODUCT CONCENTRATIONS AT VARIOUS FRACTIONS OF THE TIME STEP 


FRACTION 

XE 135 

SMI 49 

I 135 

PM149 


0.01 

0.262346E-10 

0.454001E-12 

0.5525826-09 

0.131849E-09 


0.10 

0. 196862E-09 

0. 41 6535E-10 

0.442475E-08 

0. 127826E-08 


0.50 

0.446576E-09 

0. 724495E-09 

0.101 191E-07 

0* 559089E-08 


EQUILIBRIUM 

0.482488E-09 

0. 597372E-08 

0. 1093B IE-07 

0. 19091 9E-07 


TIME AFTER SHUTDOWN UNTIL MAXIMUM XENON 

CONCENTRATION 

IS REACHED * 10.80 HOURS 

THE XENON- 

135 CONCENTRATION (PEAK) AT THIS TIME IS 0.487247E-08 


ISOTOPE 

NIT) 

N(T*1) 

N(T+n/N(T) 

POWER 

POISON 

DFSCR . 




FRACTION 

FACTOR 

HYO - 1 

0.23 643 OE- 01 

0.236430E-01 

0.1000006 01 

0. 

0.128201E 00 

OXY - 2 

0. 2546406-01 

0.254640E-01 

0. 1 OOOOOE 01 

0. 

0.2527976-01 

U-2 35- 3 

0. 

0. 

0. 

0. 

0. 

U-2 38- 4 

0.6771006-02 

0.677014E-02 

0.9998726 00 

0.1181506 00 

0.462777E 00 

U- 236- 5 

0. 

0. 

0. 

0. 

0. 

PM 1 49- 6 

0. 

0. 954454E-C8 

0. 

0. 

0. 

1-135- 7 

0. 

0. 108767E-07 

0. 

0. 

0. 

SMI 49- 8 

0. 

0. 1 923 06E-08 

0. 

0. 

0.2283146-02 

XE 1 35- 9 

0. 

0.4 79799F-09 

0. 

0. 

0. 3076Q5E- 01 

FP2 -10 

0. 

0. 943 6 8 8E— 06 

0. 

0. 

0. 145650E-02 

2ND25-11 

0.491500E-04 

0. 4817026-04 

0. 980065 E 00 

3.859464F 00 

0. 5 19940E 00 

BIO 1-12 

0. 

0. 

0. 

0. 

0. 

BIO 2-13 

0. 

0. 

0. 

0. 

0. 

PU239- 14 

0. 

0. 7353556-06 

0. 

0 • 22383 8€ -01 

0. 172339E- 01 

PU 240- 1 5 

0. 

0. 542 73 9E-08 

0. 

0. 

0. 480696E- 04 

PU241-16 

0. 

0. 4 3 783 3 E- 1 0 

0. 

0.1647996-05 

0.1122806-05 

SPMAT-17 

0. 

0. 

0. 

0. 

0. 

XF IS2- 18 

0 . 

0. 

0. 

-0 . 

0. 


** * + ** ** ** +* ** ** ** «c MESH INTERVAL NUMBER 4 ** ** ** ** ** ** ** ** *♦ ** 


FISSION PRODUCT CONCENTRATIONS AT VARIOUS FRACTIONS OF THE TIME STEP 


FRACTION 
0.0 i 
0 . 10 
0.50 

EQUIL I9RIUM 


XE135 

0. 2273966- 1 0 
0. 18 8422E- 09 
O. 432886E-09 
0. 46804 3E-09 


SMI 49 

0. 365149F-12 
0. 338872E-1 0 
0. bl 59936-09 
0. 583334E-08 


I 135 

0 .^ 38646-09 
0. 3554 1 9E-08 
0.81282 2 E-09 
0. 8 786 03 E —0 8 


PM149 

0. 105908E-09 
0. 102677E-08 
0. 449090 E -08 
0.153 356 E-07 


TIME AFTER SHUTOOWN UNTIL M4 XI MUM XENON CONCENTRATION IS REACHED = 10.71 HOURS 

THE XENON-135 CONCENTRATION (PEAK) AT THIS TIME IS 0.394938E-08 


ISOTOPE 

N ( T ) 

N<T+n 

N ( T+l J / N( T ) 

POWER 

POISON 

oescr. 

HYD - 1 

0 • 23 6430E- 0 1 

0. 2 3 643 OE- 01 

0.1 OOOOOE 

01 

FRACTION 

0 . 

FACTOR 

0.131065E 00 

OXY 

- 2 

0. 2546406-01 

0.254640F-01 

0. 1 OOOOOE 

01 

0 . 

0. 227367E-01 

U-2 35- 3 

0. 

0 . 

0 . 


0. 

0 . 

U-238- 4 

0.677100E-02 

0.6 77032 E— 02 

0. 999900E 

00 

0.105724E 00 

0. 4537S1E 00 

U- 2 36- 5 

0. 

0 . 

0 . 


0 . 

0. 

PM 149- 6 

0. 

0. 76666 8E — 08 

0 . 


0 . 

0 . 

I- 1 35- 7 

0. 

0. 873678F-08 

0. 


0 . 

0 . 

SM 149- 8 

0. 

0. 1 704 60F-08 

0 . 


0. 

0. 208432E— 02 

XE 1 35- 9 

0. 

0. 4654116-09 

0. 


0 . 

0. 3Q7306E— 01 

FP 2 

-10 

0.' 

0. 75821 5E-06 

0 . 


0 . 

0. 120525E-02 

2ND25-11 

0.491500E-04 

0.4835076-04 

0. 9837 38 E 

00 

0.876022E 00 

0.532051E 00 

BIO 1-12 

0 . 

0 . 

0. 


0. 

0. 

810 2-13 

0 . 

0. 

0 . 


0. 

0. 

PU 2 39- 1 4 

0 . 

0. 58523 0E-06 

0 . 


0 • 18253 IE— 0 1 

0.141259E-01 

PU240-15 

0. 

0. 355822E-08 

0 . 


0 . 

0. 3245736-04 

PU241-16 

0 . 

0.2371016-10 

0 . 


0.914434E-06 

0.626222E-06 

SPMAT-17 

0. 

0 . 

0 . 


0 . 

0 . 

XF1S2-1B 

0 . 

0 . 

0 . 


-0. 

0 . 



** ** ** ** ** ** ** ** ** ** MESH INTERVAL NUMBER 


5 ** ** ** ** 4c* ** 4c* ** ** ** 


FISSION PRODUCT CONCENTRATIONS AT VARIOUS 

FRACTIONS 

OF THE TIME STEP 


FRACTION 

X E 1 35 


SM 1 49 

I 135 

PHI 49 


O.Ol 

0. 181 886 E- 

10 

0.2 64636E- 1 2 

0. 32 1 1 26E-09 

0.766226E-10 


0.10 

0. 179272E- 

09 

0. 2491 58E-10 

0.257139E-08 

0. 742847 E-0 9 


0.50 

0.423315F- 

09 

0. 4 8001 6E - 09 

0.5B9060E-08 

0.32490 8E— 08 


EQUILIBRIUM 

0. 45 84 1 3E- 

09 

0. 579243E-C8 

0. 635652E-08 

0. 1 10950E-07 


TIME AFT ER SHUTDOWN UNTIL MAXIMUM XENON 

CONCENTRATION 

IS REACHED = 10.53 HOURS 

THE XENON- 

135 CONCENTRATION (PEAK) AT THIS TIME IS 

0. 2910 75E-08 


ISOTOPE 

NIT) 


N(Tf 1) 

NIT+ll/NIT) 

POW ER 

POISON 

DESCR. 






FR ACT ION 

FACTOR 

HYO - L 

0.236430E- 01 


0.236430E-01 

0. 100000E 

01 

0. 

0.132028E 00 

OXY - 2 

0.25464CE- 01 


0. 2 5464 0E-01 

0. 100000E 

01 

0. 

0.220301E-01 

U-2 35- 3 

0 . 


0. 

0. 


0. 

0 . 

U-238- 4 

0 . 6771 OOE- 02 


0. 67705 1 F-02 

0. 999928 E 

00 

0.102256F 00 

0.4517C3E 00 

U-236- 5 

0 . 


0. 

0. 


0. 

0. 

PM 149- 6 

0. 


0. 554669E-08 

0. 


0. 

0. 

1-135- 7 

0 . 


0. 632 088E-08 

0. 


D. 

0 . 

SM 149— 8 

0. 


0. 140941E-08 

0. 


0. 

0. 17396 1 E- 02 

XE 1 35- 9 

0 . 


0.4 557 85E— 09 

0. 


0. 

0. 30 3706 E- 01 

FP 2 -10 

0 . 


0. 54872 9E-06 

0. 


0. 

0.880470E-03 

2ND25-11 

0. 491500F-04 


0. 485680E-04 

0. 988159E 

03 

3.884419E 00 

0. 537949E 00 

BIO 1-12 

0 . 


0 . 

0 . 


0 . 

0 . 

010 2-13 

0 . 


0 . 

0 . 


0 . 

0 . 

PU239- 14 

0 . 


0. 4240 C5E-06 

0 . 


0 • 133249F-01 

0. 103307E-01 

PU 240- 1 5 

0 . 


0. 1 881 73E-08 

0 . 


0 . 

0.173264E-04 

PU241-16 

0 . 


0. 918375F-U 

0 . 


0 . 356880E-06 

0. 244843E- 06 

SPMAT-17 

0 . 


0 . 

0 . 


0 . 

0 . 

XF IS2-18 

0 . 


0 . 

0 . 

—0 • 

0 . 


** ** ** ** 

** ** ** ** ** ** MESH INTFRVAL NUMBER 

6 ** ** ** ** ** ** ** ** ** ** 

FISSION PRODUCT CONCENTRATIONS AT VARIOUS FRACTIONS 

OF THE TIME STEP 


FRACTION 

XE 135 

SM 1 49 I 135 

PM1 49 


0.01 

0. 13 8902E- 1 0 

0. 1 86272F-12 0.225782E-39 

0 • 5 38729E-1 0 


0.10 

0. 166770E-09 

0. 1 77436F-1 0 0. 180793E-08 

0.5 22292 E-D 9 


0.50 

0.41 1354E-0 9 

0.3 585 88E - 09 0.413462E-08 

0 • 2 28 44 IE-08 


EQU IL I RRIUM 

0. 446543E-09 

0. 5760C8E-08 0.446923E-08 

0. 780084E-C8 


TIME A F T ER SHUTDOWN UNTIL MAXIMUM XENON CONCENTRATION 

IS REACHED = 10.27 HOURS 

THE XENON 

-135 CONCENTRATION (PEAK) AT THIS TIME IS 

0.2 10 290E-08 


ISOTOPE 

NIT) 

NIH-l) NIT+ll/NIT) 

POWER 

POISON 

DESCR. 




FRACT ION 

FACTOR 

HYD - 1 

0. 23 643 OE- 01 

0. 236430E-0I 0.100000E 

01 

0. 

0. 132801 E 00 

OXY - 2 

0.254640E-01 

0. 2 54640E-01 0.100000E 

01 

0. 

0.2145986-01 

U-2 35- 3 

0. 

0. 0. 


0. 

0. 

U-238- 4 

0 . 677 100E-02 

0. 677066E-02 0.999950E 

00 

0.994652F-01 

0.4500C6E 00 

U— 236- 5 

0. 

0. 0. 


0. 

0. 

PM 149- 6 

0. 

0.389985E-08 0. 


0. 

0. 

1-135- 7 

0. 

0. 44441 8E-08 0. 


0. 

0. 

SM 1 49— 8 

0. 

0. 110963E-08 0. 


0. 

0. 137985E— 02 

XE 1 35- 9 

0. 

0. 443 908E-09 0. 


0. 

0. 298086E- 01 

FP 2 - 1 0 

0. 

0.3859046-06 0. 


0. 

0. 623847E-03 

2ND 2 5— 11 

0.491 500E- 04 

0. 48738 7E- 04 0.991632E 

00 

0.891100F 00 

0. 542660E 00 

610 1-12 

0. 

0. 0. 


0. 

0. 

810 2-13 

0. 

0. 0. 


0. 

0. 

PU 2 39- 1 4 

0. 

0. 29843 1 E— 06 0. 


0.943504E-02 

0.7325656-02 

PU 240- 1 5 

0. 

0.937784E-09 0. 


0. 

0. 86995 5E-05 

PU 241- 1 6 

0. 

0.32494 EE-1 1 0. 


0. 127035E— 06 

0. 872817E— 07 

SPMAT-17 

0. 

0. 0. 


0. 

0. 

XFIS2-18 

0. 

0. 0. 


-o. 

0. 



** ** ** ** ** ** ** ** ** *♦ MESH INTERVAL NUMBER 


7 ******** ** ** ** ** ** ** 


FISSION PROOUCT CONCENTRATIONS AT VARIOUS FRACTIONS OF THE TIME STEP 


FRACT ION 

XE135 

SMI 49 

I 135 

PM149 


0.01 

0. 109043E- 1 0 

0.139640E-12 

0. 158999E— 09 

0.4032416-10 


0.10 

0. 15 1570E- 09 

0.13377CE-10 

0. 1353246-08 

0.390937E-09 


0.50 

0.391463E-09 

0.277898E-09 

0. 3D9478E-08 

0.1709896-08 


EQUILIBRIUM 

0.426035E-09 

0. 561 1C3E-08 

0. 3345236-08 

0.583895E-08 


TIME AFTER SHUTDOWN UNTIL 

maximum XENON 

CONCENTRATION 

IS REACH60 * 10 

1.0 2 HOURS 

THE XENON- 

135 CONCENTRATION CPEAKI AT THIS TIME IS 0.161651E-08 


ISOTOPE 

NIT) 

N(T*-1) 

N I Tn»/N(T» 

POWER 

POISON 

DESCR. 
HYD - 1 

0.23 64306- 01 

0. 236430E-01 

0. 100000E 01 

FRACTION 

0. 

FACTOR 

0.1360746 00 

OXY - 2 

0.2546406-01 

0.2546406-01 

0.100000E 01 

0. 

0. 184331E— 01 

U— 235- 3 

0. 

0. 

0. 

0. 

0. 

U— 238- A 

0.6771 00E-02 

0. 677Q75E— 02 

0.999964E 00 

0. 8493656-01 

0.4389796 00 

U— 236- 5 

0. 

0. 

0. 

0. 

0. 

PM 1 49— 6 

0. 

0. 2919056-08 

0. 

0. 

0. 

I- 1 35- 7 

0. 

0. 33264 8E-08 

0. 

0. 

0. 

SM 149— 8 

0. 

0.8 8770 76—09 

0. 

0. 

0. 1140236-02 

X E l 35- 9 

0. 

0 . 42 3 4 4 76 - 09 

0. 

0. 

0.2937086-01 

FP 2 -10 

0. 

0. 2 888926- 06 

0. 

0. 

0. 4823 54E— 03 

2N 025-11 

0.4915006-04 

0. 48 83 72 E -04 

0. 99363 56 00 

0.9079116 00 

0.555429E 00 

BIO 1-12 

0, 

0. 

0. 

0. 

0. 

BIO 2-13 

0. 

0. 

0. 

0. 

0. 

PU 2 39— 1 4 

0. 

0. 2203306-06 

0. 

0 . 715295F— 02 

0. 5586566-02 

PU 2 40- 15 

0. 

0.532431E-09 

0. 

0. 

0. 510182E-05 

PU 241- 1 6 

0. 

0. L42090E-11 

0. 

0. 570408E— 07 

0. 394222 E— 07 

SPMAT-17 

0. 

0. 

0. 

0. 

0. 

XFIS2-18 

0. 

0. 

0. 

— 0 • 

0. 


** ** ** ** 

** ** ** ** ** ** 

MESH INTERVAL 

NUMBER 

3 ** ** ** ** ** ** 

** ** ** ** 

FISSION PRODUCT CONCENTRATIONS AT VARIOUS 

FRACTIONS 

OF THE TIME STEP 


FRACT ION 

XE 135 

SMI 49 

I 135 

PM149 


0.0 1 

0. 106788E- 10 

0.1 37407E-12 

0. 16S463E-39 

0. 397189E— 10 


0.10 

0. 14 5596E- 0 9 

0. 131649E-L0 

0. 133293E-08 

0.38507QE-09 


0.50 

0, 373428E- 09 

0.272579E-09 

0- 3048336-08 

0. 168423E-C8 


EQU ILIBR IUM 

0 . 40 1 2506- 09 

0. 533385E-08 

0.3295036-08 

0. 575133E-08 


time AFTB* SHUTDOWN UNTIL 

. MAXI MUM XFNON 

CONCENTRAT ION 

IS REACHED = 10 

i.05 HOURS 

THE XENON 

-135 CONCENTRATION (PEAK) AT THIS TIME IS 

0. 158600E-08 


ISOTOPE 

NIT) 

NIT* 1) 

N(T*1)/N(T) 

POWER 

POISON 

DESCR. 





FRACTION 

FACTOR 

HYD - 1 

0.23 64 30E- 01 

0. 2364306-01 

0. 1000006 

01 

0. 

0 . 14 2640 E 00 

OXY - 2 

0.25464CE-01 

0.2546406-01 

o. loooooe 

01 

0. 

0. 1216606-01 

U— 235- 3 

0. 

0. 

0. 


0. 

0. 

U-238- 4 

0 . 6771 00E-02 

0. 677077E—02 

0. 99 9966 6 

00 

0 • 554 120E—01 

0.415632E 00 

U- 2 36- 5 

0. 

0. 

0. 


0. 

0. 

PM 149- 6 

0. 

0. 287524F-0B 

0. 


0. 

0. 

1-135- 7 

0. 

0. 32 76566-08 

0. 


0. 

0. 

SMI 49- 8 

0 . 

0. 86731 86-09 

0 . 


0. 

0. 118577E-02 

XE 1 35- 9 

0 . 

0. 4037936-09 

0. 


0 . 

0.2981C8E-01 

FP 2 -10 

0 . 

0.2845526-06 

0 . 


0 . 

0. 505739E- 03 

2ND25-11 

0.491500E-04 

0. 48832 76-04 

0.993544E 

00 

0.937421E- 00 

0.5787006 00 

810 1-12 

0. 

0. 

0 . 


0 . 

0 . 

BIO 2-13 

0 . 

0 . 

0 . 


0 . 

0 . 

PU 2 39- 1 4 

0. 

0. 2 098406-06 

0 . 


0.7167346-02 

0.566315E-02 

PU 2 40- 1 5 

0 . 

0.5253766-09 

0 . 


0. 

0. 5358356—05 

PU241-16 

0. 

0. 1452406-11 

0. 


0. 6134296-07 

0. 42S9C5E— 07 

SPMAT-17 

0. 

0 . 

0 . 


0 . 

0 . 

XFIS2-18 

0. 

0 . 

0 . 


- 0 . 

0 . 


25 


** ** ** ******** ** ** ** MESH INTERVALS NUMBER 


9 THRU 1+ ** ** ** ** ** ** ** ** ** *♦ 


NO FISSIONABLE ISOTOPES OR CONSEQUENT FISSION PRODUCTS OCCUR IN MESH INTERVALS 9 THRU 


ISOTOPE 

N ( T ) 

N( T+ 1 ) 

DESCR. 



HYD - l 

0.517140E-01 

0. 5I7140E-01 

OXY - 2 

0*2 585 70E- 01 

0* 2 58570E-01 

U-2 35- 3 

0. 

0. 

U-238- 4 

0. 

0. 

U— 236- 5 

0. 

0. 

PM149- 6 

0. 

0. 

1-135- 7 

0. 

0. 

SM 1 49— 8 

0. 

0. 

XE 1 35- 9 

0. 

0. 

FP2 -10 

0. 

0. 

2ND25-11 

0. 

0. 

BIO 1-12 

0. 

0. 

BIO 2-13 

0. 

0. 

PU239-14 

0. 

0. 

PU 240- 1 5 

0. 

0. 

PU241-16 

0. 

0. 

SPMAT-17 

0. 

0. 

XFIS2-18 

0. 

0. 


N(T+n /NCTI POWER 

FRACTION 

0.100000F 01 0, 0. 

0. 100000 F 01 0. 0. 

0. 0. 0. 

0. 0. 0. 

0* 0 » 0. 

0 # 0 * 0 * 

0. 0. 0. 

0. 0 * 0* 

0. 0. 0. 

0 . 0 . o. 

0. 0, 0. 

0. 0, 0. 

0. 0, 0. 

0. 0. 0. 

0* 0 * 0. 

0. 0. 0. 

0. 0, 0, 

0. 0, 0* 


ALL MESH INTERVALS HAVE BEEN TRAVERSED FOR TIME INTERVAL 1 
THE FLUX WILL NOW BE RENORMALIZED FOR THE NEXT TIME INTERVAL 


POISON 

FACTOR 


*************** ****************************************** ********************************* 


FOR THE START OF TIME INTERVAL NUMBER 2 OF 2 


FLUX NORMALIZATION FACTOR 0.101036E 01 CONVERTS TO ABSOLUTE FLUX IN NEUTRONS/BARN*SEC 
FOR A POWER LEVEL OF 0.649533E 34 WATTS 

CELL NUMBER 1 OF 2 


*************** *************************************************************************** 


NEW ATOM DENSITIES BY MESH INTERVAL FOR EACH ISOTOPE AT THE END OF TIME INTERVAL 1 
******************************************************************************** 


LAST ENTRY FOR EACH ISOTOPE IS ONE BARN-TH OF ITS TOTAL NUMBER OF ATOMS IN THE CONFIGURATION 


******************************************************************************** 

VULCAN ID = 23 

0.273800E-01 0.273800E-01 0.236430E-01 Q.236430E-01 0.236430E-01 0.236430E-C1 

0 . 236430E- 01 

0 • 236 430E-01 

0.5 17 140E-0 1 

0.517 140E- 01 

0.51 71 40E- 01 

0.517 140 E—01 

0.51 71406-3 l 

0.517 140E-01 

0.248069E 01 


VULCAN ID = 
0.1 36900 E- 0 1 

24 

0 • I 36900E- 01 

0.254640E-01 

0.254640E-01 

0.254640E-01 

0.2546406-01 

0. 25 46 4 OE “01 

0.254640E-01 

0.258570E-01 

0.258570E-01 

0.258570E-01 

0.2585 73 E-01 

0.258570E-01 

0.258 5 70E-01 

0.165617E 01 


VULCAN ID = 
0*1 387 14E- 03 

5 

0 • l 38 787E- 03 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 105737E- 0 2 


VULCAN ID = 

0. 

7 

0. 

0. 67701 4E-02 

0.6 77 032E— 02 

0.6770516-02 

0.677066E-P2 

0.6770 75E- 02 

0. 677077E— 02 

0. 

0. 

0. 

0. 

0. 

0. 

0* 206365E 00 


VULCAN ID = 
0.5C9904E-06 

6 

0 .4 c 6 C 78E“ 06 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0.381523F-05 


VULCAN ID = 
0.265219E-07 

17 

C.259210E-C7 

0.9 54454E- 08 

0.76S668E-08 

0.5 54669E-0 8 

0.389985E-08 

0.291905E-08 

0.287 524E-08 

0. 

0. 

0. 

0. 

0. 

0. 

0. 370043E- 06 


VULCAN ID = 
0.3D2238E-07 

16 

0.295390E-07 

0.1087676-07 

0. 873 678E—08 

0.632088E-08 

0 .4444I8E— 08 

0.332648E-08 

0.327656E-08 

0. 

0. 

0. 

0. 

0. 

0. 

0.421693E-06 


VULCAN ID = 
0.509474E-08 

15 

0.502907E-08 
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0.1 70460E— 08 

0.1 40941 E-08 

0.1 10963E— 0 8 

0. 8 87 707E-09 

0.86731 8E— 09 

0- 

0. 

0. 

0. 

0. 

0. 

0.79811 3E-07 
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VULCAN 10 * 

14 



0.122601E-C8 

0. 12189 3E-08 

0. 479799E— 09 

0.465 41 1 E—09 

0. 

0. 

0. 

0. 

VULCAN ID = 

13 



0 .26219 3E- 0 5 

0.256260E-05 

0. 943688E — 06 

0.75821 5 E— 06 

0. 

0. 

0. 

0. 

VULCAN 10 = 

21 



0. 

0. 

0. 461702E-04 

0.483507E-04 

0. 

0. 

0. 

0. 
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19 
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0. 

0. 

0. 

0. 

0. 

0. 
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0. 

0. 
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0. 

0. 54273 9E-08 

0. 3 5582 2 E— 08 

0. 

0. 

0. 
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VULCAN ID = 
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0. 

0. 
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0.237101 E— 1 0 

0. 
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0. 

0. 
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18 



0.87 7980E- 0 3 

0 . 8 80588E- 03 

0. 

0. 

0. 

0. 

0. 

0. 

VULCAN ID = 
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0. 138714E-03 

0.138 787E-03 

0. 

0. 

0. 

0. 

0. 

0. 


0. 455785E— 09 

0 • 4 43 90 8 F— 09 

0.423447E-09 
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0. 

0. 

0.2 30056E- 07 
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0 . 385904E— 06 

0.288892E-06 

0.2845526-06 

0. 

0. 

0.36591 8E-04 


0.485680E-04 

0 • 487387E-04 

0.488372E-04 

0. 48 832 7t— 04 

0. 

0. 

* 0. 1 48020E- 02 


0. 

0. 

0. 

0. 

0. 

0. 

Q.68641QE-Q2 


0. 

0. 

0 . 

0. 

0. 

0. 

0.686410E-02 


0. 4 2 40 0 5 E— 06 

0 • 298431 E— 06 

0. 220330E- 06 

0.2 09840E— 06 

0. 

0. 

0.1 30 59 6E -04 


0. 180173E-O8 

0.937784E— 09 

0. 532431E- 09 

0.5253 76E-09 

0. 

0. 

0.6 0805 9F- 07 


0.918375E-U 

0.324 94 8F- 1 1 

0 . 142090E- l 1 

0.L4524OF-11 

0. 

0. 

0. 44348 7E-09 


0. 

0 . 

0. 

0. 

0. 

0. 

0.6 7034 6E -02 


0. 

0. 

0. 

0. 

0. 

0. 

0. 10573 7E-02 



CONVERSION RATIO BY MESH INTERVAL AT START OF TIME INTERVAL 1 


- 0 . 

0. 


- 0 . 

0. 


0 ■ 75931 2 E 00 0.739331F 00 0.733734F 00 0.729275E 00 0.707063F 00 0.664022E 00 

0 . 0 . 0 . 0 . 


CONVERSION RATIO BY MESH INTERVAL AT END OF TIME INTERVAL 1 


-0. 

0 . 


- 0 . 

0. 


0. 749893E 00 0.732101E 00 0.728515E 00 0.725613E 00 0.704489E 00 0.661847E 00 

0 . 0 . 0 . 0 . 


BY MESH INTERVAL THE POWER IS 

• 4 627 17E 03 0.452355E 03 O.168015E 03 0.135525E 03 0.979981E02 0.688741E 02 0.515376E 02 0.507576E 02 

. 0. 0. 0. 0. 0. 0.6495336 04 

BY MESH INTERVAL THE POWER (NORMALIZED TO VOLUME-AVERAGED POWERI IS 

0.271419E OL 0.26534 OE 01 0.990226E 00 0.794954E 00 0.574833E 00 0.403998E 00 0.302307E 00 0.297731E 00 

0. 0. 0. 0. 0. 0. 0.1704816 03 

VOLUME AVERAGFD EQUILI3RIUM CONCENTRATIONS OF MAJOR FISSION PRODUCTS 
IN NUCL E I/BARN*CM IF THE POWER LEVEL OF .6495336 04 WATTS IS MAINTAINED 


XENON- 135 
SAMARIUM-149 
100 IN E— 135 
PROM FT HIUM— 149 


0.607281 E— 09 
0. 763 967E-08 
0. 111305 E-07 
0.194277E-07 
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FISSIONABLE 

ISOTOPE 

U-235 

U-238 

U-236 

PU2 39 

PU240 

PU 241 


INITIAL NASS IN 
KILOGRAMS 

0 ■ 142829E- 02 

0.815713E-01 

0 . 

0 . 

0 . 

0. 


FINAL MASS IN 
KI LOG RAMS 

0.1402 95E-02 

0. 81 5653E-01 

0.149527E-05 

0. 5 1 83 5 IE— 05 

0.2742436-07 

0. 177500E-09 


FRACTIONAL 

POWER 

0.944316E 00 
0.483206E-01 
0 . 

0.7 36294 E-02 
0 . 

0.370968 E-06 


THE FOLLOWING EDIT IS AVERAGED OVER ALL MESH INTERVALS IN WHICH EACH ISOTOPE OCCURS 


ISOTOPE 

NIT) 

NIT+l) 

NlT+ll/NlT) 

POISON 

DESCP . 





FACTOR 

HYD 

-23 

0.3642 72E- 0 L 

0. 3642 72E-01 

0* 100000 E 

01 

0 • 563034E-01 

OXY 

-24 

0 • 243197E- 0 1 

0.24319 7E— 01 

0. 100000 E 

01 

0 . 837446E— 02 

U-235- 5 

0.14 l B50F- 03 

0. 138763E-03 

0.978237E 

00 

0 • 196461 E 00 

U-238- 7 

0.6771 OOE- 02 

0.677051E-02 

0, 99992 7 E 

00 

0 • 126868E 00 

U-236- 6 

0. 

0. 500686E-06 

0. 


0 . 5998 1 1F-05 

PM 149-17 

0. 

0.971242E-08 

0. 


0. 

1-135-16 

0. 

0.U0681E-07 

0. 


0. 
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0. 

0.2094 79E-08 

0. 


0 • 129656E-0 2 

XE 1 35- 14 

0. 

0.603 82 0E-09 

0. 


0 • 186447E-0 1 

FP2 

-13 

0. 

0. 96041 3E-06 

0. 


0 • 807945F-03 

2ND2 5-21 

0 . 491 50QE- 04 

0.485632E-04 

0.988060E 

00 

0 . 149980E 00 

BIO 1-19 

0. 1000C0E-02 

0. 900801F-03 

0. 9008 01 E 

00 

0 .605550 E 01 

BIO 2-20 
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0. 900801E-03 

0. 9008 01 E 

00 

0 . 605550E 01 

PU 239- 9 

0. 

0.428466E-06 

0. 


0.343608E-02 

PU 2 40- 1 0 

0. 

0 . 22 5741 E- 08 

0. 


0 • 764434F— 0 5 

PU 24 1-1 1 

0. 

0. 1 45 501 E-l 0 

0. 


0. 153195F-06 

SPMAT-18 

0 . 1 000 OOE- 02 

0. 87971 9E-03 

0.879719E 

00 

0.725433E 01 

XF IS2-22 

0.141 850E- 03 

0. 138763E-03 

0.978237E 

00 

0 . 196461 F 00 



INI TIAL 

FINAL 

TOTAL FUEL INVENTORY 

0. 829996E-01 

0 • 8 29789 E-01 

IN K IL 3 GRAMS 



TOTAL FUEL INVENTORY 

0. 182 981 E 00 

0.182935E 00 

IN POUNDS 



CONVERSION RATIO 

0. 202256E 00 

0.205006E 00 

TOTAL POWER IN WATTS 

0.649533E 04 

0 • 642 872 E 04 


MEGAWATT-DAYS/METR 1C TON 


0.654201E 06 
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VULCAN DESCRIPTION 


VULCAN is a one -dimensional point (mesh interval) depletion or burnup code 
primarily designed to use the GAM-GATHER-TDSN calculational system and its associ- 
ated input and output. The input for VULCAN requires geometry and material specifica- 
tions, macroscopic (atom density in atoms/ (b) (cm)) specifications, microscopic cross 
sections in barns, and the space -energy flux, which is normalized, within VULCAN, to 
units of neutrons per barn per second. Variable input format options are available for 
the cross sections and the flux making VULCAN readily adaptable to any space -energy 
calculation sequence. Thus any reactor configuration, for which geometry and materials 
can be specified, microscopic cross sections obtained, and fluxes provided, may be 
treated. Input shortcuts are provided if several configurations, differing only in isotopic 
densities and flux shape, are considered. In particular, power factors may be used in 
order to adjust the flux levels. If, in a transverse dimension, a similar material con- 
figuration exists, then several transverse dimension regions may be treated sequentially 
within one VULCAN problem. Thus some measure of two-dimensional depletion is 
available. The program can treat up to 40 nuclides and 20 energy groups and can ac- 
commodate 190 mesh intervals. The mesh intervals may be grouped into 190 zones 
comprised of 190 distinct materials. These restrictions result because VULCAN is 
limited to storage in the fast memory of an IBM 7094 computer. However, minor re- 
programming can greatly increase the range of problems that may be treated. The 
restrictions on the number of mesh intervals, zones, and materials NRM and energy 
groups NG may be varied subject to the following condition: 

52 * NRM + NRM * NG + 112 * NG ^ 16 040 

The COMMON blocks /SET 1/, /SET 2/, and /SET 3/ must be suitably altered in each 
subroutine and the source decks recompiled. The storage limitation (16 040) is dependent 
on the storage allocation for internal systems routines. 


Description by Subroutine 

Main program: VULCAN. - This is the main program and it performs virtually all 
input-output operations. It controls the flow of the program and performs special edits 
of power and conversion ratios in the form of plots at the end of each time interval T^. 
The conversion ratio CR m in any mesh interval m is the ratio of the creation rate 
C m to destruction rate D m of fissionable material in that mesh interval. The equation 
for CR m is 
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The total conversion ratio CR over all mesh intervals NRM is 

NRM 

£ C m • (Volume) m 

CR = JBsl 

NRM 

2 D m • (Volume) m 
m=l 

New geometry and material specifications are punched by this routine (discussed in the 
section Subroutine OUTPUT). 

Sub routine TABLE. - This subroutine contains tabular information such as, fission 
product yields, decay constants, atomic weights, thermal v values (thermal v = 
average number of neutrons per thermal fission), and the VULCAN internal identification 
assignments. It contains a constant, 3. 1x10^ fissions per thermal watt-second, to be 
used for absolute flux normalization. The tabular information is taken, for the most 
part from references 10 to 12. Fission product yields of Xe^*^, 1*^, and 

Pm^^ from fast and thermal fission of U^^, U^^, Pu^^, and Pu^^ are included to 
the extent they are available. These tabulated yields will be used unless different values 
are specified in the input. The FORTRAN deck of this subroutine is short; thus any of 
the tabulated values may be changed and the deck recompiled with little effort. 

Subr outine ABSPHI . -Within ABSPHI, prior to the solution of any depletion equations, 
the arbitrary level input flux is normalized to a specified power. At the end of the de- 
pletion time interval the flux is renormalized to the same power before calculating the 
final conversion ratios and before traversing another depletion time interval. 

Subroutine XENON. - The equilibrium concentrations of the four specific fission 
products (XE*^, Sm^^, 1*^, and Pm^^) are calculated here as well as the time after 
shutdown (at the end of the current time interval Tj) of peak Xe 1 0 concentration and its 
value at that time. 

Subroutine SSF . - If isotopic self -shielding factors by energy group are desired, they 
are read in by this subroutine. The absorption and fission cross sections are multiplied 
by their respective self -shielding factors. 
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Subroutines BU RNUP and FISEQ. - These two subroutines perform the actual deple- 
tion calculation. For each isotope in each mesh interval, the new atom density, frac- 
tional change in atom density, fractional power supplied, and a poison factor are calcu- 
lated. The poison factor for an isotope is the ratio of its macroscopic absorption cross 
section to the total macroscopic absorption cross section of the fuel. For each mesh 
interval, the power is calculated; after normalization to the volume-averaged power, a 
power distribution plot is produced in the main program VULCAN. A fuel inventory is 
carried out, volume -averaged atom densities obtained, and the conversion ratio deter- 
mined before and after the depletion time step. 

The subroutine allows the concentrations of the four specific fission products, 
within any mesh interval, to be determined at an arbitrary number of fractional times 
within any time interval Tj. For these specific fission products, equilibrium concentra- 
tions are determined for each mesh interval and volume averaged over all mesh intervals. 
For nondepletable regions, the atom densities for the first mesh interval within such a 
region are carried through as representative of the full region. 

Detailed depletion information may be obtained for each mesh interval separately in 
addition to the information averaged over all mesh intervals. Upon completing the de- 
pletion calculation for all mesh intervals, the flux is renormalized (using new atom 
densities) to the specified power and a new time step T^ within VULCAN is initiated. 

Subroutine OUTPU T. - In a point (mesh interval) depletion program each point 
within a zone (constant material) will usually be depleted individually. Consequently, each 
point becomes a distinct material and zone in the next spatial calculation. It is desirable, 
then, to have some means of minimizing the number of materials or zones. Once all the 
desired time steps T^ are completed within VULCAN, the macroscopic absorption 
cross sections are subjected to two comparative tests. One test determines whether the 
macroscopic absorption cross section (macs) of each isotope in a mesh interval is, 
within a predetermined limit, the same as the corresponding isotopic macs in a preceding 
mesh interval. The other test determines whether the macs of each isotope in a mesh 
interval constitutes a certain fraction of the total fuel macs in that mesh interval. Each 
test is performed for each energy group. The mathematical representation of these 
tests is 

Test 1: 
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Test 2: 



NF 

< 10 €^ 

i=l 



( 18 ) 


where 

macroscopic absorption cross section for isotope k, energy group j, and mesh 
•* m interval m 

e predetermined limit (see DIFFER in the section INPUT INSTRUCTIONS) 
n designates which preceding mesh interval is being compared with mesh interval m 
i designates fissionable isotope 

If all energy groups for all isotopes satisfy both tests for two or more mesh intervals , 
the atom densities for these intervals are volume averaged and thus constitute a single 
new material and zone. After performing these two tests through all mesh intervals, new 
atom densities for each distinct material are punched. These atom densities are used to 
obtain new material macroscopic cross sections which will then be used in the spatial cal- 
culation (see the section Depletion cycle). 

This subroutine, on option, will print out macroscopic absorption cross sections by 
group, mesh interval, and isotope. Upon return to the main program the geometry and 
material specifications for input to TDSN are punched corresponding to the new region and 
material structure. The corresponding output for nondepletable regions maintains the 
atom densities and mesh structure unaltered and alters, not the material, but its external 
identification number. 


Nuclide List 

VID numbers 1 to 17 may be used only for the indicated isotope. If any of these 
isotopes is not to be considered, its VID number may not be used by another isotope. 
Nuclides other than the 17 specifically mentioned are assigned VID numbers in unbroken 
numeric order, beginning with VID = 18. Table II lists the VID number assignments 
which VULCAN uses. 
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TABLE H. - NUCLIDE LIST AND VID ASSIGNMENT 


Nuclide 

VID number 

Th 232 

1 

Pa 233 

2 

U 233 

3 

U 234 

4 

n 235 

5 

^36 

6 

^38 

7 

Np 239 

8 

Pu 239 

9 

Pu 240 

10 

Pu 241 

11 

Pu 242 

12 

Fission product aggregate 

13 

Xe 135 

14 

Sm 149 

15 

jl35 

16 

Pm 149 

17 

Special Materials including additional 

a 17 + 1 - 17 + J 

fission product aggregates 


Burnable Poisons 

a 17 + J+ 1 — 17+J + K 

Repeated Fissionable Nuclides 

a 17 + J+ K+1 — 17+J+K+L 

Nondepletable Nuclides 

a 17 + J+ K+ L+1 — 17 + J + K+ L+ M 


a For a configuration containing J Special Materials, K Burnable Poisons, L 
Repeated Fissionable Nuclides, and M Nondepletable Nuclides where 
J + K+L+M< (40- 17). 
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INPUT INSTRUCTIONS 


This section contains the input instructions and explanations of the input parameters. 
The symbol * after a card number means to use as much of the card or as many cards 
as necessary. 


Card Format 

Variable 

Description 

1* 11, IX, 14A5 

TITLE 

Title, label, or problem description containing 
any alphameric information. The number 1 in 
card column 1 will signify the last title card. 

2 7110 

NG 

Number of energy groups ^ 20 


NFG 

Number of fast energy groups s 20 


NISOT 

Number of isotopes whose depletion and/or buildup 
is to be considered plus those classified as non- 
de pie table ^ 40. 


NMAT 

Number of separate materials < 190 (A material is 
any distinct macroscopic cross section set. ) 


NTINC 

Number of equal time intervals (DELHR) for which 
flux normalization and subsequent depletion will be 
performed. 


NFRACT 

Number of fractional steps to be taken within a 
major time interval (DELHR) for intermediate 
poison calculations s 10. NFRACT = 0 if 
MOREOUT = 0 or 1. 


KCELL 

= 0 or 1 Normal single depletion 


2 < KCELL s 20 A series of KCELL successive 
problems with distinct atom densities and fluxes 
will be depleted. 

>21 A series of KCELL - 20 successive problems 
differing only in flux level will be depleted 
(treated by power factors to be supplied). 

3 7110 NZONR Number of zones < 190 - within each zone the 

material is constant. 

NRM Number of mesh intervals ^ 190 
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Card Format Variable Description 

MOROUT = 0 Only averaged information (always obtained) at 

the end of depletion sweep will be printed. 

= 1 Macroscopic absorption cross sections (with self- 
shielding factors included) will be printed. 

= 2 Individual mesh interval information will be 
printed. 

= 3 Both macroscopic absorption cross sections (with 
self -shielding factors included) and individual mesh 
interval information will be printed. 

NUF = 0 Both p<jf and v will be read in {v is average 

number of neutrons per fission - by group and by 
isotope). 

= 1 Only erf will be read in. 

4 7110 NSPMAT Number of special materials (see eq. (14)). If more 

than one fission product aggregate is to be 
considered, the number must be included in 
NSPMAT. 

NBPOI Number of burnable poisons considered 

NFIS Number of fissionable isotopes for which fission 

product yields are to be used (indicated by Y in 
the output) ^ 10 (NFIS does not include any 
Repeated Fissionable Nuclides. ) 

NYDNUC Number of fission product nuclides for which yields 

are to be read in ^ 20 (Total number whose yields 
are to be considered is equal to NYDNUC plus the 
number whose yields are to be taken from 
subroutine TABLE.) 

NONDPL Number of nondepletable isotopes (carried through 

to maintain continuity of the calculational 
sequence) 

NFIRPT Number of times that one or more of the twelve 

fissionable isotopes is repeated. For example, if 

235 

three distinct cross section sets for U and two 
distinct sets for U are to be used in a single 
problem, then NFIRPT = (3 - 1) + (2 - 1) = 3. 
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Card 


Format 


Variable 


5 7110 


6 7E10. 5 


KGEO 


KFLUX 


KSSF 


KHAIN 


KFAST 


KFORM 


POWER 

DISTI 


DIST2 


Description 

Geometry indicator (one -dimensional): 

= 1 Slab 
= 2 Cylinder 
= 3 Sphere 

Format for flux input: 

= 1 Binary -TDSN output 

= 2 Decimal -NRM values for group 1, NRM values 
for group 2, and so forth (NG sets of NRM values) 

Self -shielding factors: 

= 0 If none 

= 1 If factors are to be used 

Fission isotope chain applicable: 

= 0 Both chains 
= 1 Th 232 - U 236 chain 
= 2 U 238 - Pu 242 chain 

= 0 For thermal reactor 
= 1 For fast reactor 
(This option chooses fission product 
yields for two general reactor classes.) 

= 0 KFLUX option is used. 

= 1 Permits optional FORMAT statement to be read 
in for cross section input 

= 2 Permits optional FORMAT statement to be read 
in for flux input 

= 3 Permits optional FORMAT statement to be read 
in for both cross section input and flux input 

Total power in watts (thermal power of 
1 W = 3. 1x10^ fissions/sec) 

First transverse dimension equals height for cy- 
linder or slab in centimeters (used to calculate 
power density) 

Second transverse dimension equals thickness for 
slab in centimeters (used to calculate power 
density) 
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Card Format Variable Description 

DELHR Time increment in hours Tj (Total depletion time 

is NTINC *DE LHR . ) 

XREMP Xenon removal rate £ (other than by decay and 

depletion, e.g. , diffusion). 

REMIOD Iodine removal rate 

DIFFER Averaging criterion for output materials. 

If DIFFER = 0, DIFFER will be set equal to 0. 001 
(0. 1 percent). See e in inequalities (17) and (18). 

Card 7 is included only if NFRACT ^ 0. 

7* 7F10. 8 FRACT(I) Fractional values of DELHR to be used - no restric- 

tion on values except that NFRACT (^ 10) values 
are necessary. 

Isotope Identification 

One card 8 must be included for each of NISOT isotopes. Any isotope may have any 
sequential ID number (I), but cards 8 must be in numerical order by (I). 

8* 415, A5 I Sequential ID number 

VID(I) VID number (see the section Nuclide List) 

NRPT(I) = 0 If not a repeated fissionable nuclide 

= 10 + VTD If a repeated fissionable nuclide 
235 

If U were used twice in addition to the initial 
occurrence, each repeated time NRPT would be 
10 + 5 = 15. Nuclides for which NRPT ^ 0 must be 
in ascending NRPT number. 

KFIS(I) = 0 If no fission product yield is to be used or if 

NRPT (I) ^ 0 (Note that, upon completion of the 
input, the KFIS(I) array will contain NFIS nonzero 
values. ) 

= 1 If a fission product yield is to be used 

DESCR Arbitrary alphameric description of VTD(I) 
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Card 


Format Variable 


Description 


Geometry and Material Specifications (NZONR Cards 9 Required) 

NZONR cards 9 are required. The first card is for zone K = 1; the second card is 
for zone K = 2, and so forth; the last card is for zone K = NZONR. 

9* 15, F10. 6, 15 NMIR(K) Number of mesh intervals in zone K 

RAR(K) External dimension of zone K (cm) 

MIDRZ(K) Material number corresponding to zone K 

If NYDNUC = 0, skip card 10. If NYDNUC / 0, NYDNUC cards 10 are required. 

10* 215, 6F10. 8 KRON = 0 If this isotope is not produced by decay of 

preceding isotope in NYDNUC series (always 0 for 
first card in series) 

= 1 If isotope is capture product of preceding nuclide 

IDFP VID number of the particular isotope (e. g. , if only 

Xe 135 yields are to be read in IDFP = 14) 

YELD(J) Fractional fission yield of isotope with VID = IDFP. 

YELD(l) from first nonzero KFIS isotope, YELD(2) 
from the second nonzero KFIS isotope, and so 
forth - supplied in the same sequence as (I) on 
cards 8*. If NFIS ^ 7, continue the yields on the 
next card (format - 7F10.8). 

If NUF = 1, skip cards 11. If NUF - 0, NFIS cards 11 are required. 

For each fissionable isotope considered and if NUF = 0, the values of v are read in 
for all fast groups. These groups are ordered from high to low. Note that cr^ rather 
than Wj is used in VULCAN. The card sequence is ordered according to VID number. 

11* 7F10. 8 XNU(IG) NFG values of u in neutrons per fission. Do not 

include for the NFIRPT isotopes. 

If KFORM = 0 or 2, skip card 12. 

12 12A6 FMT(l) Optional FORMAT statement for cross section input. 

For the specific format of cards 14, card 12 
would read (36X, 2F12. 8). Note omission of word 
FORMAT but the inclusion of parentheses. 
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Card 


Format 


Variable 


Description 


Microscopic Cross Sections 

Include NISOT sets of cards 13 and 14 in the same sequential order as (I) on cards 8. 

13 12A6 DUMMY Not retained in VULCAN (may be title card for cross 

sections) 

If KFORM = 1 or 3, read in absorption and fission cross sections for NG groups in 

specified format (card 12) on cards 14. For each card 13, there must follow NG cards 14. 

14* 36X, 2F12. 8 ABSIG(IG) Absorption cross section for group IG in barns 

(or format 
indicated on 
card 12) 

FUNSIG(IG) ctj or vo^ cross section in barns (depending on 
availability) with appropriate values of XNU(IG) 
on cards 11 and NUF on card 3. 

Material Specifications (Atom Densities) 

There must be NMAT sets of cards 15 and 16. The NMAT sets must be in numerical 

order by material number. 

15 7110 NCON Number of constituent nuclides in this material =£ 40 

(Only those nuclides present at the start of the 
depletion interval need be included. ) 

IDCON(l) Sequential ID numbers in ascending order (as 

established by (I) on cards 8) - NCON values. If 
NCON > 6, continue with 7110 format on as many 
cards as are needed. 

16* 7F10.8 DENS(I) Corresponding atom densities for this material 

(NCON values) in atoms per barn-centimeter 

If KFORM = 0 or 1, skip card 17. 

17 12A6 FMT(I) Optional FORMAT statement for flux input 
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Card 


Format 


Variable 


Description 


Fluxes 


The fluxes are supplied on cards 18. 

18* 7E10. 5 FLUX(I) Average fluxes by mesh interval for each group; 

NG sets of NHM values (normalized to units of 
neutrons/(b)(sec) within VULCAN) 


Power Factors 

If KCELL ^ 20, skip cards 19 and 20. If KCELL > 21, the power factors are supplied 
on cards 19 and 20. 

19 110 NKOREG Number of successive problems to be depleted using 

the same fluxes but adjusted by power factors 
< 19 

20* 7E10.5 CORPOW(I) Actual, fractional, or relative power for successive 

problems; NKOREG values; the normalization is 
to the input, that is, 

POWER _ * cORPOW(I) = actual power in pro- 
CORPOW(l) 

blem (I) 

Self-Shielding Factors 

If KSSF > 0, read in the self-shielding factors on cards 21, 22, and 23. 


21 

110 

NSSF 

Total number of nuclides for which self-shielding 
factors are to be considered ^20 

22* 

7110 

JSMAT(I) 

VID numbers for self -shielded nuclides - NSSF 
values 

23* 

7F10.6 

SSF(I, J) 

Self -shielding factors for NSSF isotopes (I) and NG 
groups (J); NSSF sets of NG values in the order 
established by JSMAT(l) (card 22) 


If KCELL > 21, read in KCELL - 21 sets of atom densities and self -shielding factors. 
Each set of atom densities is composed of the NMAT pairs of cards 15 and 16 and must 
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be similar to the first set to the extent that only density magnitude changes (i. e. , no 
changes in material constituents) are allowed. Different self-shielding factors may be 
used for succeeding sets. If 2 < KCELL < 20, read in KCELL - 1 additional sets of 
similar atom densities, fluxes, and self -shielding factors (cards 15, 16, 17 (if used 
before), 18, 21, 22, and 23). 


LISTING OF PROGRAM DECKS 

The FORTRAN listing of the VULCAN program which follows contains several par- 
ticular features of the Lewis Monitor System. The subroutine BCREAD (called from 
main subroutine VULCAN) may be used to read fluxes from binary cards. The subroutine 
PLOTXY forms plots of two related variables. The following printer carriage control 
characters have been used as the first character in a format statement: 

J Single space before printing rest of line 
K Double space before printing rest of line 
L Triple space before printing rest of line 
$ Punch the line instead of print 
* Print and punch this line 

For the first three control characters, the character must be counted as one of those 
being printed, but for the last two, the format is given as if the control characters were 
not a part of the format statement. The last control character * is actually followed by a 
blank to make up the complete control symbol. If these characters are not available in 
the monitor system used, then other output control statements must be used. 



SIBFTC VULCAN DECK 
C 

C VULCAN IS A POINT (MESH INTERVAL) DEPLETION CODE 

C IT HAS OPTIONAL INPUT FORMATS THAT MAKE IT ADAPTABLE TO 

C AN ARBITRARY ONE-DIMENSIONAL SPATIAL CALCULATION 

C 

C VULCAN USES NO PERIPHERAL STORAGE 

C 

C WRITTEN IN FORTRAN IV LANGUAGE 

C 

C JL ANDERSON NUCLEAR SYSTEMS DIVISION 

C LEWIS RESEARCH CENTER 

C AUGUST 9, 1966 

C 

COMMON/ SET 1/ ATDEN ( 190, AO) ,CONR AT ( 200 ) , TOPOW ( 200 ) , VOL ( 200 ) ,RR( 200) 
1 ,RATBEF ( 200 ) 

COMMON/SET 2/ FLUX ( 3800 ) , SANSNU ( AO , 20 > , ABS I G ( AO, 20 ) , SST ( 20 , 20 ) , 

1 XNU ( 12,20) 

COMMON/SET 3/ MM I R ( 200 ) , RAUM ( 200 ) , RAR ( 200 ) , DELR ( 200 ) ,MIDRZ( AOO ) , 
1NEM0I 200) 

COMMON/ S ETA /AT WGT ( 20 ) , DECAY (AO) ,THERNU(AO) , TNAME ( AO ) ,FRACT( 10) , 
lNRPT(AO) , IDVULCt AO) 

COMMON/SET 5/ VUL ABS I AO > ♦ VULF I S ( AO ) ,VULAMF(AO) , ADHREV ( AO ) ,IDYLD( 10) 
1 ,Y(20,12), KKUN ( 20 ) 

COMMON/S ET6/ NG, NM, NI SOT ♦ TOTVOL , F I WATT , POWER, AVHQW , FNORM 

COMMON/ SE T7/NFIS»NSPMAT »NBPUI ,NFRACT , NYDNUC , NF I KPT 

COMM0N/SETB/ DELHR , KHA I N , A VOGAD , NZONR , MOROUT , MOOD , D I FF ER , NONDPL 

COMMON/ SE 19/ KT I M , NT I NC , KR6LL ,KC ELL 

COMMON/ SET 10/ KOMP , KONG , KOG , LAND , DEL SEC , ME SH 

COMMON/ SET 11/ NSSF , JSMAT ( 20 ) 

C 

C 

DIMENSION COR POW ( 20 ) , RAPOF ( 20 ) , T I TL E 11 A ) , YELD (10 ) , I0C0N1 AO), 

1NUM ( 1A) ,KFI S< AO) , D E N S ( AO ) , IDFP (20 ) , P ( 1 5 ) , UNAME ( 10 ) 

DIMENSION FUNS I G ( AO , 20 ) 

DIMENSION FMT112) 

C 

EQUIVALENCE l FUNS I G , SANSNU ) 

C 

C 

C READ IN INPUT DATA 

C 

1A WRI TE (6, 107) 

2 R E AD( 5 , 100 ) I JK , ( T I T L E ( I ) , I = 1 , 1 A ) 

WR I TE ( 6 , 1 59 ) (TITLE (I), 1=1, 1A) 

IF< I JK.NE. 1 ) GO TO 2 
WRIT6(6,U0) 

READ l 5, 101 ) NG»NFG»NISOT »NMAT »NT I NC , NFR AC T , KC E L L 
WRITE (6, 101 ) NG,NFG,NISUT, NMAT , NT INC , NFRACT , KCELL 
WRI TEC 6, 111 ) 

RE AD ( 5, 101 > NZONK, NRM, MOROUT, NUF 
WR I TE ( 6, 101 ) NZONK, NRM, MOROUT ,NUF 
WR I TE ( 6 , 1 1 2 ) 

RE AD ( 5, 101 ) NS PM AT , NBPO I , NF I S , NYDNUC , NONDPL »NF IKP i 
WRITE (6, 101 ) NSPMAT, NBPO I ,NF IS, NYDNUC, NONDPL, NF I RPT 
WRlTE(6,113) 

READI5, 101 ) KG60»KFLUX»KSSF,KHAIN,KFAST,KF0RM 
WRITE (6, 101) KGEO,KFLUX,KSSF,KHAIN»KFAST,KFORM 
W R I T E ( 6 , 11A) 

RE AD (6, 102) POWER »DI ST 1 , D I ST 2 , DELHR T XK6M P , REM I OD , D I FF EK 
WRITE (6, 1A1 ) POWER, DI ST 1,DIST2, DELHR, XREMP , REM I OU , DIFFER 
IF(NFRACT.EO.O) GO TO 3 
RE AD ( 5 , 1 OA ) ( FRACT(JK) ,JK=1, NFRACT ) 

3 WR I T E ( 6 , 1 1 5 ) 

CALL TABLE ( N I SOT, F I WATT , A VOGAD ,K FAS T ) 

DECAY ( 1A)=0ECAY( 1A)+XR6MP 
DECAY < AO )=DECAY( 16) +REMIOD 

C DECAY ( AO ) IS ACTUALLY DEC AY ( 1 6 ) FOR THE 1-136 EQUATION 

DECAY ( 16) =DECAY( 1*6) -REMIOD 
C DECAY (16) IS FOR THE XE-135 EQUATION 

DO 1 K I SS= 1 »N I SOT 

C I REFERS TO MACGG SEQUENCE BUT CARDS (NISOT) MAY BE IN ANY ORDER 

RE AD ( 5 , 1 A3 ) I, IDVULC( I ) , NRPT ( I ) , KF I S ( I ) , DE SCR 
I SM= I DVULC ( I ) 

TNAME < ISM) =DESCR 
I F ( KF I S ( I ) . EQ . 0 ) GO TO A 

WRI T E ( 6 , 1 A 5 ) I, 1 DVULC ( I ) , NRPT I I ) , TNAME < ISM) 

GO TO 1 

A WR I TF ( 6 , 1AA ) I , I DVULC ( I ) , NRP T ( I ) , TNAME ( I SM ) 

1 CONTINUE 

RE AD (5, 103) (NMIR( I ) » RAR ( I) , Ml DRZ ( 1 ), 1=1, NZONR) 

WRITE(6,116) 

WRITE(6, 119) 

I GO=NRM+ 1 

00 21 1=1, NZONK 

1 END= I GO+NM IR ( I )-l 
00 22 J= I GO , I END 

22 MIDRZt J )=MIDRZ< I ) 

21 WR ITE(6,138) NM I R ( I ) , R AR ( I ) , ( M I DRZ ( K ) , K = I GO , I END ) 

C 

C MESH INTERVAL MANIPULATION 

C 

11 P I = 3« 1A15927 
OELR ( 1 )=0.0 

DE LR ( 1 ) =RAR ( 1 ) 

DO 12 1=2, NZONR 
DE LR ( I )=0.0 

12 DELRl I )=RAR( I )-RAR( 1-1 ) 

DO 13. 1 = 1, NZONR 

VM I =0 • 0 
VMI=NMIR( I ) 
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13 DELRU )=DELRU )/VMI 
RUN=0.0 
RR ( 1 )=0.0 
K-l 

DO 17 J = 1 , NZONR 
NMES=NHIR< J ) 

DO 17 1=1,NMES 
K= K+ 1 
RR t K ) = 0 • 0 

17 RR ( K ) =RR ( K-l ) +DELK ( J ) 

RUN=RR(NRM+1 ) 

9 CONTINUE 
C 

C VOLUME CALCULATION 

C 

NM=NRM 

DO 20 JK=1,NM 
20 VOL ( JX ) =0. 0 
NRMP 1 =NRM+ 1 
GO TO t 30 , 40 , 50 ) , KGEO 

30 V0LT=DIST1*DIST2*RUN 
RAUM 1 )=RAR C 1 ) -D I ST2*D I ST 1 
DO 23 N I L = 2 , NZONR 

RAUM ( N I L ) =0 • 0 

23 RAUM(NIL)=(RAK(NIL)-RAR( NIL-1 ) ) *D I ST 2*D I ST 1 
DO 31 K = 2 ? NRMP 1 

31 VOL ( K-l ) = ( RR ( K ) -RR I K- 1 > ) *D I ST2*D I ST1 
GO TO 60 

40 V0LT=PI*RUN*KUN*D1ST1 

RAUM ( 1 )=pI*RAR( 1 )*KAR( 1 )#0IST1 
DO 35 M I L = 2 » NZONR 
RAUM ( M I L ) =0 . 0 

35 RAUMIMI L)=HI^IKAR( MIL )^RAR( MIL )-RARt MIL-1 ) *R AR < M I L“1 ) ) *DI ST 1 
DO 41 K-2tNRMPl 

41 VOL (K-l )=PI*(RR( K)*RR(K )-RR(K-l )*RR(K-1 ) )*DIST1 
GO TO 60 

50 V0LT=l4.*PI/3.)*(RUh**3) 

RAUM ( 1 )=(4.*PI/3. )*(RAR( 1 )#*3) 

DO 57 J I L = 2 , NZONR 
RAUM ( J I L >=0.0 

57 RAUM( Jl L)=(4.*PI/3.)*(RAR( JIL )**3-RAR( JIL-1)**3) 

DO 51 K=2tNRMPl 

51 VOL (K-l ) = (4./3. )*PIM RR(K )**3-RR(K-l )**3) 

60 TOTVOL =0 • 0 

DO 52 LBJ=1»NM 

52 TOTVOL=TOTVOL+VOL(LBJ) 

WK I T fc < 6 , 1 17 ) 

WRITE (6, 139) ( VOL(K) ,K=1 ,NM) , TOT VOL 

WR I TE ( 6, 150 > ( RAUM ( KUK ) »KUK = 1 » NZONR ) 

QUO T = VOL T / TOT VOL 

IFIABSI l.-QUOT) -LT..001 ) GO TO 215 
WRITEI6,133> TOTVOL » VOLT 
C 

C READ IN MORE INPUT DATA 

C 

C FISSION INFORMATION 

C 

215 WKI TE ( 6, 107 ) 

DO 74 I = 1 » NYDNUC 
IF(NYDNUC-EO.O) GO TO 230 
IF(NFIS.GT.6) GO TO 5 

READ (5*131 J KRONin , IDFP( I ) » ( YELD1J ) , J=1 tNFIS) 

GO TO 6 

5 READ ( 5 » 1 3 1 ) KKON ( I ) T I DFP ( I ) , < Y EL D ( J ) , J = 1 , 6 ) 

READ< 5, 104) (YELD(K),K=7,NFIS> 

6 I JK = I DF P ( I ) - 1 2 

C I DFP IS IN VULCAN SEQUENCE 

230 DO 70 J=1,NFIS 
I F I S=0 

DO 73 ISEQ=1tNIS0T 
J F I S = I F I S + KF I S ( ISEQ) 

I F ( IFIS-NE.J) GO TO 73 
I DY LD ( J ) = I DVULC ( ISEQ) 

IDT=IDYLD( J) 

IF(NYDNUC.NE.O) Y(IJK,IDT)=YELD(J) 

QNAME ( J ) =TNAME ( IOT) 

GO TO 70 
73 CONTINUE 
70 CONTI NUF 
/ 4 CUNTINUE 

75 WRITE (6, 124) l I DY L 0 ( J A ) , J A= 1 , NF I S ) 

WRITE (6, 147) ( QNAME < JA) , JA = 1 ,NFIS ) 

NUG = NSPMA T + 17 
DO 55 NASA = l * NI SOT 
NAV=I DVULC( NASA ) 

IF (NAV.LE. 12.0R.NAV.GT . NUG ) GO TO 55 
NUT=N AV— 1 2 
DO 56 JET=ltNFIS 
ID=IDYLD( JET ) 

56 YELDt JET )=Y(NUT, ID) 

WRITE (6,119) 

WR I T E ( 6 » 1 46 ) NAV»TNAME( NAV) » (YELD(JJ ) ,JJ = 1,NFIS) 

55 CONTINUE 

IF(NUF.EQ.l) GO TO 42 
WRITE(6, 109) 

C 

C NU IN VULCAN IU ORDER 

C 

DO 79 J= 1 » 12 
DO 221 MIG=1 »NFG 
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221 XNU ( J , M I G ) =0 ■ 0 
NUM ( J ) = J 

00 220 NAX = 1 , N I SOT 

IF ( J . EO. I DVULC I NAX ) ) READ(5,104) ( XNU ( J , I G ) , I G= 1 , NFG ) 

220 CONTINUE 
79 CONTINUE 

WRITE(6»125) { NUM ( J ) » J = 1 » 12 ) 

WR I TE ( 6 » 1 1 9 ) 

DO 76 I G= 1 » NFG 

76 WR I TE I 6 , 14-0 ) I G, { XNU( J , I G ) , J=1 , 12 > 

IFIKFAST.EU.O) WRITE<6,142) ( T HERNU ( J ) , J= 1 » 1 2 ) 

42 WR ITEI6»1Q7) 

C 

C CROSS SECTIONS 

C 

IFIKF0RM.EW.1.0R.KF0RM.EW.3) READ(5,135) ( FMT I I ) » I = 1 , 1 2 > 
DO 67 I SQT = 1 » N I SOT 

C MACGG SEQUENCE 

I D= I DVULC < ISOT) 

K I 0 = I D 

READ! 5» 135) DUMMY 

IFIKF0RM.E0.0.0R.KF0RM.EW.2) GO TO 200 
DO 201 I G= 1 * NG 

201 RE AD ( 5 » FMT ) ABSIGI I D» IG) , FUNS I G ( IDt IG) 

GO TO 202 
200 DO 65 I G = 1 » NG 

READ! 5 1 105) - ABSIGI 1 D, I G > , FUNS I G ( 10, IG) 

READ(5,135) DUMMY 
65 CONTINUE 
’02 WR I TE ( 6 » 1 Ob ) 

WRITE (6,121) ISOT, 10 
WR I TE ( 6, 122 ) 

WRITE (6, 139) I ABS IGI ID, IG) , IG=1,NG) 

I F ( ID.GT. 12.AND.NRPTI ISOT) .EQ.O) GO TO 67 

IFlNUF.EQ.l) GO TO 300 

WRITE(6,123> 

WR I TE ( 6 , 1 39 ) ( FUNS 1 G ( I D , I G ) , I G= 1 , NG ) 

300 GNU=1. 

DO 68 I G = 1 » NG 

IFINUF.EO. 1 ) GO TO 68 

IFIFUNSIGI ID,IG) .EO.O.) GO TO 68 

1 F I NRPT I ISOT) • EW* 0 ) GO TO 69 
K I D = NRP T l I SOT ) — 1 0 

69 CONTINUE 

I F ( IG.GT.NFG) GO TO 72 
GNU = XNU (KID, IG) 

I F (GNU.NE.0- ) GO TO 68 
GNU = 2 . 5 

WRITEI6,153) I D , I G 
GO TO 68 

72 GNU=THERNUIKIDJ 

IF IGNU.EO.O. ) GNU= 1 • 

6 8 SANSNUI ID, IG)=FUNSIG< ID, IG)/GNU 
WRI TE(6, 134) 

WRITE (6, 1 39) I SANSNUI I D , I G ) , I G-l , NG ) 

67 CONTINUE 
C 

C ATOM DENSITIES 

C 

LUNE =0 
LUX = 0 
KANT=0 

IF IKCELL.LE.20) GD TO 207 
KANT=1 

KCELL=KCELL-20 
207 DO 208 N$EL = 1 » KCELL 
KTIM=1 
LOX =LOX+ 1 
KRELL=NSEL 

IF INSEL.GT. 1 ) WRITE(6,157) NS6L, KCELL 

61 DO 6 2 J = 1 ,NM^AT 
DO 64 K= 1 , 4 0 

64 DENS(K>=0.0 

READ! 5, 101) NC ON , ( I DCONl I ) , 1 = 1 , NCON) 

C MACGG SEQUENCE I DCON 

READ!5,15b> ( DENS < I ) , I = 1 , NCON ) 

L ADD=0 

DO 53 K= 1 , NZ ONR 

ME SH=0 

MAT =M I DR Z I K ) 

LADD=LADD+NMIR ( K ) 

IFIMAT.NE.J) GO TD 53 
MESH=LADD-NM IRIK) 

NMES=NMIR I K ) 

DO 53 L= 1 , NMES 
ME SH = ME SH + 1 
DO 235 NE 0= 1 , 40 
235 ATDEN ( MESH, NEW ) =0.0 
DO 54 KKK= 1 , NCON 
I AGO= I OCON ( KKK ) 

I D= I DVULC ( I AGO) 

54 ATDENIMESH, ID) =DENS ( KKK) 

53 CONTINUE 

62 CONTINUE 
NG0T0=0 
NST=0 

WRI TE(6,107) 

MUSS=0 
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36 NST=NG0T0+1 

IFtNMAT.GT. (NGOTO+7) ) GO TO 37 

NGOTO=NMAT 

GO TO 45 

37 NGOTO-NGOTO+7 
WRITE(6,119) 

45 DO 38 KIX=1,7 
NUM(KIX )=0 
NUM(K1X )=NST+KIX-1 
IF(NUM(KIX) .GT.NMAT ) GO TO 330 
KRAD=0 

DO 310 JKT = 1 , NZONK 
KRAD=KRAD+NMIR( JKT ) 

IF(MIDRZ( JKT).EO.NUM(KIX) ) GO TO 320 
310 CONTINUE 

320 NUM(KIX+7)=KRAD-NMIR( JKTJ+l 
GO TO 38 
330 NUM ( K I X ) =0 

38 CONTINUE 
LARAY=NG0T0-NST+1 
IF(NG0T0.GT.7) WRITE(6,119) 

WR I TE (6,118) ( NUM( K ) »K= 1 » 7 ) 

DO 66 MV= 1 *N I SOT 
JV=IDVULC(MV) 

WR I TE ( 6 » 1 19 ) 

DO 340 NVT=1,LAKAY 
LOCUS- NUM ( NVT+7 ) 

340 CONKAT(NVT ) -ATDEN ( LOCUS r JV) 

WRITE! 6,120) JV,TNAMt( JV) , ( CONRAT ( J ) , J=1 » LARAY ) 

66 CONTINUE 

IFtNMAT.GT. NGOTO) GO TO 36 
WR l TE ( 6 » 107 ) 

C 

C FLUX 

C 

LST=0 

LSP=0 

IFCNSEL.GT.l. AND.KANT.EO.il GO TO 88 
I ST = 0 
I SP = 0 
KST = 0 
KSP = 0 

204 IF (KFORM.LT. 2) GO TO (81,82), KFLUX 
READ(5,100) (FMT( I ) , 1=1, 12) 

DO 203 I G=1 , NG 
I ST= I SP+ 1 
I SP= I ST+NM-1 

203 READ ( 5 , FMT ) ( FLUX ( I CON > , ICON= I ST , I SP ) 

GO TO 87 

81 DO 77 I G= 1 » NG 
I GO^NM* ( 1 G-l )+l 
I ST0P=IG0+NM-1 

CALL BCREADt FLUX( I GO) ,FLUX( I STOP) ) 

77 CONTINUE 
GO TO 87 

8? DO 86 I G= 1 , NG 
IST=ISP+1 
ISP=IST+NH-1 

86 READ (5, 102) (FLUX( ICON) , ICDN= I ST t I SP ) 

87 CONTINUE 

IF (LUNF.EU.O) WKl T E ( 6 , 1 ?6 ) 

IF(KANT.EU.l) WR I Tt ( 6, 1 29 ) 

lFIKANT.EU.O.AND.KCELL.GT.O) WRITE(6,130) NSEL 

DO 80 I G= 1 , NG 

WkITF(6,127) IG 

KST=KSP+1 

KSP=KST +NM-1 

80 WR I TE ( 6 ♦ 1 39 ) ( FLUX ( KCON ) , KCON=KST , KSP ) 

88 CALL ABSPHI 
KOOL=0 

IF(KANT.EU.O) GO TO 206 
IF(NSEL.GT.l) GO TO 212 
LUNE=LUNE+1 
READ(5,106) NKOREG 

C ACTUAL tFRACT IONAL OR RELATIVE POWER MAY BE READ IN 

READ (5, 102) ( CURPOW (KIP) , K IP= 1 , NKOREG ) 

WRITE (6, 154) ( CORPOW (KIP) ,K I P=1 , NKOREG ) 

DO 214 KRAT=1 , NKOREG 
214 RAPOF ( KRAT ) =CORPOW ( KRAT ) / COR PUP ( 1 ) 

WRI TE (6, 128) (RAPOF ( JOE) , J 0E= 1 , NKOREG ) 

212 MG=NMvNG 
V I E T=i . 0 

IF(NSEL.GT.l) V1ET=RAPUF(L0X-1» 

DO 205 KK- 1 , MG 

205 FLUXIKK )=FLUX(KK)*KAPOF(LOX)/VIET 
WR I TE ( 6 , 160 ) RAPOF ( LOX ) 

WRI TE (6, 126) 

WR I TE ( 6 , 1 30 ) NSEL 
DO 43 I GY = 1 , NG 
WRI TE ( 6, 127 ) IGY 
LST'=LSP + 1 
LSP=LST+NM-1 

43 WR ITE(6,139) ( FLUX < NAM ), NAM=L$T , LSP ) 

206 DO 209 NTM- 1 , NT INC 

IF (KSSF.GT .0 ) CALL SSF ( NG , NTM ) 

KT I M-NTM+1 
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89 CALL BURNUP 

IF (NTM.NE.NTINC ) GO TO 90 
99 CALL OUTPUT 
C 

C CALCULATION OF DIMENSION SPECIFICATIONS 

C 

NET=1 

WR I TE ( 6 » 148 ) 

DO 92 K= 1 » MOOD 
NET=NEMO(K )+NET 
RAR(K>=RR I NET ) 

PUNCH 103, NEMO(K) tRAR(K) ,K 
92 WR I TE ( 6 , 149 ) NEMO ( K ) , RAR ( K > ,K 

90 PI1)=NH 

CALL PLOTXYlRRfCONKAT » 64, P ) 

WRITEI6, 136) 

CALL PL0TXY(RR,T0PUW,64,P> 

WRITE(6, 137) 

WR I TE ( 6 , 1 07 ) 

WR I TE ( 6 » 156 ) NTM » NT I NC 
209 CONTINUE 
208 CONTINUE 

91 GO TO 14 

100 FORMAT I I 1, IX, 14A5) 
in l FORMAT (71 10 ) 

102 FORMA I ( ZE 1 0. 5 ) 

103 FORMATt I5,F10.6, I 5) 

104 F0RMAT(7F10.8) 

105 FORMAT (36X,2E12»6) 

106 FORMAT ( I 1Q,6F10.7) 

107 FORMAT ( 1H 1 ) 

108 FORMA T ( 1HK ) 

109 FORMAT ( 1HL ) 

1 10 FORMAT ( 1 HL » 7X » 2HNG, 8X, 3HNFG, 6X , 5HN I SOT , 5X , 4HNMAT , 6X , 5HNT I NC , 4X ♦ 
16HNFRACT » 4X » 5HKCELL ) 

1 1 1 FORMAT ( 1 HL » 6X , 5HNZ0NK , 5X , 3HNKM , 6X , 6HM0R0UT , 6X , 3HNUF ) 

112 FORMAT ( 1HL, 5X,6HNSPMAT , 5X, 5HNBP0I , 5X ,4HNF I S » 5X , 6HNYDNUC » 4X , 6HN0NDP 
1L,4X,6HNFIRPT) 

113 FORMAT ( IHLt 6X ,4HKGtO»6X ,5HKFLUX»5X,4HKSSF ,6X ,5HKHAI N, 5X , 5HKFAST , 

1 5X » 5HKF0RM ) 

114 FORMAT < 1HL,9X,5HP0WER»6X,5H FIRST , 8X , 6HSEC0NQ ,8X ,4HT I ME , 1 IX , 5HXEN0N 
1 , 1 OX , 6H I DDINF/1 IX,2MIN, 7X»9HDI STANCE , 4X , 8H0I S T ANC E , 4X , 1 OH INCREME 
2NT ,7X, 7HREM0VAL »9X, 7HREMQVAL , 5X , 6HD I FFER/ 10X ,5HWATTS ,6X , 4H ( CM ) , 

3 1 OX » 6H ( CM ) ,5X, 

4 8H (HOURS) ,4X, 13H PROB AB I L I T Y , 5X , 1 1 HPKO B A8 I L I T Y ) 

115 FORMAT ( 1HL»4X» 10HSE0UENTIAL»7H VULCAN, 8H NRPT ,/HISOTOPE/ 
18X,2HID,8X,2HID»11X,5HLAB£L) 

116 FORMAT ( 1HL * 12H NMIH(I)«12H RAR(I),23H MIORZ(I)- 

1 MAP ) 

117 FORMAT ( 1HL , 15H MESH VOLUMES//) 

118 FORMAT ( 8H VULCAN, 11H ISOTOPE ,7(13H MATERIAL )/4X»2HID,5X, 

1 5H L AB E L » BX , 7 ( 12, 1 1 X ) ) 

119 FORMAT(IHJ) 

120 FORMAT! I5»6X,A5,3X,7E13.5) 

121 FORMAT { 1 7H SEQUENTIAL ID = , 12 ,15H VULCAN ID = ,12) 

122 FORMAT ( 1HK»34H ABSORPTION BY GROUP (HIGH TO LOW)) 

123 FORMAT ( 34H NU*FISSION BY GROUP (HIGH TO LOW)) 

124 FORMAT ( 1 OH F I SS I ON , 6X , 10 ( 1 OHY I ELD )/15H PRODUCT ,10( 

1 4HFR0M » I 3 , 3X ) ) 

125 FORMAT ( 7H GROUP, 12(8H NU FOR ) / 7X , 1 2 ( 5H VI0,I3)) 

126 FORMAT ( 3 1H FLUXES BY GROUP AND MESH POINT) 

127 FORMAT ( 1H J , 8H GROUP ,12) 

128 FORMAT ( 1 HK , 35H THESE ARE THE RADIAL POWER F AC TORS/ < 8 E 14 . 6 ) ) 

129 FORMAT ( 24H FUR THE BASIC CELL) 

130 FORMAT ( 22H FOR CELL NUMBER ,12) 

131 FORMAT (2I5,6FL0.8) 

133 FORMAT ( 1HL , 1 3H VOLUME SUM , F 1 3 . 4 , 2 5H NOT E U U A L F ( 0 1 ST 1 , D I ST 2 ) , F 1 3 
1,4) 

134 FORMAT ( 33H FISSION BY GROUP (HIGH TO LOW)) 

135 FORMAT ( 1 2A6 ) 

136 FORMAT (2HPL,55X,41HCONVERSlON RATIO VS. MESH INTERVAL RADIUS) 

137 FOKMAT(2HPL»55X»53H NORMALIZED POWER IN INTERVAL VS. MESH INTERVAL 
1 RADI US ) 

138 FORMAT ( I 10,7X,F8.4, 9X , 2 5 ( 1 3 ) / 34X , 2 5 ( I 3 ) ) 

139 F0RMAT(8E14.6) 

140 FORMAT 1 1 5 , 2X , 1 2F8. 3 ) 

141 FORMAT ( 7X » E 1 0 • 5 » F 1 1 • 5 » F 1 2 • 5 » F 1 4. 2 ♦ F 1 6 • 7 » F 16 • 7 , F 1 1 • 5 ) 

142 FORMAT ( 7H T ,12F8.3) 

143 FORMAT ( 4 1 5 , A5 ) 

144 F0RMAT(2I 10, I7,4X,A5) 

145 FORMAT! 21 10, 2H Y,I5«4X,A5> 

146 FORMAT ( I 5, 2H ( , A5 , 1 H ) , F9 . 5 , 9 ( FI 0 . 5 ) ) 

147 FORMAT (15X,10(1H(,A5,1H),3X)) 

148 FORMAT ( 1 HL , 41 H DIMENSION SPECIFICATION INfUT TO T DSN// ) 

149 FORMAT ( I 1 5 , F10.6 » 15) 

150 FORMAT! 1HL» 15H ZONE VOLUMES// ( 8E 14.6 > ) 

152 FORMAT ( F 1 5 • 4 ) 

153 FORMATClHLtlO! 1H+) ,61H A NON-ZERO FISSION CROSS SECTION HAS BEEN 
iSUPPLIED FOR VID , I 2/ 1 X , 10 ( 1 H+ ) , 8H GROUP ,I2,54H BUT NO NON-ZERO 
2NU VALUE IS AVAILABLE. THEREFORE THIS/1X, 10( IH+) ,41H PARTICULAR N 
3U HAS BEEN SET EQUAL TO 2.5) 

154 FORMATt 1HK , 3 1 H THE POWER FOR EACH CORE REG ION/ ( 8E 14 . 6 ) ) 

156 FORMATt 1H2, 1 0 1 1H*) ,31H THIS COMPLETES TIME INTERVAL ,I2»4H OF ,12 
1 ) 

157 FORMAT ( 1HL, 10 ( 1H* ) »62H THE FOLLOWING ATOM DENSITIES AND FLUXES AR 
IE FOR CELL NUMBER , I2,4H OF ,12) 

158 FORMAT t 7E 1 0 • 6 ) 

159 FORMAT ( 2X , 1 4A 5 ) 

160 FORMAT ( 1HK , 38H THE POWER FACTOR FOR THIS CELL IS ,F10.7> 

END 
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ilBFTC TTABLE DECK 

SUBROUTINE TABLE t NI SOT , F I WATT , AVOGADt KFAST ) 

COMMON / SET 1/ AT DENI 190,40) ,CONRAT ( 200 > , TOPOW ( 200 ) , VOL < 200 ) , RR ( 200 ) 
1 , R A TBEF ( 200 ) 

COMMON/ SET2/ FLUX ( 3800 ) , SANSNU ( 40 , 20 > , ABS IG ( 40 , 20 > , SST ( 20 , 20 ) , 

I XNU l 12 t 20 ) 

COMMON/ SET3/ NMI R ( 200 ) , RAUM I 200 > ,RAR(200) , DELR t 200 ) ,MIDRZ ( 400 > , 
INEMO ( 200 ) 

COMMON / S6T4/ A TWGT 120) » DECAY ( 40 ) » THERNU ( 40 ) ,TNAME ( 40 ) , FR ACT ( 10) * 
1NRPT (40) , I DVIJLC ( 40 ) 

COMMON/ SETS/ VULABSI 40) ,VULF!S<40 ) ,VULAMF (40) , ADPREV < 40 ) , I DYLO ( 10) 
I ,Y(20,12)»KRQN(20) 

C 

C 

C SUPPLIES MOST COMMONLY USED TABULAR INFORMATION 

C 

c 

c VULCAN 10 NUCLIDE 
C 1 TH232 

C 2 PA233 

C 3 U233 

C 4 U234 

C 5 U235 

C 6 U236 

C 7 U238 

C 8 NP239 

C 9 PU239 

C 10 PU240 

C 11 PU241 

C 12 PU242 

C 13 FISSION PRODUCT AGGREGATE (OR INDIVIDUAL PRODUCT) 

C 14 X6135 

C 15 SM149 

C 16 1135 

C 17 PM 1 49 

C 18 SPECIAL MATERIALS INCLUDING ADDITIONAL F.P. AGGREGATES 

c * • 

C 

C 

c . BURNABLE poisons 

c • • 

c • • 

C • • 

C . REPEATED FISSIONABLE NUCLIDES 

C • • 

c • • 

c • • 

c . non-depletable nuclides 

c • • 

c • » 

c 40 

c 

c 

DO 10 JFK=1,20 
DO 10 LB J= 1 » 1 2 
10 Y ( JFK » LB J ) =0. 0 

C Y(PROOUCT »SOURCE) = FR ACTIONAL YIELD 

C EACH Y ( 1 » X ) IS (l.-SUM OF INDIVIDUAL YIELDS) 

C PRODUCT NUMBER 2 = XE-135 

C PRODUCT NUMBER 3 = SM-149 

C PRODUCT NUMBER 4 = 1-135 

C PRODUCT NUMBER 5 = PM-149 

C Y ( 1 » 3 ) = .9481 

C Y ( 1 » 5 ) =.9246 

C Y(l,9>=.9113 

C Yd.ll )=.9220 


Y(2t3)=.0020 

GUESS 

Y ( 2 * 5 ) =• 003 1 

GAFKAT C 

Y ( 2 , 9 ) = .001 7 

GAFKAT C 

Y ( 2 , 1 1 ) =.0020 

GUESS 

Y(3,3)=.0012 

GaROOS C 

Y(3»5)=. 00 1 3 

GAROOS C 

Y(3»9 )=.007Q 

GAROOS C 

Y ( 3 » 1 1 ) - « 006 

GUESS 

Y(4, 3)^.0426 

NSE25 TN 

Y ( 4 » 5 ) = ♦ 06 1 0 

GENERAL 

Y ( 4 » 9 ) = *0700 

GAFKAT 

Y ( 4 », 1 1 )=.0600 

GUESS 

Y ( 5 » 3 ) = .0050 

GAROOS 

Y ( 5 » 5 ) =• 0100 

GAROOS 

Y ( 5 » 9 ) = * 0100 

GAROOS 

Y(5, 11 >=.0100 
IF(KFAST.EO.O) GO TO 30 

GUESS 

20 Y ( 2, 5 )=• 0025 

GAFKAT C 

Y ( 2 » 9 ) “ • 0090 

GAFKAT C 

Y ( 3 » 5 ) = • 00 1 5 

GAFKAT C 

Y ( 3» 9 ) = .0035 

GAFKAT C 

Y ( 4 , 5 ) = • 0600 

GAFKAT 

Y (4,9)=. 0600 

GAFKAT 

Y ( 5 , 5 ) = • 0 100 

GAFKAT 

Y t 5 ,9 ) = *0 100 

GAFKAT 
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C Y(l,5)=.9260 

C Y<1»9)=.9175 

30 ATWGT ( 1 ) =232* 111 
ATWGT(2)=233.114 
ATWGT t 3 > = 2 33. 1 12 
AlWGU4)=234.114 
ATWGT (5 ) = 235. 117 
ATWGT(6)=236. 120 
ATWGT ( 7 >=238. 125 
ATWGT(8)=239.123 
ATWGT(9)=239.127 
ATWGT l 101=240. 129 
ATWGT ( 11 >=241. 131 
A T WGT ( 12 1=242.134 
DO 90 LBJ=1 tNIS0T 
DECAY ( LB J > = 0 • 0 
90 THERNU(LBJ>=0.0 


c 

THERNU (11=2-6 
THERNU<3)=2.503 
THERNU l 3)=2. 494 

ATOMIC 

ENERGY 

REV) 

[EW 

VOL 

3 =2 

( I960) 

ANL5B00 

GATHER 

c 

THERNU ( 5 ) =2. 43 
THERNUm=2-60 
THERNU!9)=2.882 
THERNU!9)=2.871 

ATOMIC 

ENERGY 

REV] 

[ E W 

VOL 

3 =2 

! I960) 

ANLB800 

c 

THERNU! 1 1 ) =3-98 
THERNU! 1 1 >=2.969 

ATOMIC 

ENERGY 

REV] 

EEW 

VOL 

3 =2 

! 1960 > 



FIWATT=3.10 ElO 
AV0GAD=6.0247 E23 
DECA Y ( 2 ) =2 • 9 28 E-7 


DECAY l o ; = 3 • 4 1 U E-to 

NUD 

V 1N5 

DECAY 1111 = 1.630 E-9 

NUD 

V1N5 

DECAY ( 14) =2. 10 E-B 



DECAY! 16)=2. 87 E-5 



DECAY!17)=3.63 E-6 




RETURN 

END 



SIBFTC ABFLUX DECK 

SUBROUTINE ABSPHl 

COMMON/ SET l / ATDEN ( 190,40) , CONK AT ( 200 ) ,TOPOW( 200 ) , VOL ( 200 I , KR { 200 ) 
1 » R ATBEF (200) 

COMMON/S ET2/ FLUX ( 3800 ) , SANSNU ( 40 , 20 ) , ABS I G ( 40 , 20 ) , SST < 20 , 20 ) , 

1 XNU ( 12,20) 

COMMON/SET 3/ NM I R ( 200 ) , RAUM ( 200 ) ,RAR(200) ,DELR(200) ,MIDRZ(400) , 

IN E MO ( 200 ) 

CQMM0N/SET4/ATWGT( 20) , DECAY (40) ,THERNU(40) , TNAME ( 40 ) , FRACT ( 10), 
1NKPT ( 40 ) , I DVULC ( 40 ) 

COMMON/ SET 5/ VUL ABS ( 40 ) , VUL F I S ( 40 ) , VUL AMF ( 4 0 ) , ADPKE V ( 40 ) , I OY LD ( 1 0 ) 
1 »Y(20»12)» KRON ( 20 ) 

COMMON/ SET6/ NG , NM , N I SOT , TOT VOL , F I WATT , POWE R , AVPUW , FNOKM 
COMMON/ SET')/ KT I M , NT I NC , KRELL , KC ELL 
C 

C CALCULATES ABSOLUTE FLUX NORMALIZATION FACTOR 

C 

107 FORMA T ( 1H 1 ) 

132 FORMATt 1HL,45H FOR THE START OF TIME INTERVAL NUMBER ,12, 

14H OF ,12) 

134 FQkM AT ( 1HK , 1 7H CELL NUMBER ,I2»4H OF ,12) 

135 FORMAT ( 1HL,29H FLUX NORMALIZATION FACTOR ,E12.6, 48H CONVERTS T 
10 ABSOLUTE FLUX IN NEUTRONS/BAR N^S EC /20X,22H FOR A POWER LEVEL OF 
2,E12.6,7H WATTS) 

136 FORMAT ( 1HL, 30 (3H***1 ) 

IF(KTIM.EO.l) WRITE(6,107) 

V0LPGR = 0 » 0 

F I SPH I =0 . 0 
DO 10 I G = 1 » NG 
FISPGR=0.0 
DO 9 I M= 1 , NM 
FISMAC=0*0 
LOC=NM*( I G— 1 ) + I M 
DO 8 I MAT = 1 , N I SOT 
I D= I DVULC ( IMAT) 

IF(SANSNU( ID, IG) . fcO.O. .OR.ATDENt IM, ID) .ECJ-O. ) GO TO 8 
FISMAC=FISMAC+ATDEN( IM, I D ) *S ANSNU ( I D , I G ) *FLUX { LOC ) 

I F ( VOL ( IM) .LT.O. ) GO TO 8 
VOL ( IM)=-VOL( IM) 

VOLPGR=VOLPGR-VOL( IM) 

8 CONTINUE 

I F ( VOL ( I M ) • GT • 0 * ) GO TO 9 
FISPGR=FISPGR-FI SMAC^VOL ( IM) 

9 CONTINUE 

FISPHI=FISPHI+FISPGR 

10 CONTINUE 
FISPHI=FISPHI/VOLPGR 
DO 6 K=L,NM 

IF { VOL (K ) .GE.O. ) GO TO 6 
VOL ( K ) =-VOL ( K ) 

6 CONTINUE 

FNORM=( FI WAT T* POWER/ ( VOLPGR*F I SPH I ) >* 1 0 .**( -24 ) 

WR I TE ( 6 , 1 36 ) 

WKITEI6,132) KT I M , NT I NC 
WRITE(6,135) FNORM, POWER 
IFIKCELL.GT.O) WRITE(6,134) KRELL, KCELL 
WR I TE ( 6 , 1 36 ) 

NGM=NG«NM 
DO 11 LOC= 1 , NGM 

11 FLUX ( LOC ) =FLUX ( LOC ) *FNORM 
AVPOW=POWEK/ VULPGK 
RETURN 

END 
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SIBFTC TXENON OECK 

SUBROUTINE XENON ( GA I NXE * GA I NSM , GA I N I U » GA I NPM , AVGX E , A VGSM, A VG 1 0, 
1AVGPM) 

COMMON/ SET 1 / ATDENf 190,40) ,CONRAT<200) ,T0P0W(200) , VOU 200 ) , RR ( 200 ) 
1 t KATBEF ( 200 ) 

COMMON/SET 2/ FLUX ( 3800) , SANSNU ( AO * 20 ) , ABS I G ( AO , 20 ) , SS T {>20 , 20 ) T 
1 XNU ( 1 2 » 20 ) 

COMMON/SET 3/ NMIR( 200 ) ,RAUM( 200 ) ,RAR (200 ) , DELK ( 200 ) , M I DK2 ( 400) , 
1NEM0 ( 200 ) 

COMMON/S ET4/ATWGT! 20) » DECAY! 40 ) t ThERNU(40 ) ,TNAME (40) ,FRACT ( 10) T 
1NRPT ( 40 ) » 1 DVULC ( 40 ) 

COMMON/SET 5/ VUL AB S ( 40 ) , VULFIS! 40) , VULAMF ( 401 , AUPKEVt 401 * I UYLUl 10 ) 
1 » Y ( 20 1 12 ) » KRON ( 20 ) 

COMMON/ SET6/ NG » NM , N I SOT t TOT VOL , F I WATT , POWER » AVPOW « FNORM 
C0MM0N/SET7/NFIS»NSPMAT,NBP0I ,NFRACT ,NYONUC»NFIRPT 
COMMON/ SET 8/ DELHK »KHA IN* A VOGAOtNZONR, mOKOUT, M00l)» 01 FFEK * NUNDPL 
COMMON/ SET 10/ KOMP, KONG, KOG, L AND , DELS EC , MESH 
C 

C CALCULATES EQUILIBRIUM POISON CONCENTRATION AND TIME AFTER 

C SHUTDOWN OF MAXIMUM XENON CONCENTRATION 

C 

101 FORMAT ( 7 3H EQUILIBRIUM CONCENTRATIONS OF MAJOR FISSION PRODUCTS 
1 I N NUCLEI/BARN V CM/33H IF THE CURRENT POWER LEVEL UF .E10.5.21H 
2WATTS IS MAINTAINED) 

102 FORMAT ( lHKt 71H TIME AFTER SHUTDOWN UNTIL MAXIMUM XENON CONCENTR 
1ATI0N IS REACHED = ,F6.2,7H HOURS) 

103 FORMAT ( 1HK, 14X, 14HXEN0N-135 .E20.6/ 1 5X , 1 4HS AMAR I UM- 1 4SJ f E20. 

16/ 1 5 X » 14H IODINE-135 ,E20.6/ 1 5X t 1 4HPR0METHIUM- 14V » E20.6 ) 

104 FORMAT ( 1 4h EUUILIBKIUM ,4E15.6> 

105 FORMAT ( 1HL » 5 1 H ++ IT SEEMS THE FLUX LEVEL ISN.T HIGH ENOUGH FOR 

144H ++ A VALID TIME AFTER SHUTDOWN CALCULATION) 

106 FORMAT ( 55H THE XENON-135 CONCENTRATION (PEAK) AT THIS TIME 

IIS » E 1 2 • 6 ) 

C 

C 

£01 l)D=GA INI 0/ (DECAY ( 1 6 ) + VUL A BS ( 1 6 ) ) 

£OXEN=<GAINXE+DECAY< l 6 > #EQ I OD ) / ( DECAY ( 14J+VULABSI 14) ) 

FQ PKO = GA I NPM/ ( DECAY ( 1 7)+VULABS (17)) 

EOSAM=(GAINSM+DECAY( 1 7 ) *E«PRO ) /VUL ABS (15) 

IFIrtOROUT.LT.2l GO TO 60 
IF(NFRACT.EO.O) GO TO 30 
WR I TE ( 6» 104 ) EUXEN, EUSAM, EQIOO, FQPRO 
GO TO 40 

30 WRI TE(6» 101 ) POWER 

WR 1 TE ( 6 » 1 03 ) EuXEN»EQSAM»EOIUD?EQPRO 
40 TEMP=l.-( DECAY! 14) -DECAY! 16) ) *EUXEN/ < DECAY ( 16)«E0I0D) 

IF ( TEMP.LT.O. > GO TO 10 

TMAX = ALOG( TEMP*DECAY( 14) /DECAY ( 16) )/ (DECAY! 14) -DECAY! 16)1/36 00. 

GO TO 20 
10 TMAX=-1 • 

20 WR I TE ( 6 ♦ 102 ) TMA X 

IF ( TMAX.LT.O. . AND. TMAX.NE.-l . ) WR I TE (6, 10 5 ) 

T AF SEC = 3600. 4TMAX 

CONMAX= (DECAY ( 1 6 ) * L- 0 1 OD/ ( DEC A Y ( 14)-DECAY(16)))*(EXP(-TAFSEC*DECAY( 
116))— EXP t-TAFSEC* DECAY (14)) ) + EUXEN*EXP ( -T AFS tC *DECAY ( 14) ) 

WR I TE ( 6. 106 ) CONMAX 
60 AVGXE = AVGX£ + EDXEN*VC)L(MESH) 

AVGSM=AVGSM+EUSAM*VOL (MESH) 

AVGIQ=AVGIO+EOIDD w VOL (MESH) 

AVGPM= AVGPM+EOPRO^ VOL (MESH) 

50 RETURN 
END 
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SIBFTC SEISM DECK 

SUBROUTINE SSP(NG»NTM) 

C 

C CALCULATES SSF AND ADJUSTS CROSS SECTIONS 

C 

COMMON/ SET1/ ATDENI 1 90 , 40 ) , CONR AT < 200 J , TOPOW ( 200 } , VOL ( 200 ) * KR ( 200 ) 
1 , R ATBEF ( 200 ) 

COMMON/ SET 2/ FLUX! 3800) , SANSNUI 40, 20 ) , ABS I G ( 40, 20 ) , SST ( 20, 20 ) , 

1 XNU ( 1 2 » 20 ) 

COMMON/ SET 3/ NMIK(200) ,RAUM<200) ,RAR ( 200 ) , D6LR ( 200 ) ,MIDK2<400), 
l NEMO (200 ) 

COMMON/S ET4/ATWGTI 20) ,OECAY( 40) , TH6RNU < 40 ) , TNAME < 40 ) , FRACT < 10 ) , 
1NRPH40), I DVULC ( 40 ) 

COMMON/SETS/ VUL ABS ( 40 ) , VULFISI 40) , VULAMF (40) ,A0RREV(40) » I DYLD( 10) 
1 , Y ( 20, 12 ) ,KRON( 20) 

COMMON/ S ET 1 1 / NSSF » JSMAT ( 20 ) 

100 FORMAT (7110) 

101 FORMAT (7F10«6) 

103 FORMAT ( 1HL»56H SELF-SB I ELD I NG FACTORS USED ARE (HIGH TO LOW ENE 
1 RGY ) // ) 

104 FORMAT ( 1HK, 14H VULCAN ID = , 1 2 r2 X , 9F 1 0 .6 / ( 1 5X , 9F 1 0 . 6 ) ) 
IF(NTM.GT.l) RETURN 

WRITE(6,103) 

READ ( 5 t i 00 ) NSSF 

READ (5,100) { JSMAT( I ) , 1 = 1 ,NSSF) 

200 DO 201 ISF*1,NSSF 

201 READ(5t 101) ( SST ( ISF,J) ,J = l ,NG) 

500 DO 502 NSF*l,NSSF 

I SM= JSMAT ( NSF ) 

DO 503 M I G= 1 , NG 

ABSIGl ISM,M1G)=ABSIGI I SM , M I G ) *SST ( NSF , M I G ) 

503 SAN SNU ( ISM,MIG)=SANSNU( I SM , M I G ) *SS T ( NSF , M I G ) 

WRITE (6, 104) JSMAT < NSF) , ( SST (NSF, J ) , J = l »NG ) 

502 CONTINUE 
RETURN 
END 
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SIBFTC BURNT DECK 

SUBROUTINE BURNUP 
C 
C 

C PERFORMS ACTUAL DEPLETION* POWER, CONVERSION RATIO, AND 

C FUEL INVENTORY CALCULATIONS 

C 
C 
C 

COMMON/ SET l / AT DENI 190,40 > ,CONRAT (200 ) ,TOPOW ( 200 ), VOL I 200 ) ,RR( 200) 
1*RATBEF( 200 ) 

CQMMON/SE T2/ FLUX! 3800) , SANSNU ( 40 , 20 ) , ABS IG I 40 , 20 > , SST I 20 , 20 > , 

1 XNU ( 12,20) 

COMMON/ SET3/ NMIRI 200) , RAUM ( 200 ) , KAK ( 200 ) , DELK ( 200 ) , M I DKZ ( 400 ) , 
1NEM0I200) 

COMMON/ SET4/ATWGT 120) , DECAY (40) ,THERNU(40) ,TNAME ( 40 ) , FRACT ( 10) , 

1 NRP T ( 40 ) , I DVULC ( 40 ) 

COMMON/ SET 5/ VUL ABS ( 40 ) , VUL F I S ( 40 ) , VUL AMF ( 40 ) , ADPKEV ( 40 > , I DYLD ( 10) 
1 ,Y(20, 12) ,KR0N(20) 

C0MM0N/SET6/ NG, NM, N I SOT ,TOT VOL , F I WATT , POWER , AVPOW , FNORM 

CGMM0N/SET7/NFIS,NSPMAT»NBP0I,NFRACT ,NYDNUC ,NF IRPT 

COMMON/ SET8/ DELHR , KHA I N , A VOGAD , NZONR ,MOROUT »MGOD » DIFFER, NONDPL 

COMMON/ SET9/ KT I M , NT I NC , KRELL ,KC ELL 

COMMON/ SET 10/ KOMP , KONG ,KOG , LAND , DEL SEC »MESH 

C0MM0N/SET1 1/ NSSF, JSMATI 20) 

DIMENSION FP I M ( 40 ) , FAD I N ( 40 ) , BEGMAS ( 1 2 ) , F INMAS < 1 2 ) .POISON (40) 
DIMENSION XNUMERt 10 ) , GAMSP ( 20 ) , TFRAP ( 40 ) , TAX < 40 ) 

DIMENSION PRFDENI 40) ,P I ZEN I 40 > ,VOLL (40) 

C 

C INITIALIZATION 

C 

A VGX E =0 • 0 
AVGSM=0.0 
AVG 10 = 0 * 0 
AVGPM=0.0 

DO 315 J AU= 1 , N I SOT 
TFR AP ( J AO J =0 . 0 
PREDEN( JAO>=0.0 
PIZEN(JA0)=0.0 
315 VOLL ( J AU ) =0 • 0 
DO 50 ME SH = 1 , NM 
RATBEF(MESH)=0.0 
CON RAT ( MESH) =0.0 
TOPQW( MFSH)=0.O 
50 CONTINUE 
M t S H = 1 
PETK0L=0-0 
I GLQ0=0 
P0UR=0.0 
DO 2 L-1,12 
BEGMAS ( L ) =0.0 

2 F INMAS I L ) =0.0 
BMAS5=0. 0 
FMASS=0.0 
TC0NSB=0.0 
TDESTB=0.0 
TRHEF=0.0 

MUF =0 

KOOL=KOnL+l 
neLStC = 36U0.*L>ELHK 
SECOND =DELSEC 
C 

C BEGINNING OF MtSH ITERA1IUN 

C 

1 CONTINUE 
KOiMG = U 
KUG = 1 
K U M P = 1 
L AND=0 
KLON-O 
J A L P = 0 
JUMb=0 

11 DO 13 I DU= 1 , N I SOT 
1 0=1 UVULC ( I DO) 

VULABS ( I D) =0.0 
VULF I S ( ID ) =0.0 
VUL AMF ( I D ) =0 . 0 
DO 14 I G = 1 , N(, 

LOC=NM* ( I G— 1 ) +MtSH 

VULABS ( I D ) =VULABS ( 1D)+ABS1G( I 0 ♦ 1 G ) ? FL U X ( LUC ) 

14 VULF 1 S ( I D) = VUL FI S( I UJ+FLUX ( LUC) -SANSNU ( I D, IG) 

13 VUL AMF ( IU)=VULABS( ID)-VULFIS< ID) 

OEN0M=O.O 
DU 3 10=1,12 

IF(ATDEN(MtSH, IDJ.NE.O.) GO TO 4 

3 CONTINUE 

DO 365 KUKAN=1,NISUT 
ID=IOVULC(KURAN) 

I F < ATDEN(MESH, ID) ,NE.O..AND.NRPT(KORAN) .NE.O) GO TU 4 
365 CONTINUE 

IFlNBHOI.EU.O) GO TU 3H0 

JALP=NSPMAT+1B 

JUME= JALP+NBPU I — 1 

D(J 381 KOCH= JALP, JUME 

IF(ATDEN(MESH»KOCH) . FO. 0 . ) GO TO 3B1 
C K LUN= 1 INDICATES THAT A BURNABLE POISUN WILL BE CONSIDERED 

C IN THE CURRENT MESH 

KLUN=1 
GO TU 35 
381 CONTINUE 
380 MSUM=0 

Dll 10 K AT = 1 » NZUNR 
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MSUM-HSUrt+PlM 1 K 1 KAT ) 

Ir(MSUM.E(J. (MESH-1 ) ) GU TU 3b 
10 CONTINUE 

WKITE<6,151) rtSUM* MESH- 
GO Tli 3b 
4 CONTINUE 

IF(NBPOI.EU.O) GU TU 19 
JAA=NSPMAT+lB 
JZZ=JAA+NBF0I-1 
DO 413 LOUK=JAA,JZZ 

IF< ATDENImESH, LOOK) .£0.0. ) GO TU 413 
KLON-1 
GO TO 19 
413 CONTINUE 

19 IFIKHAIN.GT..1 ) CALL FISEU2(i9> 

CALL FISEUllS>9) 

C 

C FISSION YIELDS (GAmMA) 

C 

9 DO 21 J0=1,NFIS 
ID=I DYLU ( JU ) 

IF(NFlKPT.EU.O.UR.ADPKEV( ID) .NE.O. ) GU TU 361 
DO 360 K AD= 1 » N I SUT 

IF( (NKPT(KAO)-IO) .Nfc. ID) GO TO 360 
KONE=IUVULC(KAU) 

XNUMEKi J0)=( AUPREVt K(JNE)+AT DEN MESH, KONE) ) *VUL F I S ( KONE ) /2 . 

GO TO 21 

360 CONTINUE 

361 XNUMEK ( JW) = ( ADPKEVI I D ) +ATDEN ( ME SH t I D ) )*VULFIS< ID)/2. 

21 OENOM=PENUM+XNUM£k( JO) 

C 

C DENOM® T UT AL FISSION KATE PER BAKN-CM FUK ISUTUFES WITH 

C CUN S 1 DEK ED YIELDS 

C 

IF( DENUM.EU.O. ) GO TO 35 

G AMF P = 0 • 0 

GAMXE=0. 0 

GAMS M =0.0 

GAM I 0D=0 . 0 

G AmPM=0 , 0 

DO 26 J I G= 1 » NSPMAT 
26 GAMSH( JIG)=0.0 
DO 20 I 0= 1 1 NF I S 
KU=II)YLD( ID) 

GAMXE=GAMXE+XNUM£K ( I 0 ) ^Y ( 2 t KU ) /DENOM 
GAMSM=GAMSM+XNUMEK( IO)*Y< 3,KU) /DENOM 
GAM I OD = G AMI OD+XNUMEK ( I 0 ) * Y ( 4 , KO ) /DENOM 
GAMPM=GAMPM+XNUm 6R< IU ) *Y ( 5,KU ) /DENOM 
IF (NSPMAT. EU.U ) GO TU 20 
DO 2 7 J I G= 1 1 NSPMA T 
JOG® J I G+b 

2 7 GAMSP( JIG)=GAMSP( JIGJ+XNUMERl I u ) * Y ( JOG , KO ) / DtNUM 
20 GAMFP=GAMFP+XNUMEK( ID) *Y( 1 ,KD ) /DENOM 
GAIIMXe = GAMXE*UfcN0M 

gainiu=gamiod-denum 

GAINSM=GAMSM*UENOm 

GAINPM=GAMPM*DENQM 

c 

C FISSION PRODUCT AGGREGATE (OK INDIVIDUAL PRODUCT) 

C 

ADPKEV( 13)=ATU6N(MESHt 13) 

A T Dt N ( ME SHt 1 3 )= ( GAMFP* DENOM/ VULA6S ( 1 3 > ) * ( 1 . -EXP ( -DEL SEC*VUL ABS ( 13) 
1 ) )+AUPKEV(l3)*EXP(-DELSEC*VULARS< 13) ) 

I TbM=0 

IF (MOKmjT.LT. 2 ) GO TU 44 

WR I 1 E (6, 10/ ) 

WRITM6tll2) MESH 
WKI TP (6, 1 10) 

44 IF(NFKACT.EU.O) GO TO 34 
3/ ITEm=ITEM+1 

IF ( I (EM.GT.l ) GO TO 40 
DELSEC=36UO.*UbLHKvPKACT ( ITEM) 

GO ro 41 

40 DELSEC = 3600.*DELHR*(FRACT< I T EM ) -FKACT ( 1 T EM-1 ) ) 

41 IF ( ITEM.GT. 1 > GO TU 34 
SAVXE=ATDEN(MESHtl4) 

SAVI0=ATDEN(MESH,16) 

SAVSM=ATDEN( MESHt 15) 

SAVPM=ATDEN ( ME SHt 17) 

IF < MOROUT.lt .2 ) GO TO 34 
WR I TE (6» 102 ) 

C 

C EXPLICIT YIELD AND DECAY POISONS 

C 

C 

C PM 149 

C 

34 ADPKEVU7 )=ATDEN( MESH»17) 

GAMTOS=0.0 
GAMTC)X = 0.0 

ATDEN(MESH,17)=ADPKEV( 1 7 ) *EX P ( -DEL S EC* ( DEC A Y ( 1 7 ) + VUL ABS ( 17) ) ) 
i+(GAINPM/(DECAY( 17)+VULABS( 17) ) ) * ( l . -E X P ( -DEL S EC* ( DECAY ( 17)+ 
2VULABSU7) ) ) ) 

C 

C SMI 49 

C 

ADPKEV ( 15)=ATDEN(mESH, 15) 

I F ( ATD-EN(MESHtl 7) .NE.O. > GO TU 500 
GAMTOS=GAINPM+GAlNSM 

ATUENlMESH, 15) = ( AOPREVI 1 5 ) * EXP ( -DEL$EC*VUL AB S ( 15) )+GAMTUS* 

1 ( l.-EXP(-DELSEC*VULABS( 15) ) ) )/VULABS (15) 

GO TO 520 
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500 ATDEN(MESH, 15)= DECAYU7)* <ADPREV(17)*(0ECAY117)+VULABS(17))- 

lGAINPM)/( 1 DECAY ( 1 7 > + VUL ABS ( 1 7 ) ) * ( VUL ABS ( 1 5 ) -DECAY (17 ) -VUL ABS ( 17 > ) ) 

2 *( EXP(-DELSEC*( DECAY! 17) +VULABSI 17) ) )-EXP(-DELSEC*VULABS( 15) ) ) 

3 +(GAINSM*(UECAY< 1 7 ) + VUL ABS ( 1 7 ) ) +GA INPN*0ECAY < 17 ) )/ 

A! (DECAY! 17 ) +VUL ABS ( 17) ) 

5 *VULABS( 15) )*( 1 .- EXP ( -DELS EC*VUL ABS ( 15) ) ) + A0PREV(15J* 

6 EXp(-DELSEC*VULABS( 15) ) 

C 

C 1135 

C 

520 ADPREV! 16)=AT0EN(MESH, 16) 

ATDEN(MESH»16)=ADPREV{ 16)*EXP(-DELSEC*(D6CAY(40)+VULA8S( 16) ) ) 

1 + (GAIN 10/ (DECAY! 40 ) + VULABS( 16) ) )*( 1 .-EXP ( -DEL SEC* ( DECAY ( 40 ) 

2 +VUL ABS ( 16) ) ) ) 

C 

C XE135 

C 

30 ADPREVt 14)=ATDEN( MESH, 14) 

I F < ATDEN(MESH» 16I.NE.0. ) GO TO 510 
GAMTOX=GA I NI O+GA I NXE 

ATDEN (MESH, 14) =( ADPREV( 1 4 ) *EXP ( -DEL SEC * ( DEC AY ( 14 ) +VUL ABS ( 14)))+ 
IGArtTOX*! i.-EXP( -DELSEC*( DECAY ( 14)+VULABS( 14) ) ) ))/ (DECAY ( 14)+VULABS 
2(14) ) 

GO TO 530 

510 AT DEN (MESH, 14) =DEC A Y ( 1 6 ) * ( ADPRE V ( 1 6 ) * ( DEC A Y ( 1 6 > + VUL ABS < 16) )-GAINIO 
1) /( (DECAY( 16)+VULABS( 16) ) * ( D ECAY ( 1 4 ) +VUL ABS ( 14) -DECAY(16) 

2-VULABS( 16) ) )*( EXP (-DEL SEC* ( DECAY! 16)+VULARS( 16) ) ) -EXP ( -DEL SEC 
3* (DECAY! 1 4 ) +VUL ABS (14) ) ) ) + ( GA I NXE* < DECAY ( 16 ) +VULABS ( 16 ) J+GAINIO* 

4 DECAY (16) )/< (DECAY ( 16)+VULABS( 16) )* 

5 (DECAY! 14 ) +VULA8S 1 14) ) ) *( 1 .-EXP ( -DEL SEC* ( DECAY ( 14 ) +VULABS ( 14 

6) ) ) J+ADPREV! 14 )*EXP( -DEL SEC* (DECAY! 14)+VULABS( 14 ) ) ) 

530 IF(DELSEC. 60. SECOND) GO TO 38 
IF (NERACT.EO.O) GU TO 38 
IF (MOROUT.LT. 2) GO TO 300 

WRI TE(6, 101 ) FR ACT ( ITEM) , ATDEN( MESH, 1 4 ) , ATDEN ( ME SH , 1 5 ) t ATDEN (MESH, 
116 > , ATDEN! MESH, 17) 

300 IF ( I1EM.L1 .NFKACT) GO TO 37 
OELSEC=StCOND 
A10EN(MESH, 14) =SAVXE 
ATOENIMESH, 15)=SAVSM 
A T DEN ( ME SN * 16 ) = SAVIO 
ATDEN! MESH* 17)=SAVPM 
GO TO 34 

38 CALL XENON ( GA I NX E , GA I NSM , GA I N 1 0 , GA I NPM , AVGX E , AVG SM , A VG I 0 , A VGPM ) 

GO TO 39 

C 

C SPECIAL MATERIALS INCLUDING FISSION PRODUCT AGGREGATES 

C 

35 DO 42 NOD= 1*17 

42 ADPRE V ( NOD ) =0 . 0 
IFIKLON. EO.O) GO TO 382 
WRIT E (6, 107 ) 

WR ITE(6»112) MESH 
WRITE(6t 164) 

GO TO 39 

382 LUSH=MESH+NM I R ( KAT+1 )-l 
IF (MOROUT.LT. 2 ) GO TO 39 
WRITF(6* 107 ) 

I F ( NM IR(KAT+1 l.EO.l) GO TO 43 
WR ITE(6*15?) MESH, LUSH 
WR I TF ( 6, 1 53 ) MESH, LUSH 
GO TO 39 

43 WRITE(6,112) MESH 
WR I T E ( 6 , 1 00 ) MESH 

39 DO 362 KOT= 1 , N I SOT 
JOT=IDVlJLC(KOT) 

IF ( JOT.LE. 17 ) GO TO 362 

I F (NRPT (KOT ) .NE.O.AND.ATDENIMESH, JOT ) .NE .0. ) GU TO 362 
ADPRE V( JOT ) ~AT DEN( MESH, JOT ) 

362 CONTINUE 

IF ( DENOM. EO.O. ) GO TO 33 
IF(NSPMAT.EO.O) GO TO 33 
DO 31 I V = 1 , NS PMAT 
DACRUN=0.0 
DACRON=KRON( IV) 

NON= I V+ 1 7 

ADPRE V ( NON )= AT DEN (MESH, NON) 

ATDENtMESH, NON ) = ( (GAMSP( I V ) *DENOM+ ( ( ATDEN ( ME SH , NON-1 ) + 

1 ADPREV(N0N-1 ) )^VULABS( NON-1 )/2. )* DACRON) / VUL ABS ( NON ) ) * 

2 ( l.-EXP(-DELSEC v VULABS( NON ) ) )+ADPREV(NON)*EXP( —DEL SEC* VUL ABS ( NON ) ) 

31 CONTINUE 
C 

C KRON(I) IS THE KRONECKER DELTA — EOUAL TO 1 IF SPECIAL MATERIAL! I) 

C IS A DECAY PRODUCT OF THE SPECIAL MATERIAL! 1-1 ) AND EOUAL TO 

C ZERO OTHERWISE 

C 

33 DO 96 I S 60 = 1 , N I SOT 
PO I SON ( ISEO)=0.0 
96 FP I M ( ISEO)=D.O 

IF (NBPOI .EO.O) GO TO 48 
IF(KLON.EO.O) GO TO 48 
403 DO 32 LUP=1,N0POI 
L00P=LUP+NSPMAT+1 7 
ADPRE V ( LOOP )= AT DEN (MESH, LOOP) 

ATDEN ( MESH, LOOP )=ADPREV( LOOP ) *E XP ( -DEL SEC*VUL ABS ( LOOP ) ) 

32 CONTINUE 

48 IF! DENOM. EO.O- ) GO TO 90 
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C CONVERSION RATE AT THE START OF THE TIME INTERVAL 

C 

CONSB6*0,0 
OE STBE®()*0 

CONS BE® AOHREVm^VULAMFl I ) +AOPREV ( 7 ) *VULAMF ( 7 H-AOPRE V (4 > *VULAMf < 4 
l ) 4-AOPREVf LimVULAHH 10 ) 

OFSTBE* A0PREV(3)*VULABS(3) +AOPR£V( 5 ) *VULAt)S < 5 ) +ADPKE V 1 9 ) *VULABS t 
14)4- ADPREV t II) OVULARS < II ) + ADPREV ( 2 ) *VUL ABS ( 2 > *OECAY < 2 ) f t DECAY I 2 
2 ) +VUL ABS ( 2 ) )+ADPREV(8)*VULABS<8)*DECAY(8>/<0gCAY(8) +VULA8S(8) ) 

IF (NFIRPT-EO.OI GO TO 340 
00 541 NWT® 1 » N I SOT 
I F ( NRP T ( NWT ) .EO.O ) GO VO 341 

ig r*i ovuic (Nwn 

KROG=NRPT(NWTi -10 

GO TO ( 342 » 344 t 343 » 342 » 343 » 343 , 342 , 344, 343 , 342 , 343 , 343 ), KROG 

342 CON$ftE=C0NSBE + ADPKEV< 1 GT ) *VUL AMF < IGT ) 

GO TO 341 

343 DhSTBE=DESTBE+ADPRtV< IGT ) -VULAB5 t IGT ) 

GO TO 341 

344 DES TBE = DE STB E + ADPREV l IGT ) +VULABS t IGD^DEC AY l KROG ) / ( DECAY t KROG) FVUL 
lAttSUGr) 3 

341 CONTINUE, 

340 TCUNSB-rCONSBFCQNSBE* VOL (MESH) 
roesrB® rOFSTf5 + OeSTBE*VOL<NESH) 

TK8EF* fCONSB/TOESTB 

81 IF (DESTBE.EW.Q. ) GO TO 93 
RA TBEF { MESH ) “CONSBE/DFSTBE 

93 CONTINUE 
C 

C FUEL INVENTORY CALCULATION 

C 

(X) 15 J INX= 1 ,12 

BEGMA$(JINX)=6EGMAS( JINX) +VnU MESH) *ADPREV( JINX) *ATWGT( JINX)/ 
l { AVUGAO-* 10 ■ ## { -21)) 

15 F INMAS ( J I NX ) ®F l NMAS ( J r NX ) + VOL ( MESH ) #ATDEN ( MESH, J I NX ) *ATWGT C J I NX ) 

L iAVDGAD*10.**{-21) ) 

FUEL =0.0 
00 fl IU=l»l2 

r 1 FUEL® FUEL FATOEN l MESH » I D) ^ VULABS ( 10 ) *VOL ( MESH) 

PETROL =PETROL+FU£L 
IF(NFIRPT.EO.O) GO TO 51 
HO 49 Ml T® 1 , NI SOT 
IFtNRPT(MIT) .CE.O) GO FU 49 
META*IDVULC<MIT) 

MUT r = NRPT(MIT)-lO 

PUEL“FUf=l+ATDENt MESH, META ) *VUL ABS ( MET A ) *VUL (ME SH ) 

BEGMAS f MUTT ) =8EGMA5 f MUT T ) + VOL ( MESH) ^ADPREV I META ) *AT WGT (MUTT >/ 
l (A'VOGAOMO.**! -21 ) ) 

F INMAS I MUTT )=F1NHAS( MUTT) + VOU MESH) #AT DEN( M£SH,ME I A ) *ATWGT( MUTT ) / 

L lAVnGAl)*10.**(-2L)) 

44 CONTINUE 

PETROL “PETROL + FUEL 
r 

C POISON CALCULATION 

c 

51 PHI roT«0.O 

00 72 10*1 tNI SOT 
( I “ IOVULC ( ID) 

P 1 2 £ N ( Il))“Pl/ENt ID) + A TDENl MESH » I T ) OVULARS ( IT ) *VOL IMESH) 

12 POISON! IO)=ATOEN(MESH, IT)*VULAB5( IT ) *V0L ( MESH ) /FUEL 

C 

C POWER CALCULATION 

75 CONTINUE 

DO 92 ll)*l»NISOT 

1 00*1 OVULC ( 10) 

rF< IOO.Gt. 12. AND. NrtPT( 10) .LE. 01 GO TO 92 

78 SUM = 0 *0 

DO 91 IG=1,NG 
LOC“NM*( IG-D+MESH 
9 1 SUM^SUMF SANSNU I IDO, I G ) *FLUX ( t DC ) 

MINI IO)*ATOEN(MESH, 100 ) * 1 0 ( 24) *SUM*VUL( MESH ) /F l WATT 
fOPOW ( MESH ) “TUPOjfH MESH)+FPIM( ID) 

92 CONTINUE 

n IF< rOPDW(MESH) .60,0. ) GO TO 90 
IGLOO® IGLOO* l 
f)(J 89 IOTA®I,NlSOr 

TFRAP( IOTA) =TFRAP (IOTA >+FPIM( IOTA) 

89 FP I M ( IOTA)=FPIM( IOTA ) /TOPOW ( MESH) 

90 CONTINUE 

NBEC=17+NSPMAT+NBP0I+NFiRPT+l 

82 CUN T I NUE 

00 94 ISFU=1,NIS0T 
F AO I N ( ISEU)=0,0 

IDT=idvulc< iseo) 

I F ( ADPREV (IDT) »£U .0 . ) GO TO 94 

FADINI I SCO ) “ATDENf MESH, IDT )/AOPREV( IDT) 

PREDEN ( I SEU ) * AOPREV ( IOT ) *VOU MESH ) +PRE06N ( ISEO) 

94 CONTINUE 

iPtM0RDUT.LT.2J GO TO 47 

45 WRITE (6, 129) 

79 CONTINUE 

00 95 JS£Q=ltNISOT 
10*t OVULC < JSEQ) 

95 WRITE (6t 120) T NAME ( I D >, J SEO » ADPREV ( 10 ) , ATDEN ( ME SH, ID ) , FADI N ( JSEO ) , 
IPPLMt J5E0) , POISON! JSEQ) 

47 POUR = POUR + TOPUW t MESH ) 

IF(DEN0M.GT.Q*.0R.KL0N.E0.1 ) GO TO 320 

M0$T=MESH+1 

DU 321 L1TE“M0ST ,LUSH 
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DO 322 N I TE=1, N I SOT 
MI TE=I DVULC ( N I TE ) 

322 PREDEN (NITE)=ADPREV( MITE) *VOL( LITE) +PREDEN { NITE) 

321 CONTINUE 
MESH=LUSH 
320 MESH=MESH+1 

IF(MESH.LE.NM) GO TO 1 
POWIN=PQWER 
NQNE=KTIM— 1 

IFIN0NE.GT.1) WRITEC6,L07) 

WRITE(6,156) NONE 
CALL ABSPHI 
C 

C PRODUCTION RATE OF FISSIONABLE MATERIAL 

C USING NEW FLUX LEVEL AT END OF INTERVAL 

C 

SP ACE=0* 0 

DO 87 MESH=L,NM 

TQPQVl I MESH ) -TOPOW t MESH 1 -FNORM 

I F ( TOPOW t MESH ) .NE.O.) SPACE=SPAC E+VOL ( MESH } 

DO 85 IOTA=L, N ISOT 
I D=I DVULC t IOTA) 

VULABS! ID)=0.0 
VULFESt ro)=o.o 
VULAMFI I D ) =0 • O 
DO 86 I G=L tNG 
LOC=NM*( IG-D+MESH 

VULABS C ID ) =VULABS ( ID)+ABSIG( ID, IG INFLUX (LOC) 

86 VULFISt ID)=VULFISf ID J+FLUX t LOC) ^SANSNU l ID, IG) 

85 VULAMFI I D ) =VULABS I I D ) -VULF IS ( ID) 

PR0D=0.0 

PROD * ATDENI MESH, L) ^VULAMFI I )+ATOEN tMESH,7 ) ^VULAMF 1 7 ) +- 

1 ATDEN t MESH , 4) ^VULAMF t 4) +ATDEN I MESH , 10 ) 4VULAMF 110) 

C 

C DESTRUCTION RATE OF FISSIONABLE MATERIAL 

c 

OEST-O. 0 

DEST = ATDEN t MESH , 3 ) * VULABS 1 3 ) -t-ATDEN l MESH, 5 ) ^VULABS ( 5 ) 

1 + ATDEN! MESH,9)*VULABSt 9 J+ATDENI MESH, il )# VULABS ( 11 ) +ATOEN I MESH, 2 ) 

2 *VULAB5(2)-0ECAY(2)/(DECAYI2)+VULABS (2) ) +ATDEN I MESH, 8 ) 4 
3VULABS ( 8 ) ^DECAY ( 8 ) / I DECAY ( 8 ) +VULABS I 8 ) ) 

IF(NFIRPT.EO.O) GO TO 350 
DO 351 NWT=L, N ISOT 
IF(NRPTtNWT) .EO.O) GO TO 351 
IGT=I DVULC I NWT ) 

KROG^NRPT ( NWT) —10 

GO TO (352, 354,353,352, 353, 353, 352, 354, 353, 352, 353, 353), KROG 

352 PROD=PROD+ATDEN(MESH, IGT ) *VULAMF t IGT ) 

GO TO 351 

353 DEST=DEST+ATDEN(MESH.,IGT)-VULABS( IGT) 

GO TO 351 

354 DEST=DEST+ ATDEN ( MESH, IGT ) ^VULABS { IGT ) *DECAY ( KROG ) / I DECAY ( KRGG ) + 

1 VULABS I IGT ) ) 

351 CONTINUE 

350 TPROD=TPROD+PROD-VOL(MESH) 

TDEST=TDEST+DEST*VOLlMESH) 

IFtDEST.EO.O. ) GO TO S7' 

C 

C CONVERSION RATE OF FISSIONABLE MATERIAL 

C 

CONRAT ( MESH ) =PROD/DEST 
GO TO 98 

97 CONRAT! MESH) =0.0 

98 TCONV=0.0 

IF ! TDEST.EO.O. ) GO TO 87 
TCONV=TPROD/TDEST 

87 CONTINUE 
AVGXE=AVGXE/SPACE 
AVGSM=AVGSM/SPACE 
AVGIO=AVGIO/SPACE 
AVGPM=AVGPM/ SPACE 
FIGLOO=IGLOO 

IF ! NFIRPT . LE. 0 ) GO TO 370 
DO 371 JAR=1,NIS0T 
IF! NRPTt JAR) . LE-O) GO TO 371 
MATTE=NRPT( JAR) -10 

no 410 loch=i,nisot 

IF! MATTE. EO. I DVULC l LOCH) ) GO TO 411 

410 CONTINUE 

411 TFRAP! LOCH) =TFRAP( LOCH) +TFRAP! JAR) 

371 CONTINUE 

370 00 83 JAM=1,NIS0T 
83 TFRAP ( JAM ) =TFRAP ( JAM) /POUR 
DO 73 ID=L,NISOT 

73 TAX (ID) —0. 0 

DO 74 MESH=1,NM 
DO 74 KEN=1,NIS0T 
I D= I DVULC (KEN) 

IF! ATDEN (MESH, ID) ..NE.Q. ) VOLL! K EN) =VOLL IKEN ) +VOL (MESH) 

74 TAX ( KEN ) =TAX ! KEN) +ATDENC MESH, ID ) *VOL ( MESH) 

KTL=KTIM-L 

WRITE! 6, 105) KTL 
WRITE(6, 121 ) 

DO 80 MO I D=1 ,N I SOT 
I D= I DVULC ( NO I D ) 

WRITE! 6 , 104) ID 

WRITE 1 6 , 106 ) l ATDENI IMESH, ID) , IMESH=1,NM) ,TAXtNO!D) 

80 CONTINUE 

WRITE !6, 113 ) KTL, (RATBEF! ION) , I0N=1,NM) 

WR I TE (6,1 14) KTL, (CONRAT! ION ) , I0N=1 , NM ) 

WRITE (6, 109 ) 
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WRITE (6, 117) 

00 550 JME$H= 1 t NM 

550 TUPOW! JMESH)=T0P0W< JMESH)/V0l ! JMESH) 

WRITE (6» 106) (TOPOW! JMESH) ,JM£SH=1 ,NM), POWER 

on hh jerk=i»nm 

88 TQPOW! JERK ) =TOPOW( JERK)/ AVPOW 
WRITE(6»117> 

WRITEI6, 115) 

WR J TE ( 6 » 1 1 7 ) 

WRITE (6, 106) < TOPOW! JMESH) ,JMESH=1 »NM) ♦ AVPOW 

WR I TE ( 6 ♦ 1 IQ ) 

WRITE<6,119) POWER 

WRITE I 6, HR) AVGXfc,AVGSM f AVGIOrAVGPM 
WRITE (6, 111 ) 

DO 330 I D= 1 1 N I SOT 
NQW = I OVULC I ID) 

IFIN0W.GT.12) GO TO 330 

WRI TE (6, 15R ) T NAME < NOW) »BEGMAS < NOW) . F I NMA S ( NOW ) . TFRAP { I D ) 

330 CONTINUE 

DO 301 I D= 1 » N I SOT 

PI ZEN! ID) = PI 26N( ID) /PETROL 

PREDEN ( ID)~PRED£N( ID) /VOLK ID) 

301 TAX! ID)=TAX( ID) /VOLLt 10) 

WKlTF(6t 107 ) 

WRITE(6t 157) 

WRITE(6t155> 

311 CONTINUE 

DO 310 I D= 1 » N I SOT 
J S F0= I OVULC I ID) 

312 IF I PREDENI ID) . EO.O. ) GO TO 313 
XSU = TAX I ID) /PREUENI ID) 

GO TO 310 

313 XSU=0.0 

310 WR I TEI6, 120) TNAMEI JSEO) » J SEO , PKEOEN ( I D ) , T AX I I 0 ) , X SU , P I 2 EN I I D ) 

DO 84 L AWN= 1 » 1 2 
BMASS=BMASS+BEGMAS(LAWN) 

64 FMASS=FMASS+FINMAS! LAWN) 

BEGLB=2. 2046* B MASS 
FINLB=2.2046*FMASS 

DEGBUR = 2 200. SHOWER* DELHR/ I 24 . * ( B EGL B-F I NL B ) ) * 1 0 . ** ( -6 ) 

WKITEI6, 1 16) 

WR I TE I 6 » 1 40 ) BMASStFMASS 
WR II E I 6 ♦ 1 4 1 ) BEGLBfFlNLB 
WR I TE I 6 1 1 08 ) TRBEF » TCONV 
WR I T E ! 6 » 1 03 ) ROW I N » POUR 
WRITE (6tl54) DEGBUR 

100 FORMAT ( IHL ♦ 80H NO FISSIONABLE ISOTOPES OR CONSEQUENT FISSION PROD 
1UCTS OCCUR IN MESH INTERVAL rl4) 

101 F0KMAT(F9.2,5X,4E15.6) 

102 FORMAT l 7 1 H FISSION PRODUCT CONCENTRATIONS AT VARIOUS FRACTIONS 0 
IF THE TIME STEP//12H FRACT I ON , 9X » 5HXE 1 35 t 1 OX , 5HSH 149 , 10X t 5H I 13 
25»10X»5HPM149 ) 

103 FORMAT! lHK t 24H TOTAL POWER IN HA T T S , 1 OX t E 1 2 . 6 , 10X t E 1 2 .6 ) 

104 FORMAT! lHLt 15H VULCAN ID = tI2) 

105 FORMAT ( 1H1 ,84H NEW ATOM DENSITIES BY MESH INTERVAL FOR EACH I SOT 
10PE AT THE END OF TIME INTERVAL *12) 

106 FORMAT ( 8 E 1 4 . 6 ) 

107 FORMAT I 1H1) 

108 FORMAT! 1HK,24H CONVERSION RATIO , 1 OX , E l 2 . 6 * 10X , E 1 2 .6 > 

109 FORMAT ( 1HK»34H BY MESH INTERVAL THE POWER IS ) 

110 FORMA T ( 1 HK ) 

111 FORMAT! lHl t 64H FISSIONABLE INITIAL MASS IN FINAL HASS IN 

1 FRACTI0NAL/6X,7HlS0T0PE,9X,yHKlL0GRAMSt9X,9HKlL0GRAMS,9Xt5HPUWER 
2/) 

112 FORMAT ! 1 0 < 3H **),22H MESH INTERVAL NUMBER ,I4,1X,10!3H **)) 

113 FORMAT! 1H1,63H CONVERSION RATIO BY MESH INTERVAL AT START OF TI 
1ME INTERVAL»I2//(8E14.6) ) 

114 FORMAT! 1HK»63H CONVERSION RATIO BY MESH INTERVAL AT END OF TI 
1 ME INTERVALt I2//!8E14.6) ) 

115 FORMAT ( 73H BY MESH INTERVAL THE POWER (NORMALIZED TO VOLUME-AV 

1ERAGED POWER) IS ) 

116 FORMAT! 1 HL » 37X , 7H I N I T I AL ♦ 1 6X , SHF I NAL ) 

117 FORMAT (1HJ) 

1 IB FORMAT! lHLt 14Xt 14HXEN0N-135 , E20 . 6/ / 1 5X , 14HSAM AK I UM-149 »E20. 

16//15X* 14HI0DI NE-135 ♦ £20 . 6/ / 1 5X » 1 4HPR0ME TH I DM- 149 « E 20 . 6 ) 

119 FORMAT! 75H VOLUME AVERAGED EQUILIBRIUM CONCENTRATIONS OF MAJOR 

1 FISSION PRODUCTS /45H IN NUCL E I /BARN*CM IF THE POWER LEVEL 0 

2F f E 1 1 • 6 * 20H WATTS IS MAINTAINED) 

120 FORMAT! 3X» A5* 1 H-, I 2 * IX 1 6E 1 5 . 6 ) 

121 FORMAT! 1HK , 20 < 4H**** ) / / 103H LAST ENTRY FOR EACH ISOTOPE 

IIS ONE BAKN-TH OF ITS TOTAL NUMBER OF ATOMS IN THE CONFIGURATION// 
21X,20(4HW**4) ) 

129 FORMAT! 1HK»9H I SOTOPE » 10X » 4HN ( T ) *9X , 6HN ( T+ 1 ) , 7X , 1 1 HN ! T+l ) /N ( T > ♦ 7X 
1 T 5HPOWER f 9X,6HP0ISON/10H DESCK . t 53X , 8HFRACT I ON t 7X , 6HFACT0R ) 

131 FORMAT! 1HK,26H IN MESH INTERVAL NUMBER , 14// 15X , 58HTHE CONVERSION 
1 RATIO AT THE START OF THE TIME INTERVAL IS « F 1 0 . 6// 1 5X t 56HTHE CON 
2VEKSI0N RATIO AT THE END OF THE TIME INTERVAL IS ,F10.6) 

140 FORMAT! lHLt24H TOTAL FUEL I NVENTOR Y , 1 OX , E 12 . 6 , 10X , E 1 2 . 6/9X , 
112HIN KILOGRAMS) 
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141 FORMAT! 1HK,24H TOTAL FUEL I NVE NTOR Y , l OX , E 1 2 . 6 , 1 OX , E 1 2 . 6/ 10X , 
19HIN POUNDS) 

150 FORMAT ( 7E15.M 

151 FORMAT! lHLt 10( 1H+) ,95H ALL ZONES HAVE BEEN SEARCHED FOR A NON-DEP 
1 LETABLE ZONE STARTING POINT AND NUNE HAS BEEN FOUND/ IX , LO < 1H+ ) , 

?9H MS UN = » l 3 , 1 2H AND MESH = ,13) 

152 FORMAT (10 1 3H *'~),23H MESH INTERVALS NUMBER ,14, 6H THRU ,14, 

1 IX , 10 I 3H ) 

15 3 FORMAT ( 1 H L , H 1H NO FISSIONABLE ISUTOPES OR CONSEUUENT FISSIUN PROD 
1UCTS OCCUR IN ME S'H INTERVALS ,14, 6H THRU ,14) 

154 FORMAT I 1HK , 2HH MEGAWATT-DAYS/METRIC TON , 1 /X , E 1 2 .6 ) 

155 FORMAT ( 1HK, 9H I SOTOPE , 1 OX , 4HN ( T ) ,9X, 6HN< T + l ) , 7X » UHN { T + l ) /IM (T) , 7X 
1 ,6HPOISON/10H DESCR. , 54X , 6HFAC TOR ) 

156 FORMAT ( 1HL, 60H ALL MESH INTERVALS HAVE BEEN TRAVERSED FOR TIME I 
INTERVAL ,12 / 64H THE FLUX WILL NOW BE RENORMALIZED FOR THE NEXT 
2 T I ME INTERVAL) 

157 FORMAT ( 1HK, B4H THE FOLLOWING EDIT IS AVERAGED OVER ALL MESH INTEKV 
1ALS IN WHICH EACH ISOTOPE OCCURS) 

158 FORMAT ( 1HK,6X,A5,E20.6,E18.6,E16.6> 

164 FORMAT ( 1HL » 56H CUNTAINS BURNABLE POISONS BUT NO FISSIONABLE ISOT 
10PES) 

RETURN 

END 
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JIBFTC FADE DECK 

SUBROUTINE FISEQK*) 

C 

C 

C SOLVES DEPLETION EQUATIONS FOR FISSIONABLE ISOTOPES 

C 

C 

COMMON/SET 1/ ATDEN( I90» 40 ) ,CONRAT ( 200 > , TOPOWI 200 ) , VOL ( 200 ) ,RR ( 200 ) 
1 tRATBEF ( 200 > 

C0MM0N/SET2/ FLUX ( 3800 1 , SANSNUl 40 ,20 ) , ABS I G! 40, 20 ) , SST < 20, 20 > , 

1 XNU ( 1 2 *20 ) 

C0MM0N/SET3/ NMI R( 200) ,RAUM( 200 ) *RAR ( 200 ) , DELR ( 200 ) t MI DRZ ( 400) » 
1NEM0 ( 200 ) 

COMMON/ SET4/ A TWGT ( 20 ) , DECA Y< 40 ) ,THERNU ( 40 ) , TNAM6 ( 40 ) , FRACT ( 10 ) , 
1NRPT ( 40 ) * I DVULC t 40 ) 

COMMON/SET5/ VULABS ( 40 ) ,VULFIS( 40 ), VULAMF ( 40) , ADPRE V ( 40 ) , I DYLD( 10) 
1 »Y(20,12) , KRON ( 20 ) 

C0MM0N/SET6/ NG, NM, N I SOT » TOTVOL , F I WATT, POWER , AVPOW, FNORM 

C0MM0N/SET7/NFIS, NSPMAT ,NBPO I 1 N FRACT ,NYDNUC » NF IRPT 

COMMON/SETS/ DELHR , KHA I N, A VOGAD ,NZONR, MOROUT ,MOOD, DI FF ER * NONDPL 

COMMON/ SET9/ KT I M , NT I NC ,KRELL t KCELL 

COMMON/SET 10/ KOMP, KONG, KOG, LAND, DEL SEC , MESH 

Jl = l 

J 2 = 2 

J3=3 

J4=4 

J5 = 5 

J6 = 6 

GO TO 201 
ENTRY FISEQ2(*> 

22 J 1-7 
J2 = 8 
J3=9 
J4= 1 0 
J5= 1 1 
J6 = 1 2 

8 GO TO (201,202, 203 , 204*205,206) , KOMP 
C 

C ENTER THE ATOM DENSITY CALCULATIONS FOR FISSIONABLE ISOTOPES 

C 

C 

C TH232 OR U238 

C 

201 ADPREV(J1 ) =ATDEN ( MESH, J 1 ) 

ATDEN I MESH, J 1 >=ADPREV(J1 >*( 2. -VULABS (Jl )*DELSEC ) / ( 2.+VULABS ( J 1 )* 

1 DELSEC) 

IF ( LAND.EQ.l ) GO TO 7 
C 

C PA233 OR NP239 

C 

202 ADPREV( J2)=ATDEN(MESH,J2) 

ATDEN (MESH, J2)=(ATDEN( MESH, Jl )+ADPREV(Jl ) ) * VULAMF ( Jl )/ (2.*( 

IDECAYi J2 ) +VULABS ( J2 ) ) ) + ( ADPRE V ( J2 )-( ATDEN ( ME SH , J1 ) +ADPR EV ( J 1 ) > 
2*VULAMF< Jl)/(2.*( DECAY! J2 > +VUL ABS ( J2 ) ))) *EXP (-( DECAY ( J2 ) + 

3VULABS ( J2 ) ) *DELSEC ) 

IF(ATOEN(MESH, J2J.EQ.0* ) ATDEN ( MESH, J2 ) =-ATDEN ( MESH, J2 ) 
IFUAND.EQ.l > GO TO 7 
C 

C U233 OR PU239 

C 

203 ADPRE V ( J3 ) -ATDEN ( MESH, J3 ) 

DO 30 IKE-1 , NI SOT 
NIKE=I DVULC ( IKE) 

I F ( NIKE. NE* J2 ) GO TO 30 
GO TO 32 

30 CONTINUE 

ATDEN (MESH, J3)=(ADPREV( J3)*l 2 .-VULABS ( J3 ) *DELSEC ) + ( ATDEN ( MESH, Jl) 

1 +ADPREV ( J 1 ) )* VULAMF! Jl ) *DELSEC > / < 2 . +VUL ABS { J3)*DELSEC) 

GO TO 31 

32 ATDEN (MESH, J3)=( ADPREV( J3 )*( 2 .-VULABS < J3 )*DELSEC ) + ( ATDEN ( MESH, J2)+ 
1ADPREVIJ2) )* DECAY IJ2)^DELSEC)/ (2 *+VUL ABS ( J3)*DELSEC) 

31 I F ( LAND* EQ* 1 ) GO TO 7 
C 

C U234 OR PU240 

C 

204 ADPREV(J4)=ATDEN(MESH,J4) 

ATDEN ( MESH, J4)=(ADPREV( J4)*( 2. -VULABS (J4)*DELSEC)+( ATDEN( MESH, J3)+ 
1ADPREV( J3) )*VULAMF( J3)*DELSEC+( ATDEN ( MESH, J2 ) +ADPREV ( J2) )* 

2VULAMF< J2)*DELSEC)/( 2*+ VULABS (J4)*DEL SEC) 

IF(LAND.EQ* 1 ) GO TO 7 
C 

C U235 OR PU241 

C 

205 ADPREV( J5 ) =ATDEN ( MESH, J 5 ) 

IFU5.EQ.il) GO TO 29 

I F ( KONG* EQ* 1 1 ) GO TO 29 
C U235 

ATDEN (MESH, J5)=( ADPRE V( J 5 ) * ( 2 VULABS ( J 5 ) *DELSEC ) + ( ATDEN ( MESH, J4) 
1+ADPRE V ( J4 ) >*VULAMF( J4 ) *DEL SEC ) / ( 2 .+VUL AB S ( J5)*DELSEC) 

IF ( LAND* EQ* 1 ) GO TO 7 
GO TO 206 
C PU241 

29 ATDEN! MESH, J5)=(ADPREV(J5)*< 2. -(VULABS (J5)+DECAY( 11 ) )*DELSEC)+ 

1 (ATDEN! MESH, J4)+ADPREV(J4> )* VULAMF ( J4 ) *D6LSEC >/< 2* +( VULABS ( J 5 > + 

2 DEC AY (11) )# DELSEC) 

I F ( LAND. EQ- 1 ) GO TO 7 
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U236 OR PU242 


C 

206 ADPREV (J6)=ATDEN(MESH»J6) 

AT DEN (MESH,J6>=( ADPREV ( J6 ) * ( 2 •— VULABS (J6)*DELSEC)+(ATDEN(MESH»J5) 
1+ADPREV( J5) )*VULAMF< J5 ) *DELSEC > / < 2 ■ +VUL ABS < J6)*DELSEC) 

I F ( LAND* EQ» 1 ) GO TO 7 
I F ( J6.EQ.12) GO TO 25 
I F ( KHAI N. EO . 0 ) GO TO 22 
C 

C FOR DUPLICATE FISSIONABLE ISOTOPES 

C 

25 IF(NFIRPT.EQ.O) RETURN I 
L AND= 1 

7 I F ( KOG.GT *NI SOT ) RETURN I 
DO 5 K I M^KOGt N I SOT 
IF(NRPTIKIM) .EO-O) GO TO 5 
JADE = I DVULC (KIM) 

IF(ATDEN(MESH,JADE) .EG.O.) GO TO 5 

KOG=K I M+l 

KOMP=NRPT ( K I M ) — 10 

I F ( KOMP. GT • 6 ) GO TO 320 

Jl = l 

J2 = 2 

J 3 = 3 

J4 = 4 

J5 = 5 

J6 = 6 

GO TO 330 
320 J l =7 
J2=8 
J3 = 9 
J4=10 
J5 = 1 1 
J6 = I2 

330 GO TO <301*302, 303, 304, 305,306, 301, 302, 303, 304, 305, 306), KOMP 

301 J1=IDVULC(KIM) 

GO TO 315 

302 J2= I DVULC (KIM) 

IF <NRPT( KIM-1 ) .EQ.ll) J 1= I DVULC (KIM-1 ) 

I F < NRPT < K I M— 1 ) . 60 * 1 7 ) J1 = I DVULC ( K I M-l ) 

GO TO 315 

303 J3= I DVULC < K I M ) 

IF ( NRPT ( KIM— 1 ) « EO* 12 ) J2= I DVULC < K I M-l ) 

IF(NRPT(KI M— 1).EQ»18) J2= I DVULC ( K I M-l ) 

GO TO 315 

304 J4= IDVULC(KIM) 

IF(NRPT(KI M-l ) • EQ* 13 ) J3= I DVULC < K I M-l > 

IF (NRPT <K I M-l > -EQ.19) J3= I DVULC ( K IM-1 ) 

IF (NRPT( KIM-2) .E0.12) J2= I DVULC (K IM— 2 ) 

IF(NRPT(KI M-2 ) *EQ * 1 8 ) J2= I DVULC ( K IM-2 ) 

GO TO 315 

305 J5= I DVULC (KIM) 

IF (NRPT (KIM-1 ) .EQ.14) J4= I DVULC ( K I M- 1 ) 

IF (NRPT (KIM-1) .EQ.20) J4= I DVULC ( K IM-1 > 

GO TO 315 

306 J6= I DVULC (KIM) 

IF(NRPT(KIM-1).E0.15) J 5= I DVULC ( K I M-l ) 

IF (NRPT (KIM-1 ) .EQ.21 ) J 5= I DVULC ( K IM-1 ) 

315 CONTINUE 
KONG=KOMP 

IF(K0MP.LE-6) GO TO 8 
KOMP=KOMP— 6 
GO TO 8 
5 CONTINUE 
RETURN 
END 
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SIBFTC POUT DECK 

SUBROUTINE OUTPUT 

COMMON/ SET1/ ATDEN( 190,40) , CONR AT < 200 ) , TOPO W < 200 ) , VOL ! 200 ) ,RR(200) 
1 , KATBEF ( 200 ) 

COMMON/ SET 2/ FLUX! 3B00) , SANSNU ! 40 , 20 ) ♦ ABS IG(40 ,20 ) , SST ( 20,20 ) , 

1 XNU ( 12,20) 

C0MM0N/SET3/ NMIRI200) ,RAUM( 200 ) , RAR ( 200 ) , DELK ( 200 ) ,MIDKZI400), 
1NEM0 ( 200 ) 

COMMON/ SET4/ATWGT! 20) , DECAY ( 40) , THERNU ( 40 ) , TNAME { 40 ) , FKACT ( 1 0 ) , 

1 NRP T ( 40 ) , I DVULC ( 40 ) 

COMMON/ SET5/ VUL ABS ( 40 ) , VULF I S ! 40 ) , VUl AMF ( 40 ) , ADPKEV ( 40 ) , I DYLD! 10) 
1 , Y ( 20, 12 ) r KRON( 20 ) 

COMMON/ SET6/ NG,NM,NISOT,TOTVOU»FIWATT,POWER, AVPOW,FNORM 
common/ setb/ OELHK,KHAIN,AVOGAD,NZONR,MOROUT ,MOOD,DIFFEK,NONOPL 
DIMENSION LABEL! 40) »CARRY(40) 

C 

C OUTPUT EDIT - DETERMINATION OF THE NUMBER OF MATERIALS TO BE 

C NEEDED FOR THE NEXT TDSN CALC - CALCULATES AVERAGE ATOM DENSITIES 

C AFTER CHOOSING ON THE BASIS OF MACROSCOPIC COMPARISONS 

C 

100 FORMAT ( 1H1 ) 

102 FORMAT ( 1HL»46H THE ATOM DENSITIES TO BE USED FOR THE NEXT ,12, 
117H MESH INTERVALS ( ,I3,5H THRU ,13, 5H ) ARE//) 

103 FORMAT (8E14. 6) 

104 FORMAT ( 2H* ,7110) 

105 FORMAT ( 7E 10 • 5 ) 

106 FORMAT! 1HL, 10X,23HMACGG INPUT CARD IMAGES//) 

107 FORMAT! 5X,7E14. 6) 

108 FORMAT.! 1H1,88H MACROSCOPIC ABSORPTION CROSS SECTIONS BY MESH INT 
1 ERVAL , BY GROUP, AND BY VULCAN ISOTOPE) 

109 FORMAT! LHL, 24H MESH INTERVAL NUMBER ,13) 

110 FORMAT! 1HK,23H tNERGY GROUP NUMBER ,13//) 

111 FORMAT! 1HL,25H MESH INTERVALS NUMBER ,13,6H THRU ,13) 

112 FORMAT! 56H AVERAGED ATOM DENSITIES FOR NEXT SPATIAL CALCULATION 
1 ) 

151 FORMAT! 1HL, 10! 1H+) ,95H ALL ZONES HAVE BEEN SEARCHED FOR A NON-DEP 
1LETABLE ZONE STARTING POINT AND NONE HAS BEEN FOUND/ IX , 10 ( 1H+ ) , 

29H MSUM = , 1 3 , 12H AND MESH = ,13) 

LG0 = 0 
L S T 0 P = 0 

IF!MC)R0U1 .EO.0.0R.MUK0UT.E0.2 ) GO 10 57 
C 

C THROUGH STATEMENT 42 IS THE MACROSCOPIC PRINTOUT 

C 

WRITE (6, 108) 

KOKO = N I SOT -NON DP L 
L I M = 0 

55 L I M = L I M+ 1 
MSUM=0 

DO 50 MI D= 1 , KOKO 

IF! ATDENl LIM»MID) .NE.O. ) GO TO 52 

50 CONTINUE 

DO 51 K AT = 1 , NZONR 
MSUM=MSUM+NM I R ( KAT ) 

IFIMSUM.EW. ! LIM-1 ) ) GO TO 53 

51 CONTINUE 

VIR ITE16, 151 > MSUM,LIM 

53 IF(NMIR(KAT+1 ) .EO. 1 ) GO TO 52 
LUSH=L I M+NMI K l KAT+i ) -1 

WRITE (6, 111) LIM,LUSH 
GO TO 54 

52 WKITEI6,! 09 ) LIM 

54 DO 43 L I G = 1 , MG 

DO 44 M I B= 1 y N I SOT 
LID=IDVULC(NIB) 

IF ( ABSIG!LID,LIG) .EO.O. ) ADPREV ( L I D ) =ATDEN ( L I M, L I D ) 

44 ATDEN!LIM,LII))=ATDEN(L I M , L I D ) #AB S I G ! L ID,LIG) 

WRITE (6, 110) LIG 

WKI TE(6, 103 J ! ATDtN! L IM , L I D > ,L ID = 1 , N I SOT ) 

DO 45 N I B= 1 , N I SOT 
L I D= I DVULC ! N I B ) 

IF(ABSIG (LID, LIG). EO.O. ) GO TO 47 
ATUEN(LIM,LID)=ATUEN(LIM,LID)/ABSIG(L IDtLIG) 

GO TO 45 

47 ATDEN(LIM,LID)=ADPREV( LID) 

45 CONTINUE 
43 CONTINUE 

IF! MSUM. NE - 0 ) LIM=LUSH 
42 IF! LIM.LT.NM) GO TO 55 
C 
C 

57 ME SH = 0 

IFtDIFFEK.EW.O. ) DIFFER=.001 
C THE FOLLOWING LOOP IS FOR THE ZONE SWEEP 

5 00 10 IX=1, NZONR 
LG0 = LST()P + 1 
LSTOP = LGO + NMIKt I X ) — 1 
MAR S=0 
ME TK0=0 
ME SH = M E SH + 1 
J A V E = 0 

ME TRO = NM l K ! I X ) - 1 
I F ! ME TRO . EO. 0 ) GO TO 17 

C FOLLOWING IS THE MESH SWEEP WITHIN ZONE IX 

DU 11 JX=1, METRO 
KAVG=1 
MESH=MESH+1 
14 RAVG=KA VG 

COVA=DIFFER*KAVG 

MM=MESH-KAVG 

LGRUP-0 
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J 1 6=0 

C SWEEP BY GROUP 

1 JIG=J1G+1 

C THE AO LOOP IS TO CHANGE TO MACROSCOPIC CROSS SCCTIONS FOR 

C COMPARISON TEST BY MESH INTERVAL WITHIN REGION, THEN BY ISOTOPE 

DO AO J I M=MM , MESH 
F AHK=0 • 0 

DO AO NI B= 1 1 NI SOT 
JID=IDVULC(NIB) 

IF( ABSIG( JID, JIG) .GT.O. ) GO TO 56 
ATDENt JIM,JID)=“AfOEN{ JIM, J ID) 

GO TO AO 

56 ATDENt JIM, JID)=ATOEN( J IM, JID)*ABSIG( JID, JIG) 

JF( JIM. NE. MESH) GO TO AO 

IF(JID.LE.12.0R.NKPT(NIB).NE.O> FAHR=FAHR+ATDEN( MESH, JID) 

AO CONTINUE 

FAHR=10.*FAHR*DIFFER 

C THE 12 LOOP IS THE COMPARISON TEST BY ISOTOPE FOR GROUP JIG 

3 DO 12 KQZ = 1 , N I SOT 
IDZ=IDVULC(KDZ ) 

IF (ATOENCMM, IDZ ) -GT.O. ) GO TO 15 

GUOI=1.0 

GO TO 16 

15 DUG I=ATDEN(MESH» IDZ)/ATDEN( MM , I DZ ) 

16 IF ( ABS tOUUl-1. ) .GT.COVA.AND.ATDENtMESH, 10Z ) .GT.FAHK) GO TO 19 

12 CONTINUE 
LGRUP=LGKUP+1 
GO TO 19 

A IF(KAVG.GT.l) GO TO 13 
C 

C JAVE AT ANY TIME IS EQUAL TO THE NUMBER OF ADJACENT COMPARISONS OR 

C MATCHES THAT HAVE ALREADY BEEN FOUND FOR ANY ONE AVERAGED SET. 

C 

JAVE=JAVE+1 

IF(JAVE.GT.l) GO TO 13 
C 

C KAKY IS THE MESH INTERVAL THAT CONTAINS THE AVERAGED ATOM DENSITY 

C (NEGATIVE) FOR THE NEXT MESHM1-KARY MESH INTERVALS. 

C 

K AK Y = ME SH- 1 

13 I F {MM.EO. KAKY ) GO TO IB 
C 

C KAVG IS THE ORDER OF THE COMPARISON AT HAND, E.G., =1 IF CUMPAR- 

C ISON IS WITH AN ADJACENT MESH INTERVAL, =2 IF THE COMPARISON IS 

C WITH THE MESH INTERVAL SECOND PRECEDING THE CURRENT ONE, ETC. 

C 

K A VG=K A VG+ 1 
GO TO 14 

IB IFUX.EO. METRO) GO TO 22 
GO TO 11 

C THE 46 LOOP CHANGES BACK TO ATOM DENSITIES FROM MACKOSCOPICS 

19 DO 46 J I M=MM, MESH 
DO 46 NOK = 1 , N I SOT 
J ID=I DVULC ( NOK ) 

IF ( ABS I G l JID, JIG) .GT.O.) GO TU 90 
ATDEN(JIM,JID)=— ATDEN(JIM,JID) 

GO TO 46 

90 ATDEN(JIM,JID)=ATDEN( J I M , J ID)/ABSIG(JID,J1G) 

46 CUNT INUE 

IFUIG.LT.NG) GO TO 1 
IF ( LGKUP. EU.NG) GO TO 4 
IF(KAVG.EUUJAVE) GO TO 26 
IF (JAVE.NE.U) GO TO 26 
K AR Y = ME SH— 1 
ME SHM1 = K AR Y 
IFUX.EO. METRO) MARS=1 
GO 10 21 

26 MESHM1=MESH— 1 
IFUX.EO. METRO) MAKS=1 
GO TO 21 

22 ME SHM1 = ME SH 
21 VOLUME=O.D 

DO 25 MEU=KARY,MESHMl 
25 VOLUME=VOLUME+VOL( MED) 

00 23 K DA= 1 , N I SOT 
I UA= I DVULC ( KDA ) 

AXE^O.O 

DU 24 MEU=KAKY,MESHM1 
24 AXE = AXE + ATDEN(MED, I DA ) *V0L ( MED ) 

23 ATDENtKARY, 1 DA ) =-A X E/ VOLUME 
11 CONTINUE 

IFIMARS.EO. 1 ) GO TO 17 
GO. TO 10 

17 DO 27 JMP = 1 , N I SOT 
I MP= I DVULC ( JMP) 

27 ATDEN1 MESH, I MP ) =-AT DEN < ME SH, I MP ) 

10 CONTINUE 
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C CALCULATIUN OF MACGG INPUT 

C 

WR I T£ ( 6 » 100 ) 

WRIT£(6, 112 ) 

NUDE=0 

MOOD=0 

DO 32 MOUT=l»Nrt 
DO 31 KOUT = 1 » N I SOT 
I OUT = I DVULC t KOUT ) 

IF ( ATDEN(MOUT» IOUT) .GT.O. ) GO TO 30 
ATDENtMOUT, I OUT ) =-ATDEN ( MOUT , I OUT ) 

31 CONTINUE 

I F l NUDE . 60.0 ) GO To 30 

33 NEND=MQUT- 1 
NAKT=MUUT-NUDE 
WKITE(6tl02) NUDEtNARTtNEND 
MOUD=MOOD+1 
NEMO(MOOD)=NUDE 

LUG=0 

C 

C REARRANGE ATOM DENSITIES IN MACGG ORDER 

C 

DO 36 MX= 1 » N I SOT 
I D= I DVULC ( MX ) 

CARRY t MX ) =0 . 0 

36 CARRY! MX )=ATDEN( NARTt ID) 

DO 35 LX=1 T NIS0T 
IFfCARRYlLX ) . LE.O. ) GO TO 35 
LUG=LUG+1 

CARRY ( LUG )=CARRY( LX) 

LABEL t LUG ) =LX 
35 CONTINUE 

WRITE(6t 106) 

WRITE (6 t 104) LUGt ( LABEL t K J ) t KJ=l ,LUG) 
PUNCH 1 05 t ( CARRY ( KM) ,KM=1 tUJG) 

WRITE (6t 107) I CARRY (KM) T KM=1 ,LUG) 

NUDE=1 
GO TO 34 
30 NUDE=NUDE+1 

34 IF ( MOUT.NE.NM ) GO TO 32 
M0UT=M0UT+1 

GO TO 33 

32 CONTINUE 
RETURN 
END 
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tIBFTC PLOTXY 

SUBROUTINE PLOTXY ( X* Y ,K ,P ) I 

C REFERENCE NASA TN D-2174 « APRIL»1964, BY LOIS T . DELLNER 

C AND BETTY JO MOORE 

C 2 

COMMON/ JOLO/N,F,DX,XYX,FORY,STUG,L ABOUT, T0NLY,KSW64,KPWR,KF0,TLINX 3 

LOGICAL X YX» FORY, stug,tonly,xgl,ls 4 

DIMENSION XI 1) ,Y( 1) ,Pll ) 5 

DIMENSION FLS< 3) ,FLAB(4) ,FYLAB<4) ,YlABEL I 11 ) »A( 104) ,ELS< 3) 6 

EQUIVALENCE I FLAB ( 3 ) , I FLAB3 ) , ( FYLAB (3) » I Yl AB ) 7 

DATA MASK1 , MASK2, MASK4, MASKS , MASK16,MASK32 ,MASK64 / 8 

1 01, 02, 04, 010, 020, 040 , 0100 / 9 

DATA FYLAB(l) »FYLAB< 2) , FYLAB3, FYLAB < 4 ) / 10 

1 6HI2HP , , 6H20 X ,11 ,6H F10.0,1H) / 11 

DATA BLANK, XGRID,YGRID /1H ,1H-,1H1 / 12 

DATA RMARK, PCSTD/0726060606060 , 1H*/ 13 

DATA FLSI 1 ) * FLS28, FLS264, FLS2B*FLS38, FLS3B / 14 

1 6HI2HP+»,4H11X,» 3H5X,,3HA5,,3HA6) ,6H2F6.3) / 15 

DATA FLAB(1),FUAB(2),FLAB3,FLAB(4) / 16 

1 6H(2HP+,,4H17X, ,6H F9.0,IH) / 17 

C 18 

100 WRITE (6,500) 19 

500 FORMAT 1 2HPT ) 20 

102 KODE=K 21 

N=P ( 1 ) 22 

LAB0UT=1 23 

FLAB ( 3) = FLAB3 24 

LS = .FALSE. 25 

FYLAB(3)=FYLAB3 26 

KSW0=O 27 

K$W64=0 28 

110 PC=PCSTD 29 

112 IF! (AND(K0DE»MASK1 > ) .GT.O. )PC - P t 2 ) 30 

114 M=10 31 

116 I F ( ( AND! KODE , MASK2) ).GT.0.)M=P(3) 32 

117 IF (M.EQ.O)M=1000 33 

118 NY=10 34 

120 IFt ( AND(K0DE,MASK4) ) .GT.O. > NY = P(4) 35 

121 IF(NY.EQ.O) NY=1000 36 

122 I F I lAND(K0DE,MASK64) ) .GT.O. )KSW6 4=2 37 

124 IF( (AND! KODE, MASK 8 > ) . GT .0. )KSW8 = 1 38 

C 39 

125 KB6 4=K5W8+K SW64 40 

126 I F ( K864-2 ) 132,128,138 41 

128 FLS ( 2 )— FLS264 42 

130 GO TO 139 43 

132 FLSI 2 ) =FLS28 44 

134 FLSl3)=FLS38 45 

136 GO TO 140 46 

138 FLS I 2 ) = FLS2B 47 

139 FLS I 3 )=FLS3B 48 

C 49 

140 XYX=. FALSE. 50 

142 FORY=. TRUE. 51 

144 STUG=. FALSE. 52 

146 TONLY=. FALSE. 53 

C 54 

148 IF! I ANO I KODE, MASK 32) .LE.O.) ) GO TO 172 55 

151 STUG=.TRUE. 56 

152 KSY=P I 9 ) 57 

154 PWRi0Y=10.**(KSY-6) 58 

156 FY =Pl 10)*PWR10Y 59 

158 F - FY 60 

C 61 

160 IFIPI5J.GE.2. ) GO TO 172 62 

162 TONLY=.TRUE. 63 

164 DY= Pill )*PWR10Y 64 

166 DX= DY 65 

172 CALL PISTUGIY) 66 

173 I F I DX.EQ.O. ) GO TO 700 67 

174 FY=F 68 

176 DY =0X 69 

180 IF I KSW64. EQ.2 ) KPWRY=KPWR 70 

190 IYLAB=IYLAB-KFD 71 

C 72 

200 XYX =.TRUE. 73 

202 FORY=. FALSE. 74 

204 $TUG=. FALSE. 75 

206 TONLY=. FALSE. 76 

208 TLINX=55=M l+N/35 > 77 

C 78 

210 IF! I ANDIK0DE,MASK16) .LE.O.) ) GO TO 232 79 

213 STUG=.TRUE. 80 

214 KSX = P ( 6 ) 81 

216 PWR10X=010.**(KSX-6) 82 

218 FX= PI7)*PWR10X 83 

220 F=FX 84 
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C 85 

222 IF(MODUFIX(Pt 5) ) 1 2).EQ.l) GO TO 232 86 

224 TONLY=.TRUE. 87 

226 DX =P(8)*PWR10X 88 

232 CALL PISTUG(X) 09 

IF(DX.EQ.O.) GO TO 700 90 

234 FX=F 91 

240 I F ( KSW64. EQ. 2 ) KPWRX=KPWR 92 

248 IFLAB3=IFLAB3— KFD 93 

C 94 

250 IF (KSW64.EQ.0JG0 TO 264 95 

252 KOUTX=—KP WR X 96 

254 KOUTY =— KP WRY 97 

256 F10X=10.**KPWRX 98 

258 F10Y=10.**KPWRY 99 

260 WRITE (6,502) KOUTX, KOUTY 100 

502 FORMAT ( 2HP T , 6X , 3HX*E,I2,4H Y*E,I2) 101 

C 102 

264 DO 278 1=1,11 103 

266 TEMP = FY+FLOAT( I- 1 ) #DY*10. 104 

268 ATEMP= ABS(TEMP) 105 

270 IF UTEMP.LT.l.E-7) TEMP = 0. 106 

272 IF < ATEMP.GE.1.E+7)LAB0UT=2 107 

278 YLABEL ( I )=TEMP 108 

300 KSYLAB =1 109 

302 WRITE ( 6 , FYLAB ) ( YLABEL ( I) , 1=1 , 1 1 ) 110 

304 GO TO ( 306, 700) , KSYLAB 111 

C 112 

306 KSYLAB =2 113 

310 LCTR=0 114 

NCTR=1 115 

KOUT=l 116 

KQU I T= 1 H7 

c 118 

320 1F(M0D(LCTR,M) >328,322,328 119 

322 AF I LL= XGRID 120 

324 XGL = . TRUE* 121 

GO TO 330 12 2 

328 XGL =. FALSE. 123 

AF I LL = BL ANK 124 

330 DO 332 1=2,104 125 

332 A ( I ) = AFILL 126 

334 DO 336 1=2, 104, NY 127 

336 A ( I ) = YGRID 128 

All) =BLANK 129 

338 GO TO ( 340,400 ),KOUT 130 

C 131 

340 KX = ( X ( NCTR ) — FX ) /DX +.5 132 

342 IF(KX-LCTR)630, 350,600 133 

350 KY= ( Y(NCTR)-FY)/DY+.5 134 

351 LS= .TRUE. 135 

352 TPC = PC 136 

353 KYL = KY+2 137 

354 I F ( KY.LT.O ) GO TO 360 130 

356 IF(KY.GT.lOl) GO TO 364 139 

358 GO TO 370 140 

360 KYL= 1 141 

362 GO TO 366 142 

364 KYL= 1 04 143 

366 TPC=RMARK 144 

C 145 

370 A ( KYL ) =TPC 146 

372 J=1 147 

374 IFIKSW8.EO.O) GO TO 300 148 

376 ELS(J)=P( NCTR+ 11) 149 

378 J= J + l 150 

380 IF(KSW64.E0.0) GO TO 386 151 

382 ELS( J)=X(NCTR)/F10X 152 

384 ELS ( J+l ) = Y ( NCTR ) /F10Y 153 

C 154 

386 IF(NCTR.GE.N)GO TO 392 155 

388 NCTR=NCTR+1 156 

390 GO TO 340 157 

C 158 

392 KOUT = 2 1 59 

394 M= 10 160 
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c 161 

400 IFIXGL.AND. (MODI LCTR , 10 ) . EQ.O >) KQUI T=2 162 

C 163 

600 WRITE (6,504M Milt 1 = 1,104) 164 

504 FORMAT ( 2HP ,26X,104A1) 165 

602 IF(K864.NE.0) GO TO 620 166 

C 167 

604 IF(M00(LCTR,10)>614»606,614 168 

606 XLABEL =FX+FL0AT ( LCTR ) *DX 169 

608 TEMP =ABS ( XLABEL I 170 

610 IF(TEMP.GE.1.E+7)LAB0UT=2 171 

612 IF(TEMP.LT.1.E-7)XLABEL=0. 172 

613 WRITE(6,FLAB> XLABEL 173 

C 174 

614 LCTR = LCTR+1 175 

616 GO TO (320,302 ) » KQU I T 176 

C 177 

620 IF(.NOT.LS) GO TO 604 178 

622 LS = .FALSE. 179 

624 WRITE 16»FLS) IELSI I>t I=1,K864> 180 

626 GO TO 604 181 

630 LABOUT= 4 182 

700 GO TO ( 710,702,706,704) ,LABOUT 183 

702 WRITE (6,506) 184 

506 F0RMAT(2HPL,3X,10HBA0 LABELS ) 185 

GO TO 720 186 

C 187 

704 WR I TE ( 6 , 520 ) NCTR 188 

520 FORMAT ( 2HPL, 18HX OUT OF DRDER . I=,I5> 189 

706 WRITE (6,508 H X(I ) ,Y( I ) , 1=1,2) ,K,P< 1) 190 

508 F0RMAT1 2HPL»5H N. G. , 4G20 .8 , 1 6 , F8 .2 > 191 

708 GO TO 720 192 

710 WRITE(6,510) 193 

510 FORMAT ( 2HPL ) 194 

720 RETURN 195 

END 196 



SIBFTC PISTUG ALTIO 

SUBROUTINE P I STUG ( ARRAY I 1 

DIMENSION ARRAY! 1) 2 

COMMON/ JOLO/N, F, DX, XYX, FORY, STUG, L ABOUT, T0NLY»KSW64»KPWR, KFD, TL I NX 3 

LOGICAL XYX»FORY, STUG * TONL Y A 

126 XI = ARRAY ( 1) 5 

128 IFIXYX) GO TO 133 6 

130 DO 132 J = 2»N 7 

132 XI = AMIN1(X1»ARRAY( J) > 8 

133 I F ( STUG ) X1=F 9 

134 XN = 0.0 10 

136 DO 146 J = 1 » N 11 

138 DIF = ABS(X1“ARRAY( J) ) 12 

140 IFIDIF.LE.XNIGO TO 146 13 

142 XN=DI F 14 

144 IHOLD= J 15 

146 CONTINUE 16 

147 XN = ARRAY! IHOLD) 17 

148 I F ( KSW64* EQ • 2 ) KPWR = KHAR ( AMAX1 ( ABS ( XI > , ABS ( XN > ) ) 18 

149 IF(TONLY) GO TO 240 19 

150 TL I N= 1 01. 20 

152 IF!. NOT. FORY) TLIN = TL I NX 21 

154 C5 = ! XN-X1 ) /TL 1 N 22 

156 C6 = ABS ! C5 ) 23 

158 IF(C6.EQ.O.) GO TO 300 24 

159 K7 = KHAR ( C6 ) 25 

160 C8 = 10.**K7 26 

162 C9 = C6/C8 27 

164 IF! !2.5-C9).LE.0.0) GO TO 172 28 

166 0=2. 29 

168 IF! (2.0-C9) .LE.O. I D=2.5 30 

170 GO TO 176 31 

172 D=5 • 32 

174 IF! !5.-C9).LE.0.0> D=10. 33 

176 Cll = D*C8 34 

178 DX = S I GN ( C 1 1 » C5 ) 35 

179 HUND = 100. *DX 36 

240 K7 = KHAR ( ABS ( DX ) ) 37 

250 KFD =0 38 

252 I F ! K7 ) 260,270,254 39 

254 IF(K7.GE.5) LAB0UT=2 40 

256 GO TO 270 41 

260 KFD =6 42 

262 IF(K7.LT. ( -7 > ) LABOUT = 2 43 

264 IF!K7.GT. (-6) ) KFD = -K7 44 

270 IF! STUG ) GO TO 230 45 

182 KC 1 2 = INT!ABS(X1 > / C 1 1 ) 46 

184 JJ = 1 47 

186 1FIXU 168,192,190 48 

188 JJ = 3 49 

190 IFIDX.LT.O.) JJ = JJ+1 50 

192 GO TO ! 193, 194, 195, 196) ,JJ 51 

193 KC 1 4 = KC 1 2 52 

GO TO 198 53 

194 KC 14 = KC12+1 54 

GO TO 198 55 

195 KC 14 = — KC 12—1 56 

GO TO 198 57 

196 KC 14 = — KC l 2 58 

198 KC 1 3 = MOO ! KC 1 2 , 10 ) 59 

KC 1 5 = KC12-KC13 60 

199 KC1 8 = KC 1 5 61 

200 GO TO 1212,202,202,210) »JJ 62 

202 KC 1 8 = KC18+10 63 

204 IFIKC13.NE.9) GO TO 210 64 

206 KC 1 8 = KC 1 4 65 

208 GO TO 212 66 

210 IFUl.LT.O. )KC18 = -KC18 67 

2 12 F=C11*FL0AT!KC18) 60 

214 IF!. NOT . FOR Y ) GO TO 230 69 

220 TEMP = F+HUND 70 

222 GO TO (224,228, 224, 228), JJ 71 

224 IFITEMP.GE.XN) GO TO 230 72 

226 GO TO 229 73 

228 IF l TEMP.LE.XN) GO TO 230 74 

229 F=C11*FL0AT ( KC 14) 75 

230 CONTINUE 76 

RETURN 77 

300 DX=0. 78 

LAB0UT=3 79 

GO TO 230 80 

END 81 

MBFTC KHAR ALTIO 

FUNCTION KHAR(XMAX) 1 

KHAR = INTI ALOGt XMAXI/2. 302585+40. 0)-40 2 

RETURN 3 

END 4 
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